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ABSTRACT 

Quasi-periodic brightness oscillations (QPOs) with frequencies ranging from ~ 300 Hz 
to ~ 1200 Hz have been discovered in the X-ray emission from fourteen neutron stars in 
low-mass binary systems and from another neutron star in the direction of the Galactic 
center. These kilohertz QPOs are very strong, with rms relative amplitudes ranging up 
to ~ 15% of the total X-ray countrate, and are remarkably coherent, with frequency to 
FWHM ratios as large as ~ 200. Two simultaneous kilohertz QPOs differing in frequency 
by ~ 250-350 Hz have been detected in twelve of the fifteen sources. 

Here we propose a model for these QPOs. In this model the X-ray source is a neutron 
star with a surface magnetic field ~ 10^-10^° G and a spin frequency of a few hundred 
Hertz, accreting gas via a Keplerian disk. Some of the accreting gas is channeled by 
the stellar magnetic field but some remains in a Keplerian disk flow that penetrates to 
within a few kilometers of the stellar surface. The frequency of the higher-frequency 
QPO in a kilohertz QPO pair is the Keplerian frequency at a radius near the sonic point 
at the inner edge of the Keplerian flow whereas the frequency of the lower-frequency 
QPO is approximately the difference between the Keplerian frequency at a radius near 
the sonic point and the stellar spin frequency. The difference between the frequencies of 
the pair of QPOs is therefore close to (but not necessarily exactly equal to) the stellar 
spin frequency. The amplitudes of the QPOs at the sonic-point Keplerian frequency and 
at the beat frequency depend on the strength of the neutron star's magnetic field and 
the accretion rate and hence one or both of these QPOs may sometimes be undetectable. 
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Oscillations at the stellar spin frequency and its overtones are expected to be weak but 
may sometimes be detectable. 

This model is consistent with the magnetic field strengths, accretion rates, and scat- 
tering optical depths inferred from previous modeling of the X-ray spectra and rapid 
X-ray variability of the atoll and Z sources. It explains naturally the frequencies of the 
kilohertz QPOs and the similarity of these frequencies in sources with different accretion 
rates and magnetic fields. The model also explains the high coherence and large am- 
phtudes of the kilohertz QPOs, and the steep increase of QPO amphtude with photon 
energy. The increase in QPO frequency with inferred accretion rate seen in many sources 
is also understandable in this model. 

We show that if the frequency of the higher-frequency QPO in a pair is an orbital 
frequency, as in the sonic-point model, the frequencies of these QPOs place interesting 
upper bounds on the masses and radii of the neutron stars in the kilohertz QPO sources 
and provide new constraints on the equation of state of matter at high densities. Further 
observations of these QPOs may provide compelling evidence for the existence of a 
marginally stable orbit, confirming a key prediction of general relativity in the strong- 
field regime. 

Subject headings: accretion, accretion disks — stars: neutron — stars: oscillations — 
stars: rotation — X-rays: stars 
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1. INTRODUCTION 

Observations of accreting neutron stars in low- 
mass X-ray binaries (LMXBs) with the Rossi X-ray 
Timing Explorer have revealed that the persistent X- 
ray emission of at least fifteen show remarkably co- 
herent quasi-periodic brightness oscillations (QPOs), 
with frequencies z^qpo ranging from ~ 300 Hz to 
~ 1200 Hz. These kilohertz QPOs are the highest- 
frequency oscillations ever seen in any astrophysical 
object. 

Eight of the fourteen identified sources in which 
kilohertz QPOs have been detected are "atoll" sources 
(4U 0614-F091, 4U 1608-52, 4U 1636-536, 4U 1728-34, 
KS 1731-260, 4U 1735-444, 4U 1820-30, and Aql X-1). 
Kilohertz QPOs have also been detected in all six of 
the originally identified "Z" sources (Sco X-1, GX 5—1, 
GX 17+2, GX 340+0, GX 349+2, and Cyg X-2). (For 
the definitions of atoll and Z sources, see Hasinger & 
van der Klis 1989.) Highly coherent brightness os- 
cillations with frequencies ranging from 360 Hz to 
~ 580 Hz have been detected during type I (thermonu- 
clear) X-ray bursts from four kilohertz QPO sources. 
Another neutron star in the direction of the Galac- 
tic center shows burst oscillations with a frequency of 
589 Hz, but this source has not yet been positively 
identified. The frequency ranges and rms amplitudes 
of the currently known kilohertz QPOs and burst os- 
cillations are listed in Table 1, with references. 

Two simultaneous kilohertz QPOs have so far been 
seen in six of the eight atoll sources in which kilohertz 
QPOs have been detected (all except 4U 1735-444 
and Aql X-1) and in all six of the Z sources. The dif- 
ferences Ai' between the frequencies of the two QPOs 
seen in these sources all fall in the range ~ 250- 
350 Hz. In the atoll sources 4U 0614+091 (Ford 
et al. 1996, 1997a), 4U 1636-536 (Wijnands ct al. 
1997b), and 4U 1728-34 (Strohmayer et al. 1996a, 
1996b, 1996c) and in the Z sources GX 5—1 (van der 
Klis et al. 1996e), GX 17+2 (Wijnands et al. 1997c), 
and Cyg X-2 (Wijnands et al. 1998) Av is constant in 
time, within the errors. In the Z source Sco X-1, Av 
changed from ^ 250 Hz when the accretion rate was 
probably near and sometimes even slightly greater 
than the Eddington critical rate to ^ 300 Hz when the 
accretion rate was significantly lower (van der Klis et 
al. 1997b). 

In 4U 0614+091, a brightness oscillation with a fre- 
quency consistent with the difference Av = 328 Hz 
between the frequencies of the two kilohertz QPOs 



Hcxai simultaneously in this source was detected with 
marginal significance during one 30-min interval (Ford 
et al. 1996, 1997a). In 4U 1728-34, a strong, rela- 
tively coherent brightness oscillation was detected in 
six of the twelve X-ray bursts so far observed from this 
source, with a frequency that has remained constant, 
within the errors, for more than a year and is con- 
sistent with the difference Aj^ = 363 Hz between the 
frequencies of the two kilohertz QPOs seen simultane- 
ously in this source (Strohmayer et al. 1996b, 1996c; 
Strohmayer, Zhang, & Swank 1997; Strohmayer 1997, 
personal communication). Burst oscillations have also 
been seen in 4U 1636—536, with a frequency ^580 Hz 
(Zhang et al. 1996, 1997), and in KS 1731-260, with 
a frequency ~520 Hz (Morgan & Smith 1996; Smith, 
Morgan, & Bradt 1997). In both cases the frequency 
of the burst oscillation is approximately twice the dif- 
ference between the frequencies of the two kilohertz 
QPOs observed simultaneously in these sources. 

The peaks in power density spectra associated 
with these kilohertz QPOs are relatively narrow, with 
FWHM SvQPo as small as 0.005 z/qpo- When in- 
tegrated over all photon energies, the rms relative 
amplitudes of the kilohertz QPOs seen in the atoll 
sources range up to ~ 15% of the total 2-60 keV 
X-ray countrate and are typically much greater than 
the < 1% amplitudes of the kilohertz QPOs ob- 
served in the Z sources. The amplitudes of the kilo- 
hertz QPOs increase steeply with photon energy, up 
to at least 15 keV, in all the atoll and Z sources 
where the photon-energy dependence of the kilohertz 
QPOs has been studied (4U 1608-52: Berger et al. 
1996; 4U 1636-536: Zhang et al. 1996; 4U 1728-34: 
Strohmayer et al. 1996c; KS 1731-260: Wijnands & 
van der Klis 1997; GX 5-1: van der Khs et al. 1996e; 
GX 17+2: Wijnands et al. 1997c; Cyg X-2: Wijnands 
et al. 1998). 

The frequencies of the kilohertz QPOs have been 
seen to vary by as much as a factor ~ 2. The frequen- 
cies of the two kilohertz QPOs observed in the persis- 
tent X-ray emission between type I X-ray bursts of the 
atoll sources 4U 0614+091 (Ford et al. 1996, 1997a), 
4U 1728-34 (Strohmayer ct al. 1996a, 1996b, 1996c), 
KS 1731-260 (Wijnands & van der Klis 1997), and 
4U 1820-30 (Smaie et al. 1996, 1997) increase steeply 
with increasing countrate (d log z//(i log CR > 1). 
The frequencies of the pair of kilohertz QPOs seen 
during two different observations of 4U 0614+091 
(Ford et al. 1997a) correlated tightly and positively 
with the energy of the peak in the keV X-ray spectrum 
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(Ford et al. 1997b), even though the tracks in the 
frequency- countrate plane made by the source dur- 
ing the two observations differed significantly. The 
frequencies of the pair of high-frequency QPOs ob- 
served in the Z sources Sco X-1, GX 5-1, GX 17-^2, 
and Cyg X-2 increase steeply with the inferred accre- 
tion rate in these sources (van der Khs et al. 1996a, 
1996b, 1996d, 1996e, 1997a; Wijnands et al. 1997c, 
1998; as discussed by van der Klis [1989] and Lamb 
[1989, 1991], there is no simple relation between coun- 
trate and accretion rate in the Z sources). 

Two simultaneous kilohertz QPOs are so common 
that attention has focused on models that might be 
able to explain the occurrence of such kilohertz QPO 
pairs. Strohmayer et al. (1996c) and Ford et al. 
(1997a) have suggested adapting the magnetospheric 
beat-frequency model to explain these QPOs. This 
model was originally proposed by Alpar & Shaham 
(1985) and Lamb et al. (1985) and developed further 
by Shibazaki & Lamb (1987) to explain the single 15- 
60 Hz horizontal-branch oscillation (HBO) discovered 
earlier in the Z sources (see van der Klis 1989). As 
we discuss in § 2, it is very difficult to understand 
several key features of the kilohertz QPOs, including 
the simultaneous occurrence of two kilohertz QPOs, 
in terms of the magnetospheric beat-frequency model. 

Here we propose a model for the kilohertz QPOs 
that is fundamentally different from the magneto- 
spheric beat-frequency model. In the model presented 
here, the X-ray source is a neutron star with a surface 
magnetic field ~ 10^^-10^° G and a spin frequency of a 
few hundred Hertz, accreting gas via a Keplerian disk. 
In this model the key characteristics of the kilohertz 
QPOs are explained as follows (details are presented 
in § 3): 

Accretion flow. — Magnetoturbulence, differential 
cooling, and radiation forces create density inhomo- 
gcneities ( "chmips" ) in the gas in the nearly Keplerian 
disk flow. These clumps are dissipated by the shear 
in the azimuthal velocity combined with gas pressure 
forces and turbulence. If the neutron star has a mag- 
netic field > 10'' G, some of the accreting gas is chan- 
neled out of the disk flow by the stellar magnetic field. 
Channeling occurs at several stellar radii, if the stellar 
magnetic field is ~ 10^-10^° G, but only very close to 
the stellar surface, if the field is ~ 10^-10^ G. Even 
if the field is as strong as ^ 10^ 10^^° G, some of the 
accreting gas continues to drift inward in a nearly Ke- 
plerian disk flow until either radiation drag forces be- 
come important, the gas reaches the innermost stable 



circular orbit, or the gas begins to interact viscously 
with the stellar surface. The fraction of the accreting 
gas that remains in the Keplerian flow close to the 
star depends primarily on the strength of the stellar 
magnetic field and the accretion rate. 

Role of radiation drag and general relativity. — For 
the accretion rates typical of the Z and atoll sources, 
the nearly Keplerian flow ends and the inward radial 
velocity of the gas in the disk increases sharply as the 
gas nears the star, either because of the azimuthal 
drag force exerted on the gas in the disk by the radi- 
ation coming from near the stellar surface or — if the 
drag force is weak and the radius of the neutron star 
is less than the radius of the innermost stable circular 
orbit — because of the general relativistic corrections 
to Newtonian gravity that produce an innermost sta- 
ble orbit. In either case, the inward radial velocity 
of the gas becomes supersonic within a very short ra- 
dial distance from the radius at which it begins to 
increase. 

Effect of inhomogeneities in the flow. — Roughly 
speaking, clumps that form outside the radius where 
the radial inflow becomes supersonic are dissipated 
before gas from them reaches the stellar surface, 
whereas gas from clumps that form near the radius 
where the flow becomes supersonic falls inward and 
collides with the stellar surface before the clumps are 
destroyed. Where the streams of denser gas from 
the clumps orbiting near the sonic point collide with 
the stellar surface, the streams create arc-shaped ar- 
eas of brighter X-ray emission that move around the 
star's equator with a frequency equal to the orbital 
frequency of the clumps, which is approximately the 
Keplerian frequency at the sonic point. 

Higher-frequency kilohertz QPO.- -The frequency 
t'QP02 of the higher-frequency of the two kilohertz 
QPOs is approximately equal to the Keplerian fre- 
quency at the sonic point. A detectable quasi-periodic 
oscillation of the X-ray flux and energy spectrum may 
be produced at this frequency by the periodically 
changing aspect (and. for most viewing directions, by 
the neutron star's periodic eclipse) of the arc-shaped 
areas of brighter X-ray emission, as they move around 
the star with a frequency equal to the orbital fre- 
quency of clumps in the Keplerian disk flow near the 
sonic point. The oscillation appears quasi-periodic 
becaiise the chimps have finite lifetimes and slightly 
different orbital frequencies. We find that v/dv can 
be as large as ^ 100 for these oscillations, consistent 
with the observed coherence of the kilohertz QPOs. 
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The frequency of the sonic-point Kcplcrian QPO is 
almost independent of the spin rate of the neutron 
star (frame-dragging causes the Keplerian frequency 
at the sonic point, and hence the QPO frequency, to 
depend weakly on the star's spin rate). 

Lower-frequency kilohertz QPO. — The frequency 
j/QPOi of the lower-frequency of the two kilohertz 
QPOs is equal to one (or possibly two) times the dif- 
ference between the Keplerian frequency at a radius 
near the sonic point and the stellar spin frequency. A 
QPO may be produced at this beat frequency if the 
magnetic field of the neutron star is weak enough that 
a Keplerian disk flow penetrates close to the star but 
strong enough to channel some of the accreting gas 
before it collides with the stellar surface. For the ac- 
cretion rates of the Z and atoll sources, this requires 
surface magnetic fields Bs in the range ~ 10^-10^" G. 

Channeling of gas out of the disk generates a 
QPO at the beat frequency because it creates slightly 
brighter regions that rotate with the star (these slightly 
brighter regions are distinct from the brighter spots 
that move around the star at the sonic-point Keple- 
rian frequency). The slightly brighter regions that ro- 
tate with the star produce a weakly-beamed radiation 
pattern that rotates at the stellar spin frequency. The 
radiation in this pattern produces a slight increase of 
the azimuthal drag force on the gas in the clumps or- 
biting near the sonic point, once or twice each beat 
period (depending on the symmetry of the slightly 
brighter region). 

This oscillation of the drag force in turn creates 
an oscillation of the inward flux of gas from a given 
clump at the sonic point once or twice each beat pe- 
riod. This modulation of the mass flux in the gas 
streams flowing inward from all clumps orbiting near 
the sonic point generates a quasi-periodic oscillation 
of the luminosity and spectrum of the X-ray emis- 
sion from the stellar surface, with a frequency equal 
to one (or possibly two) times the difference between 
the Keplerian frequency at a radius near the sonic 
point and the stellar spin frequency. The difference 
between the frequencies of the two QPOs in a pair is 
therefore close to (but not necessarily exactly equal 
to) the stellar spin frequency. 

Oscillation amplitudes. — The amplitudes of the 
QPOs at the sonic-point Keplerian frequency and at 
the beat frequency depend on the strength of the neu- 
tron star's magnetic field and the mass accretion rate, 
and hence one or both of these QPOs may sometimes 
be undetectable. Oscillations at the stellar spin fre- 



quency and its overtones arc expected to be weak but 
may sometimes be detectable. 

Similarity of kilohertz QPO frequency behavior in 
different sources. — The similar frequency ranges and 
the similarity of the frequency-accretion rate corre- 
lations of the kilohertz QPOs observed in the Z and 
atoll sources — which have accretion rates and mag- 
netic fields that differ by factors of ten or more — is 
explained by the crucial role of radiation drag, which 
is only effective within one or two stellar radii, the ap- 
proximate proportionality of the angular momentum 
carried by the accreting matter and the radiation drag 
force, the fact that the stellar magnetic fields as well 
as the accretion rates are larger in the Z sources, and 
the tendency for the vertical thickness of the disk to 
increase with increasing mass flux. 

The radius where the radiation drag or general rel- 
ativistic corrections to Newtonian gravity cause the 
radial velocity to increase abruptly plays a key role in 
the model of the kilohertz QPOs proposed here. Once 
the radial velocity begins to increase, the flow be- 
comes supersonic within a very short radial distance 
(the inward radial velocity increases so abruptly that 
the location of the sonic point is insensitive to the 
precise deflnition used), so the sonic point is a use- 
ful indicator of the point in the disk flow where the 
radial velocity begins to increase steeply. Hence, for 
convenience we shall call the point where the radial 
velocity increases sharply the "sonic point" and refer 
to this model for the kilohertz QPOs as the "sonic- 
point model". 

We show that if the frequency of the higher- frequency 
QPO in a pair of kilohertz QPOs is the orbital fre- 
quency of gas in a stable Keplerian orbit around the 
neutron star, as in the sonic-point model proposed 
here, the frequencies of these QPOs provide interest- 
ing new upper bounds on the masses and radii of the 
neutron stars in the Z and atoll sources and new con- 
straints on the equation of state of matter at high 
densities. 

Any model of how the kilohertz QPOs are pro- 
duced must be consistent with the known properties 
of the atoll and Z sources. Hence, before analyzing 
the sonic-point QPO model of the kilohertz QPOs in 
more detail, we first summarize in § 2 the previously 
known X-ray spectral and lower-frequency X-ray vari- 
ability properties of the atoll and Z sources and the 
physical picture of these sources that has been devel- 
oped based on these properties. In § 3 we analyze the 
physics of the sonic-point model and show that it is 
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consistent with tlio basic properties of the kilohertz 
QPOs. In § 4 we demonstrate that the sonic-point 
model is also consistent with the existing physical pic- 
ture of the Z and atoll sources and with many of the 
more detailed properties of the kilohertz QPOs. In § 5 
we show how to derive upper bounds on the masses 
and radii of the neutron stars in the kilohertz QPO 
sources from the frequencies of stable Keplerian or- 
bits in the kilohertz range and discuss the constraints 
on the properties of neutron-star matter that follow 
from these bounds; the bounds we derive include the 
effects of frame-dragging. Finally, in § 6 we discuss 
several specific predictions of the sonic-point model. 

2. PROPERTIES OF THE NEUTRON STARS 
IN LMXBS 

2.1. Observed Properties of the Atoll and Z 
Sources 

The atoll sources are LMXBs that, over time, trace 
atoll-shaped patterns in X-ray color-color diagrams 
(Hasinger & van der Klis 1989). They have lumi- 
nosities L ~ lO^MO^"^ ergss-\ i.e., -1-10% of the 
Eddington critical luminosity Le oi & neutron star. 
Power-density spectra of their brightness variations 
show broad, band-limited noise components at fre- 
quencies below ~ 100 Hz. No QPOs with frequen- 
cies < ICQ Hz have so far been detected in any of the 
atoll sources with the exception of Cir X-1, which has 
a QPO that increases in frequency from 1 to 30 Hz 
as the countrate increases, when the source is very 
bright (Oosterbroek et al. 1995; Shirey et al. 1996; 
Bradt, Shirey, & Levine 1998). 

Recent comparisons of models of the X-ray emis- 
sion of neutron stars in LMXBs with the X-ray spec- 
tra of atoll sources observed with EXOSAT {Vsaiiis & 
Lamb 1998a, 1998b, 1998c) suggest that these sources 
can be subdivided into two groups, the "4U" atoll 
sources (such as 4U 1636-53, 4U 1705-44, 4U 1820-30, 
4U 1608-52, and 4U 1728-34) and the "GX" atoll 
sources (such as GX 9+1, GX 9+9, GX 3+1, and 
GX 13+1), based on the strengths of their inferred 
magnetic fields; one atoll source, 4U 1735—44, has in- 
termediate spectral properties and therefore probably 
has a magnetic field of intermediate strength. 

The Z sources are LMXBs that produce a charac- 
teristic Z-shapcd track in X-ray color-color diagrams 
(Hasinger & van der Klis 1989). They have luminosi- 
ties L ~ 10^* ergss"^, i.e., comparable to Le- The 
three branches of the Z are called the horizontal, nor- 



mal, and flaring branches. When a Z source is on 
the horizontal branch, a QPO with a frequency in the 
range 15-60 Hz is observed (van der Klis et al 1985; 
van der Klis 1989). This "horizontal branch oscilla- 
tion" (HBO) is also detectable in some Z sources when 
they are on the upper part of the normal branch. The 
relative width Su /u of the HBO peak in power den- 
sity spectra is typically ~ 0.1-0.3 and its centroid fre- 
quency typically increases with increasing countrate. 
As a Z source moves down the normal branch, the 
HBO becomes weaker and eventually disappears into 
the noise continuum. 

Near the middle of the normal branch, a dif- 
ferent QPO appears. The peaks in power density 
spectra produced by this QPO have relative widths 
5v/v ~ 0.3 and centroid frequencies in the range 4- 
8 Hz (Middleditch & Priedhorsky 1986; van der Klis 
1989). The properties of this QPO do not vary ap- 
preciably on the lower normal branch. As a Z source 
moves from the lower normal branch to the flaring 
branch, the frequency of this second QPO increases 
abruptly to ~ 15-20 Hz and the QPO becomes weaker 
and less coherent, eventually disappearing into the 
noise (van der Klis 1995; Dieters & van der Klis 1997). 
This QPO is the "normal/flaring branch oscillation" 
(N/FBO). 

In addition to these two types of quasi-periodic 
oscillations, power spectra of the brightness varia- 
tions of the Z sources also show three distinct band- 
limited noise components at frequencies below ~ 
200 Hz (Hasinger & van der Klis 1989). The proper- 
ties of the QPO and noise components vary system- 
atically with the position of a source on its Z track 
(Hasinger & van der Klis 1989; see Kuulkers 1995, Di- 
eters & van der Klis 1997, and Wijnands et al. 1997a 
for detailed studies of this behavior and possible ex- 
ceptions). 

Recent analyses of archival EXO SAT data (Kuulk- 
ers 1995; Kuulkers et al. 1995, 1996) suggest that the 
Z sources can be subdivided into two groups, based 
on the morphology of their Z tracks and power spec- 
tra: the "Cyg-like" Z sources (Cyg X-2, GX 5-1, and 
GX 340+0) and the "Sco-like" Z sources (Sco X-1, 
GX 17+2, and GX 349+2). At a more detailed level, 
GX 17+2 shares some of the characteristics of the 
"Cyg-like" sources and therefore probably has inter- 
mediate properties (Hasinger & van der Klis 1989; 
Wijnands et al. 1997a). 
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2.2. Current Physical Picture of the Atoll 
and Z Sources 

The atoll and Z sources are neutron stars accret- 
ing gas from a Keplerian disk fed by a low-mass 
companion star (see van der Klis 1989). The mag- 
netic fields and the accretion rates of these neutron 
stars are thought to be the most important parame- 
ters that determine their X-ray spectral and tempo- 
ral characteristics. According to the most complete 
and self-consistent current model of these sources, 
the so-called unified model (Lamb 1989, 1991), the 
atoll sources have dipole magnetic fields < 5 x 10^ G 
and luminosities ~ 1-10% of the Eddington critical 
luminosity, whereas the Z sources have dipole fields 
~ 10^-10"'^° G and luminosities very close to (and 
sometimes slightly above) the Eddington critical lu- 
minosity. 

X-ray spectra. — The X-ray spectra of the Z and 
atoll sources and the low upper limits on the am- 
plitudes of any periodic variations of their persistent 
X-ray brightness (less than 1% in some cases; see, 
e.g., Vaughan et al. 1994) constrain the properties of 
these neutron stars. The low upper limits on bright- 
ness variations at their spin frequencies can be under- 
stood if the magnetic fields of these neutron stars are 
< 10^" G and the magnetospheres and inner disks are 
surrounded by a central coronae with electron scat- 
tering optical depths ~ 3-5, even at low accretion 
rates (Lamb et al. 1985; Lamb 1989, 1991), because 
such coronae strongly suppress the X-ray flux oscilla- 
tion produced by any radiation pattern that rotates 
with the star (Brainerd & Lamb 1987). This effect is 
discussed further in § 3.6. 

When, as in the Z sources, the total luminosity of 
the neutron star, inner disk, and central corona be- 
comes comparable to Le, the vertical radiation force 
drives gas upward, out of the disk and radiation drag 
removes its angular momentum in less than one or- 
bit, creating a region of approximately radial inflow 
that extends out to ~ 300 km (Lamb 1989, 1991). 
Oscillations in this radial flow are thought to be re- 
sponsible for the N/FBO (Lamb 1989; Fortner et al. 
1989; Fortner 1992; Miller & Lamb 1992). 

Detailed physical modeling of the X-ray spectra 
of the atoll and Z sources indicates that soft 
0.5-1 keV) photons are produced by optically-thick 
bremsstrahlung and other processes at the surface 
of the star and by self-absorbed, high-harmonic cy- 
clotron emission in the inner magnetosphere (Psaltis, 



Lamb, & Miller 1995; Psaltis & Lamb 1998a, 1998b, 
1998c). As we discuss below, the accretion rate and 
the strength of the stellar magnetic field determine 
whether emission from the stellar surface or cyclotron 
emission in the magnetosphere is the dominant source 
of soft photons. These soft photons are upscattered 
by electrons in the magnetosphere and the central 
corona that surrounds the neutron star to produce 
the X-ray spectrum that emerges from the corona. In 
the Z sources, scattering by cooler electrons in the 
region of approximately radial inflow that surrounds 
the inner disk and corona further deforms the X-ray 
spectrum. 

Cyclotron emission is the dominant source of soft 
photons if upscattered cyclotron photons are able to 
supply the full accretion luminosity (see Psaltis et al. 
1995). The range of stellar magnetic flelds for which 
cyclotron emission is dominant can be estimated as 
follows. Within the magnetosphere, the spectrum of 
self-absorbed cyclotron photons can be approximated 
by a blackbody spectrum truncated at the energy 
at which cyclotron emission becomes optically thin. 
Comptonization within the magnetosphere and cen- 
tral corona increases the energy of a typical photon by 
a factor ~ , where y = {AkBTe/rrieC^) maxfr, t^} is 
the Compton y parameter, and m^. are the electron 
temperature and rest mass, r is the electron scatter- 
ing optical depth, and fc^ is the Boltzmann constant 
(see Rybicki & Lightman 1979, pp 195-223). Hence 
cyclotron emission is the dominant source of soft pho- 
tons if e^Lcyc ~ GMM/R or, equivalently, 

Vl06cmy V15/ \1AMqJ VlO^cmy ' 

(1) 

where M^s is the mass accretion rate onto the neu- 
tron star surface, i?cyc is the effective radius of the 
cyclotron photosphere, n is the harmonic number at 
which the transition from optically thick to optically 
thin emission occurs, and R is the radius of the neu- 
tron star. Figure || shows how fj.cyc,27 depends on M 
for the range of electron temperatures expected in the 
magnetospheres and central coronae of the atoll and 
Z sources. 

As Figure ^ shows, if the star's magnetic field is 
< 5 X 10^ G and the accretion rate is ~ 0.01-0.03 Me, 
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the dominant source of photons is the thermal emis- 
sion from the surface of the neutron star. These pho- 
tons are then upscattered by the electrons in the mag- 
netosphere and the central corona. Numerical calcu- 
lations of the X-ray spectra produced by stars with 
these field strengths and accretion rates agree well 
with iJXOS'^r observations of the spectra of the "4U" 
atoll sources when they are in the so-called "banana" 
spectral state (Psaltis & Lamb 1998c). 

If instead the star's magnetic field is ~ 5 x 10^-5 x 
10^ G and the accretion rate is <0.1Me, cyclotron 
emission is the dominant source of photons. These 
photons are then upscattered by the electrons in the 
magnetosphere and the central corona. Numerical 
calculations of the X-ray spectra produced by stars 
with these field strengths and accretion rates agree 
well with iJXOS'^ T observations of the spectra of the 
"GX" atoll sources (Psaltis & Lamb 1998a, 1998b, 
1998c). 

If the stellar magnetic field is ~ 10^-10^^° G and 
the mass accretion rate is > 0.5 Me, electron cy- 
clotron emission in the magnetosphere is very efficient 
in producing soft photons. At these accretion rates 
soft photons are Comptonized not only by the hot 
electrons in the magnetosphere and central corona but 
also by the cool electrons in the approximately radial 
inflow. Numerical computations of the X-ray spec- 
tra and color tracks predicted by this physical model 
agree well with EXOSAT and Ginga measurements 
of the X-ray spectra and color tracks of the Z sources 
(Psaltis et al. 1995; Psaltis & Lamb 1998a, 1998b, 
1998c). 

Accretion flows. — If the magnetic field of the neu- 
tron star is < 10^ G and the accretion rate is > 0.01 Le, 
the Keplerian disk flow may extend inward all the way 
to the surface of the neutron star without being chan- 
neled by the stellar magnetic field. (Here and below 
we refer to accretion flows confined near the orbital 
plane as disk flows, whether or not the azimuthal ve- 
locity field is nearly Keplerian, and as Keplerian disk 
flows if the azimuthal velocity field is nearly Keple- 
rian, that is, if the flow is nearly circular.) If instead 
the stellar magnetic field is > 10^^ G, most of the 
gas in the disk will couple to the stellar field and be 
funneled out of the disk plane toward the magnetic 
poles of the star at a characteristic cylindrical radius 
vuo (Ghosh & Lamb 1979a, 1979b, 1992) that is much 
larger than the radius of the star. Hence, the disk flow 
around such a star ends far above the stellar surface. 

For stellar magnetic fields of intermediate strength. 



some of the accreting gas is likely to couple to the 
stellar magnetic field beginning at wq, which will then 
channel it toward the star's magnetic poles, but some 
is also likely to continue as a disk flow inside Wq, as 
a result of Rayleigh- Taylor instability and incomplete 
coupling of the flow in the inner disk to the stellar 
magnetic field (see, e.g., Scharlemann 1978; Ghosh & 
Lamb 1979a, 1991; Lamb 1984; Spruit & Taam 1990). 
In this case we expect some of the gas in the disk flow 
to be channeled out of the disk by the magnetic field 
over a range of radii inside wq, with the rest of the 
gas remaining in a geometrically thin Keplerian flow 
that penetrates close to the stellar surface, as shown 
in Figure [|. This last case is the one that we expect 
to be relevant to the Z sources and most of the atoll 
sources. 

Ghosh & Lamb (1992) derived expressions for the 
characteristic cylindrical radius zuq at which gas in 
the disk couples strongly to the stellar magnetic field 
and begins to be channeled out of the disk flow, as 
a function of the stellar dipole magnetic moment fi, 
the mass accretion rate Mi through the inner disk, 
and the neutron star gravitational mass M, both 
for gas-pressure-dominated (GPD) and for radiation- 
pressure-dominated (RPD) inner disks, assuming zuq 
is much larger than the radius R of the star. The GPD 
model is expected to be valid for the atoll sources. 
The disks around the Z sources are expected to be 
RPD. The scaling of -cuo with /i. Mi, and M given 
by the analytical RPD model is expected to be valid 
for the Z sources, but the vertical thickness of the 
boundary layer given by the model is not expected to 
be reliable. We have therefore reduced the vertical 
thickness of the boundary layer in the model so that 
the Keplerian frequency at the radius where the gas 
couples strongly to the magnetic field agrees quanti- 
tatively with the Keplerian frequency inferred from 
spin frequency of Cyg X-2, based on the frequencies 
of its kilohertz QPOs (Wijnands et al. 1998; Psaltis 
et al. 1998) and the magnetospheric beat-frequency 
interpretation of its HBO (Psaltis et al. 1998), for a 
magnetic moment equal to the moment inferred from 
its X-ray spectrum (Psaltis & Lamb 1998c). Specif- 
ically, the Keplerian frequency at the coupling ra- 
dius was assumed to be 380 Hz, as indicated by the 
sum of the spin and HBO frequencies when Cyg X-2 
is on the normal branch, for a magnetic moment of 
3.5 X 10^^ G cm'^ and a mass flux through the inner 
disk equal to the mass accretion rate Me onto the 
neutron star that produces the Eddington critical lu- 
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minosity. 

The resulting two expressions for n7o(/i, Mi, M) can 
be solved for the dipole moment /io that gives a Kep- 
lerian frequency i^ko at the coupling radius; this char- 
acteristic dipole moment is 



1^0,27 



1.5 (i/Ko/llOO Hz 
X (M/1.4 Mq^-^^ 



0.43 



for GPD disks 



0.88 (i^co/llOO Hz)~^-^ {ivh/Mi 



X (M/1.4 Mq 



,0.91 



0.29 



for RPD disks 



(2) 

Here /Lio,27 = ^o/lO^^ Gcm'^ . If the star's dipole mo- 
ment exceeds /ig, most of the gas in the disk flow will 
couple to the magnetic field and be channeled out of 
the disk at a radius where the Keplerian frequency is 
less than z^ko- Figure ^ shows the dipole moment that 
gives Keplerian frequencies of 500 and 1100 Hz at the 
main coupling radius, as a function of the mass accre- 
tion rate. These two curves show that for the neutron 
star magnetic fields inferred from X-ray spectral mod- 
eling, the Keplerian frequency at the gas coupling ra- 
dius cannot be made simultaneously consistent with 
the ~500-1100 Hz ranges of the kilohertz QPOs in 
the 4U and Z sources. 

Accreting gas that is not channeled out of the disk 
plane at voo is expected to continue to drift slowly 
inward in a Keplerian disk flow until (a) it is chan- 
neled out of the disk at smaller radii, (b) it begins 
to flow rapidly inward, either because it loses angular 
momentum to the radiation drag force or because it 
crosses the radius of the marginally stable orbit, or 
(c) it collides with the stellar surface. Gas that is 
flowing rapidly inward remains in a (non-Keplerian) 
disk flow unless it is channeled out of the disk by the 
stellar magnetic field. We emphasize that the accre- 
tion disk ends above the stellar surface only if all the 
gas in the disk is channeled out of the disk flow before 
it reaches the star. 

In the magnetospheric beat-frequency model of 
the 15-60 Hz HBO observed in the Z sources (Al- 
par & Shaham 1985; Lamb et al. 1985; Shibazaki & 
Lamb 1987), a substantial fraction of the gas in the 
Keplerian disk flow is channeled out of the disk at 
^ 2-3 stellar radii and flows toward the star's mag- 
netic poles. This model therefore requires that the Z 
sources have magnetic fields ^ 10^-10^° G. The atoll 
sources have not previously shown any direct indica- 
tions that they have dynamically important magnetic 



fields (e.g., no HBO-like QPOs have so far been ob- 
served in an atoll source), and hence their magnetic 
fields are thought to be < 5 x 10^ G. However, the 
existence of kilohertz QPOs at the difference between 
an orbital frequency and the spin frequency of the 
neutron star indicates that most if not all of the atoll 
sources have magnetic fields > 10^ G. This is consis- 
tent with the detailed modeling of their X-ray spec- 
tra described above, which indicates that some atoll 
sources have magnetic fields as strong as 5 x 10® G. 
These sources may produce weak HBOs, which may 
be detected in the future. 

2.3. Magnetospheric Beat-Frequency Model 
and the Kilohertz QPOs 

As noted in § 1, Strohmayer et al. (1996c) and Ford 
et al. (1997a) have suggested adapting the magneto- 
spheric beat-frequency QPO model to explain the two 
simultaneous kilohertz QPOs. As described in § 2.2, 
the magnetospheric beat-frequency model was origi- 
nally proposed by Alpar & Shaham (1985) and Lamb 
et al. (1985) and developed further by Shibazaki & 
Lamb (1987) to explain the single, variable frequency, 
15-60 Hz HBO seen in the Z sources. Strohmayer et 
al. and Ford et al. have suggested that the frequency 
of the higher-frequency QPO in a kilohertz QPO pair 
might be the Keplerian orbital frequency at the ra- 
dius where gas in the accretion disk begins to couple 
to the stellar magnetic field and that the frequency 
of the lower-frequency QPO is the beat between this 
frequency and the stellar spin frequency. There are 
several difficulties with this application of the magne- 
tospheric beat-frequency model. 

Widespread occurrence of two simultaneous, highly 
coherent kilohertz QPOs. The occurrence of two si- 
multaneous, strong (amplitudes up to ~ 15% rms), 
and highly coherent {y/5v ~ 50-100) kilohertz QPOs 
appears inconsistent with the magnetospheric beat 
frequency model, because a strong, highly coherent 
oscillation with a frequency equal to the orbital fre- 
quency where gas first couples to the stellar magnetic 
field is not expected, nor is it observed in other types 
of sources where the magnetospheric beat frequency 
mechanism is thought to be operating. 

The magnetospheric beat-frequency mechanism is 
expected to produce a single, strong, moderately co- 
herent oscillation with a frequency equal to the beat 
frequency (power is of course also expected at over- 
tones of this frequency), but it is not expected to 
produce a strong, highly coherent oscillation with a 
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frequency equal to the or})ital frc;quency at the initial 
gas coupUng radius, for two reasons. First, no effect is 
known that would select a sufficiently narrow range 
of radii near this radius. Second, even if a narrow 
range of radii were selected, no mechanism is known 
that would make the orbital frequencies at these radii 
visible. 

Gas is expected to leave the disk along field lines 
that thread a very narrow annulus where the gas cou- 
ples strongly to the magnetic field of the star (Ghosh 
& Lamb 1979a). The range of orbital frequencies in 
this annulus is likely to be very small, consistent with 
the ^ 5-10 Hz observed width of the HBO (Lamb et 
al. 1985). The magnetic field is expected to couple 
more weakly to the gas in the disk over a wider annu- 
lus (Ghosh & Lamb 1979a). However, in both annuli 
the gas-magnetic field interaction repeats each beat 
period, not each orbital period (Lamb et al. 1985), so 
there is no obvious way that this interaction would 
make the orbital period visible in X-rays. No other 
mechanism is known that would single out a narrow 
range of radii in the inner disk. 

Even if a narrow range of radii were singled out, 
little luminosity is expected to be generated in such 
a narrow annulus, as was explained in the context of 
the HBO by Lamb (1988). Generation of a QPO by 
beaming of radiation via periodic occultation of the 
emission from the neutron star by clumps orbiting at 
the initial coupling radius would require us to be view- 
ing all the kilohertz QPO sources from within a small 
range of special inclinations and would also require 
the height of the disk at its inner edge to vary with 
accretion rate in just the right way (again see Lamb 
1988). Moreover, such beaming would be strongly 
suppressed by scattering in the ionized gas that sur- 
rounds the neutron star (see Lamb et al. 1985; Brain- 
erd & Lamb 1987; Kylafis & Phinney 1989; Lamb 
1988, 1989; Psaltis et al. 1995). 

A QPO that appears to be the magnetospheric 
beat frequency QPO has been observed in several 
accretion-powered pulsars (sec, e.g., Angclini, Stella, 
& Parmar 1989; Lamb 1989; Shibazaki 1989; Shinoda 
et al. 1990; Lamb 1991; Finger, Wilson, k Harmon 
1996; Ghosh 1996). In these pulsars only a single, 
fairly coherent oscillation has been seen, with proper- 
ties which indicate that it is the magnetospheric beat- 
frequency QPO; no QPO is observed at what would 
be the Keplerian frequency at the coupling radius. 

To summarize this point, the theory of the magne- 
tospheric beat-frequency QPO mechanism predicts. 



and observations of accretion-powered pulsars con- 
firm, that this mechanism produces only a single 
strong, moderately coherent oscillation, with a fre- 
quency equal to the beat frequency. In order to ex- 
plain the two simultaneous kilohertz QPOs, the mag- 
netospheric beat-frequency mechanism would have to 
generate two strong, highly coherent oscillations, at 
the beat frequency and at the Keplerian frequency 
at the radius where gas in the disk first strongly cou- 
ples to the magnetic field. This is the most important 
reason why the magnetospheric beat-frequency mech- 
anism appears unpromising as an explanation of the 
two kilohertz QPOs. 

Simultaneous occurrence of two kilohertz QPOs 
and the HBO in the Z sources. Identification of the 
15-60 Hz HBO in the Z sources as the magneto- 
spheric beat-frequency oscillation is supported by the 
fact that (1) the magnetic field strengths required are 
in agreement with those inferred from the spin-down 
rates of the millisecond recycled rotation-powered 
pulsars and detailed modeling of the X-ray spectra 
of the Z sources (see Alpar & Shaham 1985; Lamb 
et al. 1985; Ghosh & Lamb 1992; and § 2.2), (2) the 
~ 200-350 Hz neutron star spin frequencies predicted 
by the magnetospheric beat-frequency model (Lamb 
et al. 1985; Ghosh & Lamb 1992) are consistent with 
the spin rates of the recycled rotation-powered pul- 
sars (Ghosh & Lamb 1992) and the ~ 250-350 Hz 
spin frequencies inferred from the difference between 
the frequencies of the two kilohertz QPOs seen in the 
Z sources (see § 1), and (3) only the fundamental and 
(sometimes) the first overtone of the HBO have ever 
been detected in the Z sources when they are on the 
horizontal branch, as one would expect if the HBO is 
the magnetospheric beat-frequency QPO (see van der 
Klis 1989). If the HBO is the magnetospheric beat- 
frequency oscillation in the Z sources, then neither of 
the two kilohertz QPOs can be the magnetospheric 
beat-frequency oscillation, because the HBO was al- 
ways present when the kilohertz QPOs were detected 
(van der Klis et al. 1996e). 

Steep increase of kilohertz QPO frequencies with 
increasing mass accretion rate. In at least some 
sources, the frequencies of both kilohertz QPOs in- 
crease more steeply with inferred mass flux through 
the inner disk than is expected in the magnetospheric 
beat-frequency model, which predicts vko oc M°'^ for 
a GPD disk and i^ko oc M°-2 for an RPD disk (Ghosh 
& Lamb 1992; see Figure 2). It is therefore difficult 
to explain why the frequencies of the kilohertz QPOs 



10 



increase so steeply with Mi using this model. 

The magnetospheric beat-frequency model can ex- 
plain the steep increase of the HBO frequency with 
Mi observed in the Z sources, because these neutron 
stars are likely to have been spun up to frequencies 
comparable to the Keplerian frequency at the gas- 
magnetic field coupling radius, so the HBO frequency 
in a Z source is expected to be small compared to 
either the star's spin frequency or the Keplerian fre- 
quency at the gas coupling radius (see Lamb et al. 
1985; Shibazaki & Lamb 1987). This expectation is 
consistent with the ~ 250-350 Hz spin frequencies in- 
ferred from the differences between the frequencies of 
the QPOs in the kilohcrtz QPO pairs. 

In contrast to the comparatively low frequency of 
the HBO, the frequency of the lower-frequency QPO 
in the kilohcrtz QPO pairs (which is the beat fre- 
quency in the magnetospheric beat-frequency inter- 
pretation) is up to 0.7 times the frequency of the 
higher-frequency QPO (which in this interpretation 
is the Keplerian frequency where the gas in the disk 
couples to the magnetic field). It follows that in this 
interpretation the neutron star spin frequency is in 
some cases only ~ 0.3 times the Keplerian frequency 
at the coupling radius. Hence the frequencies of both 
of the two kilohcrtz QPOs should increase as Mi in- 
creases in the same way that i^ko increases, which is 
much more slowly than observed, as Figure 2 shows. 

To summarize, whereas the inferred dependence of 
the HBO on the mass flux through the inner disk 
is explained naturally by the magnetospheric beat- 
frequency model of the HBOs, the inferred depen- 
dence of the kilohcrtz QPOs on this mass flux is dif- 
ficult to explain using this model. 

Observed range of kilohcrtz QPO frequencies. The 
frequencies of the QPOs at what would be the Keple- 
rian orbital frequency in the magnetospheric beat-fre- 
quency interpretation all range up to ^ 1000-1200 Hz. 
This is understandable if these QPOs are generated 
within a few kilometers of the surface of the neutron 
star, as in the sonic-point model, whereas in the mag- 
netospheric beat-frequency model the frequencies of 
these QPOs would be expected a priori to range from 
^ 50 Hz up ~ 1500 Hz for the wide range of neutron 
star magnetic field strengths and accretion rates in- 
ferred in the atoll and Z sources. 

The expected range of orbital frequencies could be 
reduced if the mass accretion rate and the strength of 
the star's dipole magnetic field are tightly correlated 



in just the right way. However, the correlation that 
is indicated by spectral modeling (Psaltis et al. 1995; 
Psaltis & Lamb 1998a, 1998b, 1998c) gives magneto- 
spheric beat frequencies in the range ^ 20 200 Hz and 
orbital frequencies in the range ~ 300-500 Hz, much 
smaller than the frequencies of the kilohcrtz QPOs. 

Anticorrclation between kilohcrtz QPO amplitudes 
and inferred magnetic field strengths. In the magneto- 
spheric beat-frequency interpretation, the 
kilohcrtz QPOs should have higher amplitudes in sys- 
tems with stronger neutron star magnetic fields. How- 
ever, as discussed further in § 4.4, the rms amplitudes 
of the kilohcrtz QPOs in the atoll sources, which are 
thought to have the weakest neutron star magnetic 
fields, are considerably stronger than the rms ampli- 
tudes of the kilohcrtz QPOs in the Z sources, which 
are thought to have the strongest neutron star mag- 
netic fields. 

In the next section we describe a model that is 

consistent with both the physical picture of the atoll 
and Z sources developed over the past decade and 
with the properties of the kilohertz QPOs. 

3. THE SONIC-POINT MODEL 

We now analyze the physics of the sonic-point 
model of the kilohertz QPOs. In § 3.1 we describe 
the key elements of the sonic-point model and sum- 
marize the results of the general relativistic gas dy- 
namical and radiation transport calculations that are 
described in more detail later in this section. In § 3.2 
we investigate the motion of the gas in the inner disks 
of the Z and atoll sources and show that at a char- 
acteristic angular momentum loss radius i?ami that 
is typically several kilometers larger than the radius 
of the neutron star, radiation drag extracts angular 
momentum from the gas so quickly that the gas ac- 
celerates sharply inward. As a result, density fluctu- 
ations in the gas near the sonic point create brighter 
footprints that rotate around the stellar equator at 
the sonic point Keplerian frequency uks, producing a 
QPO at this frequency. We show further that a de- 
crease in i?ami with accretion rate is to be expected. 

In § 3.3 we show that the sonic-point Keplerian 
frequency is expected to be between ~300 Hz and 
^ 1200 Hz in all sources and should increase with in- 
creasing mass flux through the inner disk. The model 
therefore explains one of the most important features 
of the kilohertz QPOs. We also show that the sonic- 
point mechanism naturally produces a second QPO at 
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the beat frequency between i^ks and the stellar spin 
frequency I'spin- In § 3.4 we show that this mechanism 
can produce QPOs with v/Sv ~ 100, comparable to 
the largest Q values observed. In § 3.5 we consider 
radiation transport in the footprints and surround- 
ing hot central corona and show that the amplitude 
of the QPO at the sonic-point Keplerian frequency is 
expected to increase steeply with increasing photon 
energy. Also, the Keplerian-frequency oscillations at 
higher photon energies are expected to be delayed by 
~ 10 /xsec relative to the oscillations at lower ener- 
gies, to the extent that the effects of Comptonization 
dominate. In § 3.6 we consider the effect of gas sur- 
rounding the neutron star and magnetosphere on the 
amplitudes of the kilohcrtz QPOs and demonstrate 
that QPO amplitudes as large as ~ 20% are under- 
standable in the sonic-point model. 

The light travel time across a neutron star is a 
fraction < 10^^^ of the time required for accretion to 
change the star's mass and spin rate, so the exterior 
spacetime is stationary to extremely high accuracy. 
The spacetime outside a steadily and uniformly rotat- 
ing axisymmetric star with gravitational mass M and 
angular momentum J is unique to first order in the di- 
mensionless angular momentum parameter j = J /IvP 
(Hartle & Thorne 1968) and is the same as the Kerr 
spacetime to this order (the spacetime outside a rotat- 
ing star differs from the Kerr spacetime in second and 
higher orders; see, e.g.. Cook, Shapiro, & Teukolsky 
1994). All, or almost all, of the atoll and Z sources 
appear to have spin frequencies < 350 Hz (see § 2 and 
§ 4), in which case j is < 0.3 for these stars (see § 5.2). 
Hence the Kerr spacetime is a reasonably accurate ap- 
proximation to the exterior spacetime. Therefore, in 
this work we use the familiar Boyer-Lindquist coordi- 
nates. 

All analytical expressions given in this paper are 
accurate only to first order in j. These expressions 
would therefore be the same if written in terms of the 
circumferential radius. Calculation of physical quan- 
tities to higher order is straightforward, given a stellar 
model and interior and exterior metrics valid to the 
higher order, but no analytical expressions are known, 
so such calculations must be carried out numerically. 
For conciseness, in this section (only) we use units in 
which c — G = ks = 1, where c is the speed of light, G 
is the gravitational constant, and fc^ is Boltzmann's 
constant. 



3.1. Physical Picture and Summary of Calcu- 
lations 

Before presenting our calculations of the most im- 
portant elements of the sonic-point model, we first 
describe the model in more detail, summarizing the 
results of our calculations and indicating where these 
results are presented in the subsections that follow. 
We shall assume that the spin axes of the neutron 
stars in the kilohertz QPO sources are closely aligned 
with the rotation axes of their accretion disks. This 
is expected to be the case in these LMXBs, because 
mass transfer is expected to produce an accretion disk 
in which the gas circulates in the same sense as the 
orbital motion of the system. The torque on the neu- 
tron star created by such a disk will align the spin 
axis of the star with the axis of the accretion disk in 
a time short compared to the duration of the mass 
transfer phase in such systems (see Daumerie et al. 
1996). 

Accretion flow near the star. — As discussed in § 2.2, 
if the stellar magnetic field is of intermediate strength 
we expect a Keplerian disk flow to penetrate inside 
the magnetosphere. What happens to the flow there 
depends on the magnitude of the radiation drag force 
(Miller & Lamb 1993, 1996) and on whether the ra- 
dius R of the neutron star is larger or smaller than 
the radius i?,ns of the innermost stable circular orbit. 
There are four possibilities (see Table || for a sum- 
mary) . 

Suppose first that R < i?ms- This will be the 
case if the equation of state of neutron-star mat- 
ter is relatively soft and the star in question has an 
intermediate-to-high mass and is not spinning near 
its maximum rate. There are then two possibilities 
(§ 3.2): 

(la) If the drag force exerted by radiation coming 
from near the neutron star's surface is strong enough 
to remove > 1% of the angular momentum of the gas 
in the Keplerian disk flow just outside Rms, gas in 
the inner disk drifts slowly inward in nearly circular 
orbits until it reaches the critical angular momentum 
loss radius Rami, where it transfers so much angular 
momentum to the radiation so rapidly that centrifu- 
gal support fails and the radial velocity of the gas 
increases sharply. Inside i?ami the gas spirals inward 
in a disk flow in which the radial velocity is supersonic 
until the gas is channeled by the star's magnetic field 
or collides with the star's surface, 
(lb) If instead the radiation drag is so weak that 
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it removes <C 1% of the angular momentum of the 
gas in the Keplerian disk flow by the time the gas 
reaches -Rms, gas in the inner disk drifts slowly in- 
ward in nearly circular orbits until it approaches Ry^s , 
where general relativistic corrections to Newtonian 
gravity cause the radial velocity of the flow to increase 
steeply (see, e.g., Muchotrzcb 1983; Muchotrzcb- 
Czerny 1986). Inside i?msi the gas spirals inward in 
a disk flow in which the radial velocity is supersonic 
until the gas is channeled by the star's magnetic field 
or collides with the star's surface. 
Whether case la or lb applies depends primarily on 
the strength of the star's magnetic field and on the 
mass accretion rate. 

Suppose now that R > i?ms- This will be the 
case if the equation of state of neutron-star mat- 
ter is relatively stiff and the star in question has an 
intermediate-to-low mass or is spinning very rapidly. 
There are again two possibilities (§ 3.2): 
(2a) If the radiation drag force is at least moderately 
strong, gas in the inner disk drifts slowly inward until 
it reaches the critical radius i?ami, where it accelerates 
rapidly inward. Inside i?ami the gas spirals inward in 
a disk flow in which the radial velocity is supersonic 
until the gas is channeled by the star's magnetic field 
or collides with the star's surface. 
(2b) If instead the radiation drag is weak, gas in the 
inner disk drifts slowly inward in a disk flow in which 
the radial velocity is supersonic until the gas is chan- 
neled by the star's magnetic field or interacts directly 
with the star's surface. 

Again, whether case 2a or 2b applies depends primar- 
ily on the strength of the star's magnetic field and on 
the mass accretion rate. 

Regardless of whether the slow inward drift of the 
gas in the disk is terminated by loss of angular mo- 
mentum to the radiation field or because the gas 
crosses the radius of the innermost stable circular or- 
bit, the inward radial velocity of the gas increases 
abruptly at i?ami or i?ms and becomes supersonic 
within a very small radial distance. As mentioned 
in § 1 and earlier in this section and described in 
detail below, the radius where either radiation drag 
or general relativistic corrections to Newtonian grav- 
ity cause the radial velocity to increase sharply and 
become supersonic plays a key role in the model of 
the kilohertz QPOs proposed here. The radius of the 
sonic point is a useful indicator of the location of this 
transition. 



Radiation drag can remove at most only a frac- 
tion of the specific angular momentum of the gas and 
therefore can create a radially supersonic inflow in the 
inner disk only within a few stellar radii (§ 3.2). The 
reason is that the specific angular momentum of gas 
in circular Keplerian orbits at radii much larger than 
the radius of the neutron star is much greater than in 
orbits just outside the stellar surface, so removal of a 
small fraction of the angular momentum of the gas at 
large radii cannot cause a large fraction of it to fall to 
the stellar surface. In contrast, the specific angular 
momentum of gas orbiting near the star is not much 
greater than the angular momentum of gas orbiting 
at the stellar surface, so radiation drag can cause gas 
in orbit near the star to plunge to its surface. 

Radiation drag can cause gas to plunge inward 
from a radius larger than one would estimate us- 
ing the Newtonian approximation, because (1) spe- 
cial and general relativistic effects on the gas dynam- 
ics and radiation transport significantly increase the 
fraction of the angular momentum of the accreting 
gas that can be removed by the radiation and (2) in 
general relativity the specific angular momentum of 
gas in circular Keplerian orbits with radii r < 3 i?ms 
varies only slowly with r. For this reason, the sonic- 
point model predicts that if radiation drag produces 
a transition to rapid radial inflow, the transition will 
occur within ~ 36 km (for a 1.4 Mq neutron star) and 
hence that the Keplerian frequency at the transition 
radius win be > 300 Hz (§ 3.3). 

Inside the sonic point, the vertical optical depth 
of the disk flow falls steeply with decreasing radius, 
usually to a value small compared to unity. The op- 
tical depth of the disk flow measured radially from 
the stellar surface to the sonic radius typically also 
becomes less than unity within a very small radial 
distance, unless the accretion rate is very high or the 
geometrical thickness of the disk in the vertical direc- 
tion is very small. When the change from "optically 
thick" to "optically thin" disk flow is caused by the 
radiation drag force, the transition is somewhat anal- 
ogous to the ionization front at the boundary of an 
Hii region. In the disk flow transition, the photon 
mean free path becomes longer because the radiation 
is removing angular momentum and the flow is ac- 
celerating inward, causing the density to fall sharply, 
whereas in an ionization front the radiation is remov- 
ing bound electrons from atoms and molecules, caus- 
ing the opacity to fall sharply. 

Generation of the QPO at the sonic-point Keple- 
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rian frequency. — Suppose first that the magnetic field 
of the neutron star is dynamically negligible. An el- 
ement of gas in the disk flow outside the sonic ra- 
dius diffuses inward subsonically. When it reaches 
the sonic radius, it falls supersonically to the stellar 
surface along a spiral trajectory like that shown in 
Figure This trajectory was computed in full gen- 
eral relativity for a nonrotating, isotropically emit- 
ting star, using the numerical algorithm described in 
Miller & Lamb (1996). (The shape of the trajectory 
depends on the luminosity, spin, and other proper- 
ties of the source; the example shown in Fig. ^a is 
only illustrative.) The surface density of the disk flow 
is much smaller inside the sonic radius than outside 
because of the sharp increase in the radial velocity 
at the sonic point. Gas falls inward from the sonic 
radius and impacts the star all around its equator, 
producing a bright equatorial ring of X-ray emission. 

The flow in the inner part of the accretion disk is 
expected to have density fluctuations ( "clumps" ) pro- 
duced by a variety of mechanisms, such as thermal in- 
stability, Kelvin-Helmholtz instability, and magneto- 
turbulence (see Lamb et al. 1985; Shibazaki & Lamb 
1987). There may be as many as several hundred 
such clumps at a given radius in the disk outside 
the sonic point. The velocity-independent radially 
outward component of the radiation force and the 
velocity-dependent radiation drag force both tend to 
cause gas in the shadow of a clump to overtake the 
clump in azimuth and radius, and hence the radiation 
force tends to increase clumping. 

As discussed in § 3.4, a clump that forms outside 
the sonic radius is dissipated by gas pressure forces, 
turbulence, and the shear in the velocity field before 
the gas from it can reach the stellar surface, so such 
clumps do not produce significant inhomogeneities in 
the infiow from the sonic radius. In contrast, a clump 
that forms near the sonic radius can persist for up to 
^ 100 times the infall time to the stellar surface, so 
such a clump generates a stream of denser infalling 
gas from the sonic point to the stellar surface that 
lasts as long as the clump survives. 

The shape of the pattern formed by the denser gas 
falling inward from a clump depends on how the an- 
gular velocity of the infalling gas varies with radius. 
For example, if the angular velocity of the gas were in- 
dependent of radius, the pattern of higher gas density 
would be a straight, radial line from the clump to the 
stellar surface. In general, the pattern formed by the 
denser gas from a clump is a fairly open curve. The 



shape of the pattern, like the trajectories that pro- 
duce it, depends on the luminosity, spin, and other 
properties of the source. Figure ^ shows the pattern 
formed by inflow of gas along spiral trajectories with 
the shape shown in Figure In this example the 
pattern of denser gas is also a spiral. It is more open 
than the spiral trajectories of the individual elements 
of denser gas that produce it, because the source of 
the denser gas (the clump) orbits the star only slightly 
slower than the infalling gas. 

We expect the time-averaged radiation field around 
the star to be nearly axisymmetric. Moreover, the gas 
from a given clump typically orbits the star ^ 5-10 
times before colliding with the stellar surface, so the 
effects of any azimuthal variations in the radiation 
drag force are averaged out. Therefore, the spiral 
pattern of denser gas produced by the infiow from 
a given clump rotates nearly uniformly around the 
star with a rotation frequency equal to the orbital 
frequency of the clump that is producing it, which 
is the Keplerian frequency at a radius near the sonic 
point. 

The time evolution that generates an X-ray oscil- 
lation at the sonic-point Keplerian frequency is illus- 
trated by the four panels of Figure ^. In this fig- 
ure a single clump is shown advancing in its orbit 
by 90° from one panel to the next. The pattern of 
higher density gas rotates uniformly around the star 
with a frequency equal to the orbital frequency of the 
clump, so it also advances by 90° from one panel 
to the next. Where the denser gas from a clump 
collides with the stellar surface, it produces an arc- 
shaped area of brighter X-ray emission. This arc- 
shaped brighter "footprint" moves around the star's 
equator with a frequency equal to the rotation fre- 
quency of the pattern, which is the orbital frequency 
of the clump and is therefore approximately equal to 
the Keplerian frequency vks at the sonic point. In 
reality, many clumps are crossing the sonic radius at 
any given time, and hence there are many bright foot- 
prints moving around the star's equator. The radia- 
tion from these footprints carries the kinetic energy 
of infall that is released when the gas that has fallen 
inward from the clumps at the sonic radius collides 
with the stellar surface. 

As seen by a distant observer whose line of sight is 
inclined with respect to the orbital axis of the disk, 
the aspect presented by a given bright footprint varies 
periodically and the footprint is eclipsed with a fre- 
quency equal to the rotation frequency of the pat- 
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tern, which is approximately the Keplerian orbital fre- 
quency at the sonic point. Footprints come and go as 
clumps form near the sonic point and then dissipate. 
Also, the orbital frequencies of the clumps that are 
producing footprints at any given time differ slightly. 
As a result, a distant observer sees a strong, quasi- 
periodic oscillation of the X-ray flux and spectrum at 
the pattern rotation frequency, which is close to i^ks- 
This is the the sonic-point Keplerian-frequency QPO. 

The spin of the neutron star is not involved in 
generating the sonic-point Keplerian frequency QPO. 
The precise frequency of the sonic-point Keplerian 
QPO does depend weakly on the spin frequency of 
the neutron star, because the frame-dragging caused 
by the star's spin has a small effect on the frequencies 
of Keplerian orbits near the neutron star. The sonic- 
point Keplerian frequency QPO mechanism described 
here, which generates a QPO with a frequency equal 
to the orbital frequency at the sonic point, differs fun- 
damentally from the suggestion of Paczynski (1987), 
who speculated that when the sonic point is located 
at the radius of the marginally stable orbit, a QPO 
might be produced by unsteady flow through it with 
a low (~ 20-50 Hz ) frequency that would have noth- 
ing directly to do with the orbital frequency of the 
marginally stable orbit. 

Generation of the QPO at the sonic-point beat 
frequency. — Suppose now that the magnetic field of 
the neutron star is weak enough that a prograde Ke- 
plerian disk flow penetrates near the surface of the 
star but strong enough that close to the star some of 
the gas is channeled by the field. The channeled gas 
produces slightly brighter spots where it collides with 
the stellar surface. If these slightly brighter spots are 
offset from the star's spin axis, the enhanced radia- 
tion from them generates a beamed pattern of radi- 
ation that rotates with the star. This rotating ra- 
diation pattern creates a periodic oscillation in the 
radiation force that is acting on the clumps of gas 
crossing the sonic radius. Hence, if conditions near 
the star are such that the nearly Keplerian motion of 
the gas in the clumps orbiting near the sonic radius is 
terminated by loss of angular momentum to the radi- 
ation, the inward flux of mass from the clumps near 
the sonic radius will increase once or twice each beat 
period, depending on the symmetry of the radiation 
pattern that rotates with the star. 

As a specific example, suppose that the radiation 
drag force experienced by the gas in a given clump or- 
biting near the sonic radius peaks each time the beam 



of radiation that rotates with the star sweeps across 
the clump. This will occur with a frequency equal 
to one or two times the difference between the orbital 
frequency of the clump, which is the Keplerian orbital 
frequency near the sonic point, and the spin frequency 
i^spin of the neutron star. We call ^ks — t'spin the sonic- 
point beat frequency and denote it i/bs- If the transi- 
tion to supersonic inflow is caused by radiation drag, 
the oscillation of the radiation drag force will cause 
the supersonic inflow of gas from the clumps orbit- 
ing near the sonic point to oscillate quasi-periodically 
with a frequency equal to one or two times the sonic- 
point beat frequency. 

The time evolution that produces this oscillation is 
shown schematically in the six snapshots of Figure ^ 
For clarity, only a single clump is shown orbiting at 
the sonic radius and the six snapshots show the se- 
quence of events as seen in a frame corotating with 
the star. In panel (a) the clump is in the beam of 
radiation that rotates with the star and hence the in- 
ward flux of denser gas from the clump is greater than 
average. In panel (b) the clump has moved out of the 
beam and hence the inward flux of gas from the clump 
is smaller than average; the pulse of gas that was 
dragged off the clump by the stronger radiation flux 
in frame (a) is now falling inward from the clump. In 
panel (c) the clump has moved ahead of the beam by 
180°; the inward flux of gas from the clump remains 
smaller than average and the pulse of denser gas that 
began to fall inward in frame (a) is further away from 
the clump. In panel (d) the clump is now ahead of 
the beam by 270°; the inward flux of gas from the 
clump remains smaller than average and the pulse of 
denser gas that began to fall inward in frame (a) is ap- 
proaching the stellar surface. In panel (e) the clump 
has now moved back into the beam of radiation and a 
new pulse of denser gas is leaving the clump; the pulse 
of denser gas that began to fall inward in frame (a) is 
now very close to the stellar surface. In panel (f) the 
clump has again moved ahead of the beam by 90° and 
the inward flux of gas from the clump is again smaller 
than average; the pulse of denser gas that began to 
fall inward in frame (a) is now colliding with the stel- 
lar surface, producing a brighter beam of radiation 
from the footprint of the stream. 

In reality, the radiation drag force, and hence the 
inward mass flux from a clump, may be greatest not 
when the clump is fully illuminated by the beam that 
rotates with the star, as assumed in this illustration, 
but at some other relative phase, because the radia- 
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tion drag force depends in a complicated way on the 
various components of the radiation stress-energy ten- 
sor (see Miller & Lamb 1996). In any case, the inward 
mass flux from the orbiting clumps oscillates quasi- 
periodically with a frequency approximately equal to 
the sonic-point beat frequency ubs, producing a quasi- 
periodic oscillation in the luminosity and spectrum of 
the X-ray emission from the stellar surface. This is 
the sonic-point beat-frequency QPO. 

Properties of the sonic-point QPOs. — In § 3.3 we 
show that the sonic-point mechanism naturally gen- 
erates Keplerian-frequency QPOs with frequencies 
^ 300-1200 Hz for a wide range of stellar magnetic 
fields and accretion rates and that it is natural for 
the frequencies of the sonic-point Keplerian and beat- 
frequency QPOs to increase steeply with increasing 
mass flux through the inner disk. The angular dis- 
tribution of the emission from each emitting point 
within a footprint at the stellar surface is broad, so 
the power at is likely to be much greater than 
the power at overtones of uks- The modulation of 
the inflow from clumps orbiting at the sonic radius 
is unlikely to be perfectly sinusoidal, so overtones of 
the sonic-point beat frequency i/bs as well as ubs itself 
may be detectable. 

In § 3.4 we argiie that the sonic-point mecha- 
nism can produce QPOs with i^/6i^ ratios as large 
as 100, i.e., similar to those observed in the atoll 
and Z sources. The oscillation at the sonic-point beat 
frequency is generated by the beating of the stellar 
spin frequency, which is periodic, against the quasi- 
periodic sonic-point Keplerian frequency oscillation. 
It is therefore natural to expect the ly/Siy values of 
the oscillations in a pair to be roughly similar, in 
the absence of other disturbing effects. In § 3.5 we 
show that the steep increase in QPO amplitude with 
photon energy observed in many kilohertz QPOs is 
understandable in the sonic-point model. 

The visibilities of the various QPOs generated by 
the sonic-point mechanism depend on the mass flux 
in the Keplerian disk flow that penetrates to the sonic 
point, the number and distribution of clumps in the 
disk at the sonic point, the brightness of the footprints 
produced by the streams from the chimps relative to 
the brightness of the rest of the stellar surface, the 
geometry and optical depth of the scattering material 
around the star, and the inclination of the system. 

If radiation forces are required to produce large 
clumps with substantial density contrasts, the ampli- 
tudes of the QPOs at the sonic-point Keplerian and 



beat frequencies are likely to become much smaller 
and may become undetectable once the angular mo- 
mentum loss radius i?ami has retreated inward to the 
radius R^s of the innermost stable circular orbit, be- 
cause once this has happened, the gas in the disk near 
the sonic point is increasingly shielded from radiation 
coming from the stellar surface by gas in the disk 
closer to the star. However, if large clumps with sub- 
stantial density contrasts are formed by other mecha- 
nisms, the sonic-point Keplerian frequency QPO may 
still be detectable if the sonic point is at Rms, even if 
the sonic-point beat frequency QPO — which can only 
be generated if radiation forces have a significant dy- 
namical effect at the sonic radius — has disappeared. 

In § 3.6 we show that scattering by the electrons in 
the central corona that surrounds the neutron stars 
in the Z and atoll sources strongly attenuates the al- 
ready intrinsically weak beaming oscillation at the 
stellar spin frequency as well as the weak beaming 
oscillations at other frequencies, making all but the 
strongest beaming oscillation — the one at the sonic- 
point Keplerian frequency — difficult to detect with 
current instruments. In contrast, scattering only 
weakly attenuates luminosity oscillations, so the lu- 
minosity oscillations at the sonic-point beat frequency 
and its overtones are likely to be detectable. Tak- 
ing into account both the generation of these vari- 
ous oscillations and their attenuation in the circum- 
stellar environment, Keplerian-frequency and beat- 
frequency QPOs with rms amplitudes as large as 15% 
appear possible. 

As discussed in § 2, the Z sources appear to have 
significantly stronger magnetic fields than the atoll 
sources. The magnetic fields of the Z sources are 
therefore likely to channel a larger fraction of the ac- 
creting gas out of the disk plane before it penetrates 
close to the stellar surface, so we expect the ampli- 
tudes of the kilohertz QPOs to be smaller in the Z 
sources than in the atoll sources. The magnetic field 
of the neutron star in some atoll sources may be so 
weak that the QPO at the sonic-point beat frequency 
is undetectable, even though a QPO at the sonic-point 
Keplerian frequency is visible. The sonic-point beat 
frequency may also be undetectable at some times if 
the scattering optical depth between the stellar sur- 
face and the sonic point is too large. In other sources 
or at other times, the luminosity oscillation at the 
sonic-point beat frequency may be detectable even 
though the beaming oscillation at the sonic-point Ke- 
plerian frequency has been suppressed by scattering 



16 



in the corona surrounding the star. 

3.2. Transition to Rapid Radial Inflow 

In both the atoU and the Z sources, we expect ra- 
diation drag to cause gas in the Keplerian disk flow to 
make an abrupt transition from slow inward drift to 
rapid radial inflow several kilometers above the sur- 
face of the neutron star. This transition occurs at the 
radius where the drag exerted by radiation from the 
star removes enough angular momentum from the gas 
quickly enough that it falls to the surface of the star 
unimpeded by a significant centrifugal barrier. We 
first describe the consequences of the transfer of an- 
gular momentum from the gas in the inner disk to the 
radiation and then present a fully general relativistic 
calculation of the gas dynamics and radiation trans- 
port in the inner disk that shows the nature of this 
transition in the atoll sources, which have luminosi- 
ties much less than the Eddington critical luminosity. 
We next give approximate analytical expressions for 
the location and width of this transition, then discuss 
why the nature and location of the transition may be 
similar in the Z sources even though their accretion 
rates and luminosities are much higher than the ac- 
cretion rates and luminosities of the atoll sources, and 
finally summarize the implications for the sonic-point 
QPOs. 

Transfer of angular momentum to the radiation 
field. — The azimuthal drag force exerted by radiation 
from the star can create a radially supersonic inflow in 
the inner disk. However, this is possible only within 
a few stellar radii if the luminosity of the star is pro- 
duced by accretion and most of the accreting gas flows 
through the inner disk. The reason is that radiation 
drag can remove at most a fraction 77i.ad(inax) r:!20% 
of the specific angular momentum of the accreting gas 
(see Fortner et al. 1989; for more detailed discussions, 
see Miller & Lamb 1993, 1996). The specific angular 
momentum of gas in circular Keplerian orbits at radii 
large compared to the radius R of the neutron star is 
much larger than at R, so radiation drag cannot cause 
more than ^--^20-30% of the gas at such large radii to 
fall to the stellar surface. In contrast, the specific an- 
gular momentum of gas in circular Keplerian orbits 
very near the star is not much larger than the an- 
gular momentum of orbits at the stellar surface, so 
radiation drag can cause gas in Keplerian orbits near 
the star to plunge to the surface. 

Radiation drag can create a radially supersonic in- 
flow at a larger radius than one would estimate us- 



ing the Newtonian approximation, because special 
and general relativistic effects significantly increase 
the fraction of the angular momentum of the accret- 
ing gas that is transferred to the radiation (Miller 
& Lamb 1993). Also, in general relativity the spe- 
cific angular momentum of gas in circular Keplerian 
orbits with radii r < 3 i?ins varies only slowly with ra- 
dius. The latter point is illustrated by Figure^, which 
shows the specific angular momentum of gas in circu- 
lar Keplerian orbit around a nonrotating and a slowly 
rotating 1.4 Mq neutron star, as a function of the ra- 
dius of the orbit. The specific angular momentum 
varies fairly steeply at r ^ i?nis , but only slowly with 
radius between r = Rms and r ~ 3i?ms- For exam- 
ple, the specific angular momentum of circular orbits 
changes by less than ^ 10% from r — 12M (25 km) 
to r = 6M (12 km). These curves are of course physi- 
cally meaningful only at radii greater than the radius 
R of the neutron star. The relatively flat shape of 
the angular momentum curve for circular orbits near 
a neutron star means that if gas orbiting there loses 
even 20% of its angular momentum, it will plunge 
to the stellar surface. 

We can estimate the outer boundary of the region 
where radiation coming from near the stellar surface 
can remove sufficient angular momentum from the 
gas orbiting in the disk so that it falls supcrsonically 
to the stellar surface by comparing T^j-ad (max) , the 
largest fraction of angular momentum that can be re- 
moved by radiation coming from the surface of the 
neutron star, with r/flow, the fraction of the specific 
angular momentum of gas in a circular Keplerian or- 
bit at radius r that must be removed in order for the 
gas to fall from r to R^s (if the radius of the star is 
less than i?ms, gas that has reached i?rns can fall to 
the stellar surface without losing any more angular 
momentum). These two quantities are compared in 
Figure |^. Here yyrad (max) was computed by assuming 
that the gravitational binding energy at the stellar 
surface is converted into radiation there, taking into 
account the relevant angular, special relativistic, and 
general relativistic factors (see Miller & Lamb 1996), 
and assuming that all of the radiation emitted in 
the direction of the disk interacts with the disk fiow 
(/ = 1 in the notation of Miller & Lamb [1993]). 
Figure ^ shows that ?7rad (max) exceeds Tyflow only in- 
side ~ 17 M (~ 36 km). This characteristic radius de- 
pends on the radius of the neutron star, because the 
luminosity is proportional to the gravitational bind- 
ing energy at the stellar surface, but this dependence 
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is relatively weak (the characteristic radius decreases 
from 21M to 13M as the stehar radius increases from 
4M to 7M). Hence, the sonic-point model predicts 
that if radiation drag produces a transition to rapid 
radial inflow, the transition will occur within ^ 15 M 
and hence that the Keplerian frequency at the transi- 
tion radius will be > 300 i/z (see § 3.3). Determining 
the actual radius of the transition requires gas dy- 
namical and radiation transport calculations. 

Model calculations. — The nature of the transition 
to supersonic radial inflow is illustrated by the follow- 
ing fully general relativistic calculations of the gas dy- 
namics and radiation transport in the innermost part 
of the accretion disk flow. In these calculations the 
azimuthal velocity of the gas in the disk is assumed to 
be nearly Keplerian far from the star. Internal shear 
stress is assumed to create a constant inward radial 
velocity in the disk, as measured in the local static 
frame, of 10^^ (here and below we take v"^ to be pos- 
itive for radially inward flow). The half-height h{r) 
of the disk flow at radius r is assumed to be er at 
all radii, where e is a constant and r is the radius, 
and the kinetic energy of the gas that collides with 
the surface of the star is assumed to be converted to 
radiation that emerges from a band around the star's 
equator with a half-height equal to eR. For simplic- 
ity, and to show the effects of radiation forces more 
clearly, any effect of the stellar magnetic field on the 
dynamics of the disk flow near the sonic transition is 
neglected in this model calculation. 

The surface density Sco = 'i.hpco of the disk flow as 
measured by a local observer comoving with the gas 
is related to the surface density Sstat of the disk flow 
as measured by a local static observer by a Lorcntz 
boost, so 

Sco = 7^"^Sstat , (3) 

where 7 = [1 — (w'')^ — [v'^Y]"^/'^ and we have ne- 
glected w^. For stationary disk accretion, Sgtat 
radius r is given by (see Novikov & Thornc 1973, 
eq. 5.6.3) 

Sstat = (Mco/2W)(l - 2M/r)-i/2 ^ (4) 

where Moo is the mass accretion rate through the disk, 
as measured at infinity. 

Close to the star, the photon mean-free-path is lim- 
ited primarily by Thomson scattering and is therefore 
l/ueCTT, where rig is the electron number density in 
the disk flow as measured in the frame comoving with 
the flow and ctt = 6.65 x 10^^''' cm^ is the Thomson 



scattering cross section. Hence, for a steady flow, the 
radial optical depth from the stellar surface at radius 
R through the disk to radius r is 

Tradiai(r) = ar ne{rP)drP = 
Jr 

aTM^\ r [l~v-{r')]-^ il-2M/r'r\ ^> 
Anmp I J "^r{f')l{f') r'h{r')v^ {r') ' 

(5) 

where Mi^oo is the mass accretion rate through the 
inner disk, as measured at infinity. The radial coor- 
dinate r'P in the first expression on the right is the 
proper radial distance in the frame comoving with 
the gas. In the second expression on the right we 
have estimated in the disk flow between r and 
the stellar surface by assuming for simplicity that the 
gas in the disk is fully ionized hydrogen and using 
the continuity equation. Then = Sco/2/irnp = 
7"-^(Afi_oo/47rmpi;'"r/i)(l — 2A//r)^-^/-^, where nip is 
the proton mass and h is the half-thickness of the 
disk flow. The differential proper radial distance dr^ 
in the frame comoving with the gas is related to the 
differential radial distance in the local static frame 
by a Lorentz boost, so dr^ = (1 — w'')^^7r~"^drstat, 
where is the inward radial velocity and 7^ = 
[1 — {v^Y]~^/'^. The differential radial distance in the 
local static frame is in turn related to the differential 
radial distance dr' in the global (Boyer-Lindquist) co- 
ordinate system by dr^tat = (1 — 2M /r')~'^/'^dr' , so 
drP = (1 - v^)-^ 7^-1 (1 - 2M/r')"^/^ dr' . 

Once the drag force exerted by the radiation from 
the stellar surface begins to remove angular momen- 
tum from the gas in the Keplerian disk, centrifugal 
support is lost and the gas falls inward, accelerat- 
ing rapidly. Radiation that comes from near the star 
and is scattered by the gas in the disk is generally 
scattered out of the disk plane and hence does not 
interact further with the gas in the disk. Moreover, 
second and successive scatterings do not contribute 
proportionately to the azimuthal radiation drag force 
on the gas because the radiation fleld is aberrated 
by the first scattering and afterward carries angu- 
lar momentum (see Miller & Lamb 1993, 1996). We 
therefore treat the interaction of the radiation with 
the gas in the disk by assuming that the intensity of 
the radiation coming from the star is attenuated as 
it passes through the gas in the disk, diminishing as 
exp(-Ti.adiai), where rradiai(r) is given by equation (||), 
and that scattered radiation does not contribute to 
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removal of angular momentum from the gas. In cal- 
culating the radiation drag force, we assume for sim- 
plicity that the differential scattering cross section 
is isotropic in the frame comoving with the accret- 
ing gas (this gives results very close to those using a 
Thomson differential cross section; see Lamb & Miller 
1995). The radiation field and the motion of the gas 
are computed in full general relativity. 

The optical depth of the disk flow near the star is 
generally much less than one would estimate by cal- 
culating it without taking into account the radiation 
forces. The reason is that the loss of centrifugal sup- 
port caused by transfer of angular momentum to the 
radiation causes the inward radial velocity of the gas 
in the disk at radii less than the radius -Rami of the 
radiation-induced transition to supersonic inflow to 
be orders of magnitude higher, and the density of the 
gas to be orders of magnitude lower, than if radia- 
tion forces were neglected (see Miller & Lamb 1996). 
Hence the mean free path in the disk flow is much 
larger, and the optical depth from the stellar surface 
to a given radius in the flow is much smaller, than 
they would be in the absence of radiation drag. In 
this way the radiation increases the transparency of 
the disk flow, so its effects are felt much farther out in 
the flow than one would estimate from the properties 
of the undisturbed flow. 

Figure || shows results obtained by solving self- 
consistently for the gas dynamics and the radiation 
field in this simple model, for a disk flow of semi- 
thickness e = h/r ~ 0.1, a nonrotating neutron 
star with gravitational mass M — 1.4 M© and radius 
R = 5M, and accretion rates through the inner disk of 
0.02, 0.03, 0.04, and 0.05 Mb, where Me is the mass 
accretion rate that produces the Eddington critical lu- 
minosity at infinity (see Lamb & Miller 1995). These 
accretion rates are typical of the lower- luminosity 
atoll sources. The calculations were carried out using 
the numerical algorithm described in Miller & Lamb 
(1996). 

The results shown in Figure ^ demonstrate that ra- 
diation forces and general relativistic effects create a 
sharp transition to supersonic inflow several kilome- 
ters above the stellar surface and that the transition 
radius decreases with increasing mass flux when the 
transition is caused by radiation drag. These results 
are not intended to represent the inflow in any par- 
ticular source. For Mi = 0.02 Mij and 0.03Mb, the 
transition to rapid radial inflow occurs at 7.3 M and 
6.3 M, respectively, and is caused by transfer of an- 



gular momentum from the gas to the radiation field. 
For ili = 0.04 Mb and 0.05 Mb, most of the radia- 
tion from the star does not penetrate as far out as 
the radius i?ms = 6M of the innermost stable circular 
orbit, so the transition to rapid radial inflow occurs 
at i?,„s and is caused by the absence at r < _Rms of 
circular orbits with the specific angular momentum 
of the orbit at r = i?ms- 

The radial velocity profiles plotted in Figure ||a 
show that at the sonic point the radial velocity typ- 
ically increases by about two orders of magnitude 
in a very small radial distance (Jr ~ 0.001-0.01 r). 
The angular velocity profiles plotted in Figure 
show that the effect of radiation forces on the an- 
gular velocity is smaller, but is still very significant. 
For Mi = 0.04 Mb and 0.05 Mb, the angular veloc- 
ities outside 6M are nearly Keplerian, because the 
gas there is shielded from the radiation by the gas 
further in. However, for M, = 0.02 Me and 0.03Mb, 
the angular velocity drops below the Keplerian value 
well outside 6M . It is this departure from Keplerian 
orbital motion, and the associated loss of centrifugal 
support, that causes the flow to accelerate radially 
inward. All the angular velocity profiles turn down- 
ward near the stellar surface because of the strong 
azimuthal radiation drag force there. The decrease in 
the kinetic energy and angular momentum of the in- 
flowing gas is exactly compensated by the increase in 
the energy and angular momentum carried outward 
by the escaping radiation (see Miller & Lamb 1993, 
1996). 

Figure ||c shows that the vertical column density of 
all four disk flows drops abruptly from ^ 10^ gcm~^ 
outside the sonic radius to < 10 gcm^^ inside. The 
radial optical depth measured outward from the stel- 
lar surface is shown in Figure ||d; it decreases by about 
two orders of magnitude at the sonic radius, but not 
as sharply as the vertical optical depth. Even when 
the transition to rapid radial inflow occurs at i?ins, 
the flow between R^s and the stellar surface is gener- 
ally strongly affected by radiation forces, as shown by 
the differences between the velocity and optical depth 
profiles for M, = 0.04 Mb and 0.05 Mb- We caution 
that these results are only illustrative. For example. 
X-rays may be emitted from a larger fraction of the 
stellar surface than we have assumed, in which case 
-Rami will be larger for a given accretion rate. 

Radius of the transition. — In the model of gas dy- 
namics and radiation transport just described, under 
conditions such that loss of angular momentum to 



19 



radiation drag is important outside i?ms, loss of cen- 
trifugal support and rapid radial inflow begins about 
five photon mean-free-paths from the stellar surface, 
i.e., the transition occurs where T,-adiai{r) ~ 5. In 
this model the inward radial velocity of the gas in 
the disk inside i?ami typically rises very sharply to 
~ 0.1 c and then changes more slowly until the gas col- 
lides with the stellar surface (see also Miller & Lamb 
1993, 1996). Hence we can estimate i?ami by setting 
Tradiai = 5 in equation (||), scaling and h/r in units 
of 10~^ c and 0.1, respectively, and solving for the 
radius. The result is 

This estimate is in rough agreement with the transi- 
tion radius found in our numerical calculations. 

Width of the transition. — Under conditions such 
that the transition to rapid radial inflow is caused 
by radiation drag, the width 5 of the radial velocity 
transition is determined by A^, the photon mean- free- 
path in the Keplerian disk flow at i?ami, because this 
is the characteristic distance over which the gas first 
becomes exposed to radiation from the star, and by 
the radial distance Jami over which the gas in the disk 
loses its angular momentum once it is exposed to ra- 
diation from the star. From expressions similar to 
those used above to compute the radial optical depth 
in the Keplerian disk flow, we find 

^ ^ 10- i " f ^ V 

r \ 0-1 ; \^3x 10-5cy ' 

(7) 

where we have scaled the quantities that enter this 
expression in terms of their approximate values in the 
transition region. This estimate is consistent with 
the results of our numerical calculations. From the 
expression for the radiation drag rate given by Miller 
& Lamb (1993, 1996), we find 

^aml ^ g ^3/2 / ]£_ \ 

r ^ ' \Me) 
where we have scaled Cg in terms of the typical value 



for a gas-pressure-dominated inner disk. These es- 
timates show that the radial velocity is expected to 
increase over a very small radial distance 6 ~ 10~^ r- 
10""* r. This expectation is in accord with the re- 
sults of our numerical calculations presented above. 
As noted in § 3.1, this transition is somewhat analo- 
gous to the transition at the boundary of a Stromgren 
sphere, except that here the radiation removes an- 
gular momentum from the gas rather than stripping 
electrons from atoms or molecules and causes a sharp 
increase in the inward radial velocity rather than in 
the degree of ionization. 

If the neutron star radius is smaller than R^^s and 
conditions are such that the transition to rapid ra- 
dial inflow is caused by general relativistic correc- 
tions to Newtonian gravity, the radial velocity in- 
creases rapidly near i?,„s because gas inside i?nis can 
fall to the stellar surface without losing any angular 
momentum. The precise location of the sonic transi- 
tion is slightly affected by pressure forces and, if the 
shear stress is large, by outward angular momentum 
transport within the disk flow (see Muchotrzeb 1983; 
Muchotrzeb-Czerny 1986). In the absence of shear 
stresses and pressure forces, the radial velocity of 
an element of gas at radius r = -Rms(l — C)^ where 
^ ^ 1, that has fallen inward from a circular orbit at 
i?ms is approximately (C/2)^/^ (Miller & Lamb 1996). 
Hence, for Cs ^ 10^'^, the inflow becomes supersonic 
over a radial distance 5 ~ [csY^^r ~ 10~^ r. 

Summary and implications for kilohertz QPOs. — 
For sources in which the radius R of the neutron star 
is larger than the radius i?nis of the innermost sta- 
ble circular orbit and accretion is under conditions 
such that angular momentum loss to the radiation 
field terminates the Keplerian flow at a radius larger 
than R, our results show that the Keplerian disk flow 
near the star is terminated by angular momentum loss 
to radiation and that the radius of the transition to 
rapid radial inflow decreases steeply with increasing 
mass flux through the inner disk. Hence the orbital 
frequency at the sonic radius, and therefore the fre- 
quency of the associated QPO, increases steeply with 
increasing mass flux. If instead conditions are such 
that radiation forces do not terminate the Keplerian 
flow, the Keplerian flow is probably terminated by in- 
teraction with the stellar surface, in which case gen- 
eration of a coherent QPO with a frequency equal to 
the orbital frequency at the stellar surface appears 
very unlikely. 

For sources in which R is smaller than i?,„s and 
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accretion occurs under conditions such that radia- 
tion forces terminate the flow outside i?nis , our results 
show that the Keplerian flow is terminated by angu- 
lar momentum loss to radiation and that the radius of 
the transition to rapid radial inflow decreases steeply 
with increasing mass flux. Hence the frequency of 
the QPO associated with the orbital frequency at the 
sonic radius again increases steeply with increasing 
mass flux. If instead conditions are such that radia- 
tion forces do not terminate the Keplerian flow out- 
side i?ms, the Keplerian flow is terminated by general 
relativistic corrections to Newtonian gravity and the 
radius of the transition to rapid radial inflow is ap- 
proximately independent of the mass flux. Hence, un- 
der these conditions the frequency of any QPO with 
the orbital frequency at the sonic point will be ap- 
proximately independent of the accretion rate. 

Although the numerical results presented here are 
for accretion rates typical of the atoll sources, which 
have very weak magnetic fields and are accreting at 
rates much less than the Eddington critical rate, we 
expect the structure of the accretion flow and the be- 
havior of the sonic-point Keplerian QPO frequency 
to be similar in the Z sources, even though they are 
thought to have stronger magnetic fields and are ac- 
creting at much higher rates. There are three reasons. 

(1) The radius of the transition to rapid radial in- 
flow that we are considering is necessarily close to 
the neutron star whatever the luminosity, because 
the fraction of the angular momentum in the Keple- 
rian disk flow that can be removed by radiation drag 
is small and general relativistic corrections to New- 
tonian gravity are important only near the neutron 
star. (Although the luminosities of the Z sources are 
much higher, the angular momentum flux that must 
be removed to create a radial inflow is correspond- 
ingly higher. In sources with luminosities very close 
to the Eddington critical luminosity, radiation drag 
can significantly affect the motion of up to ~ 20% of 
the accreting gas further away, but can only affect the 
motion of most of the accreting gas close to the star 
[see Fortner et al. 1989; Lamb 1989].) 

(2) The inferred magnetic fields of the neutron 
stars in the Z sources are stronger than those of the 
neutron stars in the atoll sources, and hence a larger 
fraction of the disk flow is likely to be channeled out 
of the disk by the stellar magnetic field, causing the 
mass flux through the inner disks of the Z sources to 
be comparable to the mass fluxes through the inner 
disks of the atoll sources. 



(3) The vertical thickness of the inner disk flow is 
likely to be larger in the Z sources than in the atoll 
sources. 

All of these effects act in the direction of making 
the radius of angular momentum loss in the Z sources 
similar to its value in the atoll sources. We do expect 
signiflcant departures from the behavior found in the 
present calculations when the luminosity equals or ex- 
ceeds the Eddington critical luminosity. These points 
are discussed in more detail in § 4. 

3.3. Frequencies of the Sonic-Point QPOs 

We consider now the frequencies that are generated 
by the sonic-point mechanism. There are two funda- 
mental QPO frequencies: (1) the Keplerian frequency 
at or near the sonic radius and (2) the sonic-point beat 
frequency, which is generated by interaction of radi- 
ation from the neutron star with the accretion flow 
near the sonic radius. QPOs may also be detectable 
at overtones of the Keplerian and beat frequencies. 
QPOs at sideband frequencies and oscillations at the 
spin frequency of the neutron star and its overtones 
are likely to be very weak. 

In analyzing the frequencies generated by the sonic- 
point mechanism, we shall for deflniteness use Boyer- 
Lindquist global coordinates, which are familiar be- 
cause they are commonly used to describe the space- 
time of rotating black holes. However, all of the ex- 
pressions we give are accurate only to first order in j 
and hence would be unchanged if we used circumfer- 
ential rather than Boycr-Lindquist coordinates. 

In Boyer-Lindquist coordinates, the circular fre- 
quency of an element of matter on the surface of the 
rotating star or of an element of gas in orbit around 
it is 

V EE {l/2TT){d4>ldt) , (9) 

where and t are the azimuthal and time coordi- 
nates of the element of matter or gas. When mea- 
sured in Boyer-Lindquist coordinates, any rotational 
frequency appears to be the same at every point in 
space and is equal to the frequency measured by a 
distant observer because of the symmetries of the 
spacetime outside a steadily rotating star^ and the 

''The spacetime outside a steadily rotating star has both timelike 
and spacelike Killing vector fields ^((j and C((^>)i which reflect 
the stationarity and axial symmetry of the spacetime (see, e.g., 
Misner, Thorne, & Wheeler 1973, pp. 892-895). In any co- 
ordinate system in which the time and azimuthal coordinates 
are based on ^((j and 5(0), such as the Boyer-Lindquist coordi- 
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fact that the Boyer-Lindquist time coordinate is the 
proper time of an observer at radial infinity. We em- 
phasize that the frequency (||) is not the frequency 
that would be measured by a local observer (i.e., in 
a local orthonormal tetrad), unless the observer is at 
infinity. 

Sonic-point Keplerian frequency. — As explained ear- 
lier in this section, the radiation drag force acting on 
the gas in a clump near the sonic radius generates a 
supersonic stream of denser gas that spirals inward 
toward the stellar surface. To the extent that the ra- 
diation field is axisymmetric, the inspiral trajectory 
of the gas from each such clump is identical and hence 
the azimuthal separation A^gtr between any two such 
streams is the same at every radius inside the sonic 
point and equal to the azimuthal separation A(j)ci of 
the two clumps at the sonic radius rg that are produc- 
ing the two streams. Therefore a collection of clumps 
distributed around the star at the sonic radius will 
generate a pattern of inspiral streams that rotates 
around the star at the orbital frequency vks = ^'k(''s) 
of the gas at the sonic radius. 

We can estimate vks by noting that, to first order 
in j, the the spacetime outside a slowly and uniformly 
rotating star is the same as the Kerr spacetime (Hartle 
& Thorne 1968; the two spacetimes are not the same 
to higher orders in j). Hence, to first order in j the 
orbital frequency of an element of gas in a prograde 
Keplerian circular orbit at Boyer-Lindquist radius rg 
in the rotation equator of a slowly rotating star is (see 
Bardeen, Press, & Teukolsky 1972) 
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The pattern of the gas streaming inward from the 
clumps at Tg rotates uniformly around the star with 
frequency vks- The pattern frequency is different 
from the orbital frequency i^orhir) of the gas inside 
the sonic radius, which varies with radius. In particu- 
lar, the pattern frequency is different from the orbital 



nate system, the time interval required for one rotation of an 
element of the star or one orbit of an element of gas orbiting 
the star is the same everywhere, when measured in the global 
time coordinate. Stated more concretely, if an element of gas 
emits a pulse of X-rays each time its azimuthal coordinate (f> 
increases by 27r, the time interval At between the arrival of 
successive pulses will be the same everywhere, when measured 
in the global time coordinate t. 



frequency of the gas just before it collides with the 
stellar surface. 

The brightness pattern made by the collision of 
the inspiraling streams of denser gas with the stel- 
lar surface rotates around the star at frequency vks, 
as measured by an observer at infinity. This is there- 
fore the frequency at which the radiation pattern pro- 
duced by the bright footprints of the streams rotates 
around the star and hence also the centroid frequency 
of the resulting quasi-periodic oscillation in the X-ray 
flux and spectrum seen by a distant observer whose 
line of sight is inclined to the rotation axis of the 
disk. Therefore the radiation pattern produced by 
the bright footprints of the streams generally does 
not corotate with the star, but instead moves around 
the stellar surface. 

Stellar spin frequency. — As explained in § 3.1, if 
the star's magnetic field is too weak to affect the ac- 
cretion flow at the sonic point but is strong enough to 
partially channel the flow close to the stellar surface, 
it will create a weakly beamed pattern of radiation 
that rotates with the star. The surface magnetic fields 
of the neutron stars in the atoll and Z sources, which 
are thought to be several hundred million years old, 
may be nearly dipolar. If the magnetic field is dipolar 
but offset from the center of the star, the radiation 
pattern it produces will have two unequal maxima 
around the star, whereas if the field is a centered but 
tilted dipole field, the radiation pattern will have two 
nearly equal maxima around the star. We therefore 
expect that there will generally be a brightness oscil- 
lation with a frequency equal to one or two times the 
spin frequency i^spin of the star. However, as we dis- 
cuss in detail in § 4.3, we expect the relatively weak 
magnetic fields of the atoll and Z sources to produce 
radiation patterns that are only weakly beamed, even 
near the star. When the substantial attenuation of 
this inherently weak beaming that is caused by scat- 
tering in the central corona is taken into account, the 
X-ray oscillation at Jospin or 2i^spin seen by a distant 
observer may be very weak or even undetectable with 
current instruments. 

Sonic-point beat frequency. — The weakly beamed 
pattern of radiation rotating with the star causes the 
radiation drag force acting on the gas in a given clump 
orbiting near the sonic radius to peak with a fre- 
quency equal to one or two times the difference (beat) 
frequency vbs between the Keplerian frequency at the 
sonic point and the spin frequency of the star, as ex- 
plained in § 3.1. The peak in the drag force will cause 
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a temporary increase in the flux of denser gas stream- 
ing supersonically inward from the clump with a fre- 
quency 



^spm 



(11) 



where fc=l or 2, depending on the symmetry (see 
Lamb et al. 1985). The inflow time from the sonic 
point is typically ~ 5-10 beat periods, so there are 
typically this many density enhancements in the stream 
of gas from a given clump to the stellar surface. These 
density enhancements do not lie along a single gas 
streamline, but instead along a sequence of stream- 
lines. As they collide with the stellar surface they 
produce a sequence of brighter impact arcs around the 
surface in the plane of the disk. As the quasi-periodic 
increase in the mass flux from a given clump arrives at 
the stellar surface, it produces a quasi-periodic mod- 
ulation in the luminosity of the star with frequency 
i^peak- Because this frequency is the difference be- 
tween an orbital and a spin frequency, the frequency 
of the resulting quasi-periodic modulation of the X- 
ray flux and spectrum seen by an observer at infinity 

IS l^pcak. 

Expected range of sonic-point QPO frequencies. — 
In § 3.2 we showed that either radiation forces or 
general relativistic corrections to Newtonian gravity 
cause an abrupt transition to rapid radial inflow and 
that the radius of this transition is bounded below 
by the radius of the marginally stable orbit, which is 
12 km for a nonrotating 1.4 M© star, and above by 
the ~ 20% upper bound on the fraction of the an- 
gular momentum of the accreting gas that can be 
transferred to the radiation field, which constrains 
the transition radius to be <30km for a 1.4 Mq 
star. As discussed in § 4.2, there are several effects 
that are likely to reduce further the allowed range 
of the transition radius, perhaps to ~ 15-25 km . 
Inserting these radii in equation ( |l0| ) gives an ex- 
pected range for the sonic-point Keplerian frequency 
of ^ 500-1200 Hz . For stars with spin frequencies 
~ 200-300 Hz , the expected range of beat frequen- 
cies is then ~ 200-1000 Hz . These ranges are similar 
to the frequencies of the kilohertz QPOs detected so 
far. 

Overtones and sidebands of the sonic-point QPO 
frequencies. — In addition to the beaming oscillations 
at i^Ks and i^spin and the luminosity oscillation at 
t'Bs, we expect accretion onto weakly magnetic neu- 
tron stars to generate brightness oscillations at several 
other frequencies. For example, the X-ray brightness 



modulation produced by the motion of the brighter 
impact footprints around the star with frequency j/ks 
is unlikely to be exactly sinusoidal, so there may be 
some power at overtones of t'Ks- However, the angu- 
lar distribution of the X-ray emission from the impact 
footprints is expected to be very broad (see § 3.1), so 
the power at overtones of z^ks is likely to be much 
smaller than the power at • The radiation pattern 
that rotates with the star is also expected to be very 
broad by the time the radiation escapes from the cen- 
tral corona, so the power at overtones of Vspin is likely 
to be much smaller than the power at Jospin, which is 
itself expected to be small for the reasons discussed 
above. However, the radiation pattern that rotates 
with the star is likely to be narrower at the sonic 
radius than outside the central corona, causing the 
modulation of the mass flux from the sonic point to 
be somewhat non-sinusoidal and giving rise to weak 
overtones of the sonic-point beat frequency. Lumi- 
nosity oscillations are not as strongly attenuated by 
scattering in the central corona (see § 4.1), so it may 
be easier to detect overtones of v^s- 

The sonic-point mechanism causes a modulation of 
the inward mass flux and hence the stellar luminos- 
ity at the sonic-point beat frequency; from the point 
of view of a distant observer, this luminosity is then 
modulated at the sonic-point Keplerian frequency by 
the motion of the brighter footprints around the star. 
As a result, QPOs are generated at some sideband 
frequencies. However, by far the strongest QPOs are 
generated at the sonic-point Keplerian and beat fre- 
quencies. To see this, consider as an example the 
X-ray flux waveform 



f{t) = [l + KiC0s{2TriyKt + (f>Ki) 

+ K2 cos(47ri/Ki -t- 0K2)] X 
X [1 + Bi cos(27rt'Bi -I- 0bi) 

+ B2 COs{4TTUBt + (t)B2)] 



(12) 



which describes these modulations. In equation (12) 
we have included for the sake of illustration only the 
fundamental and flrst overtones of the Keplerian and 
beat frequencies; higher harmonics may of course also 
be present. The frequencies that are generated are 
listed in Table ||. The only flrst-order QPOs are 
those at the sonic-point beat and Keplerian frequen- 
cies. There are also QPOs at the second harmonics 
of the beat and Keplerian frequencies. The lowest- 
order QPOs generated by the modulation of the beat- 
frequency "carrier" by the motion of the footprints 
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around the star are of second order and are at the 
stellar spin frequency I'spin and at the very high fre- 
quency 2i/Ks — Jospin- The only QPO of intermedi- 
ate frequency generated by this modulation is a QPO 
at vks — 2i/spin that is third-order and hence likely 
to be very weak. All the other QPOs generated are 
also third-order and in addition have very high fre- 
quencies. They are therefore likely to be intrinsically 
weak, strongly attenuated by scattering in the central 
corona, and difficult to detect. 

As we discuss in § 4.3, only brightness oscillations 
that are fairly large near the neutron star and are 
not strongly attenuated by the gas surrounding the 
star will be detectable by a distant observer. In 
the sonic-point model, only the sonic-point beat fre- 
quency QPO and its overtones are luminosity oscil- 
lations; the other QPOs are beaming oscillations and 
are therefore more attenuated by scattering. For this 
reason, we expect strong QPOs only at i^bs and vks- 
Nonetheless, QPOs at other frequencies, such as I'spin, 
2i'Bs, and 2z^ks may be detectable. Detection of oscil- 
lations at any of these frequencies would corroborate 
the sonic-point model. 

3.4. Coherence of Sonic-Point QPOs 

A key question any theory of the kilohertz QPO 
sources that relates a QPO frequency to orbital mo- 
tion must address is why narrow peaks are seen in 
power spectra of the brightness variations of these 
sources, rather than a broad continuum correspond- 
ing to the range of orbital frequencies in the in- 
ner disk. Indeed, the remarkably high coherence 
(i^QPo/<5^QPO ^ 100) of some of the kilohertz QPOs 
places very strong constraints on any model of these 
QPOs, because there are many physical effects that 
tend to decrease the coherence of oscillations at kilo- 
hertz frequencies. We discuss some of the most im- 
portant effects, derive the resulting constraints on the 
accretion flow, and show that the sonic-point mech- 
anism can produce narrow peaks consistent with the 
observed coherence of the kilohertz QPOs. 

Gas in the inner disk orbits the neutron star with a 
wide range of frequencies, so at first glance one might 
think that fluctuations in the gas density throughout 
the inner disk would generate a broad spectrum of 
brightness variations. For the reasons discussed in 
§ 2.3, it is very unlikely that such fluctuations can 
produce large amplitude brightness variations via di- 
rect emission or occultation. We have therefore fo- 
cused our attention on the brightness variations that 



such density fluctuations produce indirectly, as gas 
from them falls inward and impacts the stellar sur- 
face. Here we show that the accretion flow in the 
vicinity of the sonic radius acts as a filter that se- 
lects brightness oscillations at the sonic-point orbital 
frequency while suppressing brightness oscillations at 
higher and lower frequencies. In order to see how this 
filtering occurs, we analyze the X-ray emission at the 
stellar surface produced by clumps that form at dif- 
ferent radii in the inner accretion disk. 

Effects of azimuthal shear. — Let us suppose for the 
sake of argument that the thickness of the inner disk 
is infinitesimal and that the flow there is laminar (we 
discuss the effects of turbulence shortly). Suppose 
further that a very small, roughly spherical density 
fluctuation (clump) has formed well outside the sonic 
radius and consider what will happen to the gas in this 
clump and how it will affect X-ray emission from the 
star. When a small clump forms, the orbital phases 
of the elements of gas that comprise it are necessarily 
very similar, because of the small azimuthal extent 
of the clump. The elements of gas that comprise the 
clump also have very similar orbital frequencies, be- 
cause the radial extent of the clump is small. As the 
clump drifts inward, elements of gas at different radii 
are sheared relative to one another in the azimuthal 
direction. 

The frequency of the X-ray brightness oscillation 
that is generated by the inspiraling gas from a clump 
is approximately equal to the Keplerian frequency at 
the radius where the clump originally formed (see 
the discussion in § 3.3) A clump that forms very 
near the sonic radius will therefore produce an X-ray 
brightness oscillation at the sonic-point Keplerian fre- 
quency. If clumps can form inside the sonic radius, 
such clumps would produce X-ray brightness oscilla- 

*At first glance, one might think that if the inward radial drift 
time from a given radius r in the inner disk to the stellar surface 
is large, azimuthal shearing of a clump formed at r would by 
itself reduce the amplitude of the brightness variation produced 
by gas from the clump, but this is not so. The reason is that 
to first order, the azimuthal velocity shear in a Keplerian flow 
conserves the volume of the clump, because the velocity field is 
divergence-free to this order. Hence, the density of the gas in a 
clump relative to the density of the background gas is affected 
only weakly by the differential rotation. In fact, if the flow 
were incompressible as well as laminar, the angular width of 
the brighter impact arc on the stellar surface would also be 
unaffected, even if the total angular extent of the clump at all 
radii becomes large, because the angular extent of that part of 
the clump that is crossing the sonic radius at any given time is 
unchanged by azimuthal shear. 
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tions at the orbital frequencies where they formed. 

This analysis shows that if the clumps were small 
and the flow laminar, differential rotation of the gas 
would not by itself pick out any particular orbital fre- 
quency. Consequently, if clumps were to form at a 
wide range of radii and nothing besides inward drift 
and azimuthal shear were to happen, gas inspiraling 
from the clumps and colliding with the stellar surface 
would produce brightness fluctuations with a wide 
range of frequencies up to the ~ 1500-2000 Hz orbital 
frequency at the stellar surface (see, e.g., Kluzniak, 
Michelson, & Wagoner 1990). 

Effects of dump destruction. — In reality, clumps 
that form outside the sonic radius arc destroyed be- 
fore they drift inward to the sonic point and gas in 
them can reach the stellar surface. Clumps (if any) 
that form inside the sonic radius orbit the star at a 
rapidly changing frequency and collide with the stel- 
lar surface before gas from them can produce a long 
wavetrain. Hence, only clumps that form very near 
the sonic radius can produce strong, relatively coher- 
ent brightness oscillations. In this way the disk flow 
filters out brightness oscillations at orbital frequencies 
other than the sonic-point Keplerian frequency. 

In order to illustrate the effects of clump destruc- 
tion in the Keplerian flow, wc consider a simple model 
in which clumps form and are then destroyed by 
small-scale turbulence in the disk, such as that pro- 
duced by the magnctoturbulcncc that is thought to be 
responsible for angular momentum transport within 
the disk flow (Balbus & Hawley 1991, 1992, 1997; 
Brandenburg ct al. 1995, 1996; Hawley, Gammic, & 
Balbus 1995, 1996). In this simple model we assume 
that the forces that formed the clump in the first place 
no longer act to hold it together. If the clump has 
size ^, turbulent motions on scales will tend to 
advect the clump without destroying it; only turbu- 
lent motions on scales i' <C ^ will disrupt the clump. 
We assume that the effect of these small-scale tur- 
bulent motions on the clump can be described by a 
diffusion coefficient D ^ /?^f;,s, where (3 is a dimcn- 
sionless parameter that describes the strength of the 
small-scale turbulence and Cg is the thermal sound 
speed in the disk; for turbulent motions with scales 
^ -C ^, /3 is small compared to unity. As a result of 
outward angular momentum transport by the magnc- 
toturbulcncc, a clump of initial size ^ that forms at 
radius r will drift slowly inward with radial velocity 
v'^ a{h/r)cs, where a is the usual viscosity pa- 
rameter and h is the half-thickness of the disk. Such 



chimps will be dissipated by turbulent diffusion in a 
time tdiss ~ (^//3c,) ~ iC/PhnK). 

In this model, a clump that forms a radial distance 
Sr outside the sonic radius can reach the sonic ra- 
dius and generate a supersonically inspiraling stream 
of gas before being destroyed only if Sr < v'^t^iss, i-e., 
only if the clump forms within a radial distance 

Sr < {a/p){h/r,){^/r,)r, < lO'^ r, (13) 

of the sonic radius. Turbulence on small scales may be 
weak, and the effects that form the clumps in the first 
place — such as thermal instability, magnetic stresses, 
and radiation forces — may tend to continue to hold 
them together as they drift inward toward the sonic 
point, so that clumps can survive somewhat longer 
than tdiss- Even so, it is clear that clumps that form 
even a small distance outside the sonic radius will be 
disrupted by turbulence before they reach the sonic 
point. This process suppresses brightness oscillations 
with frequencies less than the sonic-point Keplerian 
frequency. 

It appears improbable that clumps will form in the 
hypersonic radial inflow between the sonic radius and 
the stellar surface, but even if clumps do form in this 
region, they will not produce appreciable brightness 
oscillations, for two reasons. First, a clump that forms 
inside the sonic radius will collide with the stellar sur- 
face after completing ~1-10 orbits. A clump that 
lives only for a time tufetimc will generate power over a 
range of frequencies (5fiifctimo ^ (Trtufctime)"^- Hence 
the power produced by a clump that forms inside the 
sonic radius will be spread over a frequency range 
<5j^iifctimc <^ 0.03 v. Second, a clump that forms inside 
the sonic radius orbits the star at a frequency that 
rapidly increases from, for example, the ~ 500 Hz or- 
bital frequency at the sonic radius to the ^ 1500 Hz 
orbital frequency at the stellar surface, as the clump 
spirals inward. For these reasons, even if clumps do 
form inside the sonic radius, any power they generate 
will be spread over a wide range of frequencies rather 
than concentrated in a narrow peak. 

Coherence of the sonic-point Keplerian frequency 
QPO. — Clumps that form outside the sonic radius but 
close enough to it to reach it before being destroyed 
have initial orbital frequencies in the relatively narrow 
range 




<10-2i/Ks, (14) 
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which is consistent with the observed narrowness of 
the kilohertz QPO peaks. Clumps orbiting at dif- 
ferent distances above or below the midplane of the 
disk will generally also have slightly different orbital 
frequencies. To estimate the resulting spread in fre- 
quencies, we assume that at each distance z above 
and below the disk plane the inward radial veloc- 
ity increases sharply from subsonic to supersonic at 
some cylindrical radius R^. We call the axisymmet- 
ric two-dimensional surface defined by Rs{z) the sonic 
surface. If the sonic surface is cylindrical, the oscil- 
lation frequency generated by the clumps orbiting in 
the midplane of the disk is slightly greater than the 
oscillation frequency generated by the clumps orbit- 
ing above and below the midplane. The spread in 
oscillation frequencies depends on the precise shape 
of the sonic surface, but to lowest order in {h/vs) 
the spread in frequencies caused by this effect is 
Siy ^ (h/rg)'^ I'Ks) which is consistent with the ob- 
served narrowness of the kilohertz QPO peaks if 
hjr, <0.1. 

The persistence of a clump orbiting at the sonic 
radius is limited by its decay as gas inspirals from it 
to the stellar surface and by its disruption by small- 
scale turbuk;nc;e, pressure forces, and other effects. 
In the simple model considered here, the clump decay 
timescale is idccay ^ S./v'' ~ {l/a){^/hs){rs/hs){l/ilK), 
which broadens the QPO peak by an amount SiMccny ^ 
ct{hs/^){hs/rs) fKs ^ 10~^ vks, which is consistent with 
the observed narrowness of the kilohertz QPO peaks. 
The dissipation of clumps by turbulent diffusion con- 
tributes a relative width ^z^diss ~ l3{hs/^) pks, which is 
consistent with the observed narrowness of the Keple- 
rian frequency QPO peaks if /3 <C 1 or the effects that 
form the clumps in the first place, such as thermal in- 
stability, magnetic stresses, and radiation forces, con- 
tinue to hold them together as they orbit at the sonic 
point. 

This simple model shows that if a clump is too 
small, it will persist for such a short time that the 
power it generates will be sprc^ad over a broad fre- 
quency range, whereas if a clump is too large it will 
persist for such a long time that its orbital frequency 
will change appreciably during its lifetime and the 
power it generates will again be spread over a broad 
frequency range. It is the clumps of intermediate size 
that generate relatively coherent oscillations. 

In addition to frequency variation and lifetime 
broadening, other effects may increase the width of 
the Keplerian frequency QPO peak. For example, the 



radiation field inside the sonic radius is not perfectly 
axisymmetric, and hence the radiation forces acting 
on the inspiraling gas will vary slightly with azimuth. 
As a result, the radius of the sonic point and the den- 
sity patterns produced by the inspiraling gas will be 
slightly different at different azimuths. A quantita- 
tive treatment of these effects is well beyond what is 
possible with current accretion disk models. 

The simple model we have discussed shows the cru- 
cial role played by the sharp transition to rapid ra- 
dial inflow at the sonic radius that was found in § 3.2. 
In the absence of such a transition, density and mag- 
netic field fiuctuations in the inner disk would produce 
weak brightness variations over a wide range of fre- 
quencies rather than in the narrow range of frequen- 
cies needed to create a QPO peak. In the presence 
of such a transition, on the other hand, clumps orbit- 
ing in a narrow range of radii near the sonic radius 
produce strong X-ray brightness variations with a cor- 
respondingly narrow range of frequencies. Simple es- 
timates of the broadening produced by the different 
orbital frequencies of those chimps and by their finite 
lifetimes appear consistent with the high observed co- 
herence of the Keplerian frequency QPOs. 

Coherence of related QPOs. — If there is a QPO at 
the beat frequency i^bs , we expect that its FWHM will 
in general be comparable to the FWHM of the QPO 
at vks- This is because the peak at vbs is at the beat 
frequency of the sonic-point Keplerian frequency with 
the stellar spin frequency, which is nearly coherent 
and therefore adds relatively little to the width of the 
QPO peak. However, we do not expect the widths of 
the QPO peaks at ly^s and vks to be identical, because 
scattering by the moving gas in the central corona 
affects the oscillations at vbs and z^ks differently. 

3.5. Photon-Energy Dependence of the Sonic- 
Point QPOs and Time Lags 

The X-ray spectra of systems powered by accret- 
ing, weakly magnetic neutron stars, such as the kilo- 
hertz QPO sources, are formed primarily by Comp- 
tonization, as described in § 2.2. Cyclotron and 
bremsstrahlung photons with energies < 1 keV are 
produced near the neutron star surface and are then 
Comptonized by hot electrons in the region where the 
accretion flow interacts with the stellar surface and in 
the hot central corona that surrounds the star and its 
magnetosphere, yielding the observed X-ray spectra. 

During a sonic-point Keplerian or beat-frequency 
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oscillation, the rotation of the patterns of denser gas 
spiraling inward from clumps orbiting at the sonic ra- 
dius causes the optical depth along the line of sight 
from the neutron star surface to the observer to vary 
quasi-periodically. The rate of production of soft pho- 
tons and the physical properties of the Comptonizing 
gas may also oscillate, especially during the sonic- 
point beat-frequency oscillation, which is mainly a 
luminosity (accretion rate) oscillation. These oscil- 
lations cause the X-ray spectrum produced by the 
system to oscillate quasi-periodically with various fre- 
quencies. A small fractional variation in the opti- 
cal depth causes a much larger fractional variation in 
the number of photons at high photon energies, be- 
cause of the characteristic way in which the spectrum 
produced by Comptonization changes as the optical 
depth oscillates (see Miller & Lamb 1992). To the 
extent that the effects of Comptonization dominate, 
the oscillations at > 10 keV will lag the oscillations 
at ~ 5 keV , because the photons above 10 keV have 
scattered more times in escaping from the source. 

An accurate, quantitative treatment of the X-ray 
spectral oscillations produced by the inhomogeneous 
and time-dependent accretion flow in the kilohertz 
QPO sources would require a three-dimensional, time- 
dependent radiation hydrodynamic calculation in full 
general relativity. This is beyond present computa- 
tional abilities. However, one can obtain a qualitative 
understanding of the dependence of QPO amplitudes 
on photon energy and the time lags to be expected us- 
ing the following simplified model of the Comptoniza- 
tion process. 

Assume that soft photons with a characteristic en- 
ergy Ein that is much smaller than the electron tem- 
perature Te are injected at the center of a static, uni- 
form, spherical Comptonizing region of radius Rc, 
electron density rig, and scattering optical depth r = 
neCT-Rc: where ctt is the Thomson scattering cross 
section. Assume further that y = ATeT^/rrie <^ 1, 
where mg is the electron rest mass, and that the elec- 
trons have a negligible bulk velocity. These are good 
approximations for modeling formation of the time- 
averaged X-ray spectra of the atoll sources (see Psaltis 
& Lamb 1998b). The region of approximately radial 
inflow that develops when the mass accretion rate be- 
comes comparable to Me (see Lamb 1989, Lamb 1991, 
and § 2) introduces additional effects on the spectra 
of the Z sources that cannot be treated in this way; 
we do not consider these effects here (but see Psaltis 
& Lamb 1998b). We treat the effects of oscillations in 



the rate of soft-photon production and in the prop- 
erties of the Comptonizing region near the star by 
varying these quantities in the spherical model. In 
particular, we mimic the effects of the quasi-periodic 
variation of the optical depth along the line of sight by 
varying the optical depth of the model. We assume 
that the oscillation period is longer than the mean 
time for photons to escape from the Comptonizing re- 
gion, an approximation that appears to be excellent 
for the kilohertz QPO sources. 

Photon- energy dependence of the sonic-point Ke- 
plerian and heat-frequency oscillation amplitudes. — 
The X-ray number spectrum that emerges from the 
Comptonizing region at photon energies E 3> -Bin can 
be approximated by (see, e.g., Rybicki & Lightman 
1979, pp. 221-222) 



(^ 



2+a 



E<^Eb 

(15) 

where Tg is the electron temperature, Ea is a nor- 
malization constant that is typically a few keV, Eb ~ 
(l-3)Te is the cutoff energy, and 
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For small-amplitude oscillations in which the varia- 
tion of the electron temperature in the Comptonizing 
region is negligible (see below), the relative ampli- 
tude 7(i?) of the oscillation is approximately (Miller 
& Lamb 1992) 
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dr 
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Here ANg/Ns is the relative change in the rate at 
which soft photons are injected into the Comptoniz- 
ing region during an oscillation and At is the change 
in the optical depth of the region. The relative ampli- 
tude j{E) is computed by considering a full oscillation 
period. 

In general, the soft photon injection rate and the 
optical depth both vary during the oscillation. If so, 
equations (p^-dT7|) give 
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where 
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in terms of the expected properties of the Comptoniz- 
ing region and we have neglected the weak depen- 
dence of the normahzation constant Ea on r. Equa- 
tions (^8|) and (|l^) show that in this case the relative 
amplitude of the QPO depends both on the relative 
change in the soft-photon injection rate, ANg/Ng, 
and on the relative change in the optical depth, At/t. 
In writing equation (jl^) we have assumed that the 
electron temperature in the Comptonizing cloud re- 
mains constant during the oscillation. If the electron 
temperature were allowed to vary, then the relative 
amplitude of the oscillation would not be constant at 
high photon energies, in disagreement with observa- 
tions. 

The quantity 7,- is negative for the spectrum ( |l5| ) 
and is small at low photon energies. If the oscilla- 
tion in the injection rate of the soft photons is in 
phase with the oscillation in the optical depth, i.e., 
if ANs/Ns and At/t have the same sign, and if 7,- 
is larger in magnitude than ANg/Ns at high photon 
energies, then the relative amplitude of the oscillation 
in the photon number will not be a monotonic func- 
tion of photon energy. For ANg/Ng comparable in 
magnitude to At/t and a Comptonizing region with 
the properties used in scahng equation (pj[), the rel- 
ative amplitude of the oscillation has a minimum in 
the energy range ~ 5-10 keV . The sonic-point beat 
frequency QPO is primarily a luminosity oscillation 
(see § 3.1), so the photon production rate is likely to 
vary significantly at the beat frequency. The relative 
amplitude of the beat-frequency QPO may therefore 
have a minimum in the 5-10 keV energy range. 

If the relative change in the photon injection rate 
during an oscillation is very small, the first term on 
the right side of equation (|l^) is likely to be negligible 
compared to the second term. In this case the rela- 
tive amplitude of the oscillation will increase mono- 
tonically with increasing photon energy, for energies 
less than i?b ~ 10 — 30 keV, but will become indepen- 
dent of photon energy for energies greater than Eb- 
We expect this to be the case for the oscillation at 
the sonic-point Kcplcrian frequency, which is mainly 



a beaming oscillation (see § 3.1) and hence involves 
a relatively small variation in the photon production 
rate. We therefore expect the relative amplitude of 
the Keplerian- frequency QPO to increase steeply in 
the 5-10 keV energy range and then flatten at high 
energies. 

When the photon injection rate and spectrum of 
injected soft photons do not change appreciably dur- 
ing a QPO cycle, the amplitude of the oscillation in 
luminosity caused by the oscillation in optical depth 
is comparable to the amplitude of the optical depth 
oscillation and is given by 



^^A, = 0.07f^ 
L ^ I 10 kcV 
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Equations ( |19| ) and (|2C|) demonstrate that a relatively 
weak oscillation in optical depth of only a few per- 
cent that is accompanied by a luminosity oscillation 
of moderate amplitude can produce naturally a much 
larger oscillation in the countrate at high photon en- 
ergies. 

Equation (|l^) also shows that the relative ampli- 
tude of the oscillation at a given photon energy is 
different for different electron temperatures and op- 
tical depths, unless the relative amplitude of the os- 
cillation in optical depth At/t varies in such a way 
as to compensate for this effect. We therefore ex- 
pect that as the mass accretion rate onto the neutron 
star changes (on timescales much longer than a beat- 
frequency period), the electron temperature and op- 
tical depth in the Comptonizing region will change, 
causing the photon-energy dependence of the beat- 
frequency and Keplerian-frequency oscillation ampli- 
tudes to change. 

Photon- energy dependence of the sonic-point Kep- 
lerian and beat-frequency oscillation phases. — At pho- 
ton energies < 30 keV, the cross section and hence the 
electron scattering mean free path of the photons in 
the central corona is nearly independent of photon en- 
ergy. Therefore, the average escape time of photons 
from the corona is also independent of their energy. 
However, the photons that stay in the corona longer 
experience more scatterings on average, and therefore 
gain more energy by scattering off hot electrons and 
emerge from the medium with a larger energy. As a 
result, the oscillation at high photon energies is ex- 
pected to lag the oscillation at low photon energies. If 
the injection rate of soft photons and the properties 
of the corona are constant in time, the magnitude 
of the time lag is determined mainly by the prop- 
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erties of the corona (see, e.g., Wijers, van Paradijs, 
& Lewin 1987; Bussard et al. 1988). If instead the 
source of soft photons or the properties of the corona 
are time-dependent, then the time lag wiU depend on 
the details of this time-dependence. 

For the simpUfied model described earlier in this 
section, the average energy after u scatterings of a 
photon injected at energy E-m is 



E(u) ~ min[£'in exp(4Teu/me), Eb 



(21) 



Therefore, for photons with energies smaller than 
~ Ef,, the ratio of the energies of two photons that 
have experienced a different number of scatterings is 
E2/E1 ^ exp{ATeAu/me), where Au is the difference 
in the number of scatterings. The photon mean- free 
time is ~ (rieCTc)^^ and hence the time lag intro- 
duced between the oscillations at the two photon en- 
ergies is 



St ^^I!^ln ^ 
crATe \Ei 
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400 
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when both Ei and E2 are smaller than ~ Et- Be- 
cause of the diffusion in energy and the systematic 
downscattering of photons with energies > Et, equa- 
tion (^2|) is not valid at energies ^ Et- 

We emphasize that, for a variety of reasons, equa- 
tion (^) gives only an upper bound to the expected 
time lag at photon energies <, Eb- First, this equation 
assumes that all photons are injected at the center of 
a spherical medium. However, if most of the photons 
are produced near the surface of the neutron star and 
escape a few kilometers above it, the time lag can be 
significantly smaller. Second, equation (p^ ) was de- 
rived under the assumption that photons slowly dif- 
fuse outward. This is not a very good approximation 
for electron scattering optical depths t < 3; if a frac- 
tion of the photons escape directly from the corona 
without interacting with the electrons, the average 
escape time could be significantly smaller. Finally, 
equation (|2^) was derived under the assumption of 
a uniform electron density. If instead most of the 
optical depth is concentrated near the center of the 
corona, this will also reduce significantly the average 
escape time. 

In summary, to the extent that the effects of Comp- 
tonization dominate, we expect the sonic-point Keplerian- 
and heat-frequency oscillations at high photon ener- 
gies should lag the corresponding oscillations at lower 



photon energies and the time lag between a few keV 
and ~ 15 keV should be a fraction of a millisecond, 
for the densities and temperatures of the electrons 
expected around the neutron stars in the Z and atoll 
sources. 

3.6. Attenuation of Kilohertz QPOs 

In the sonic-point model, gas spirals inward from 
density inhomogeneities in the accretion flow at the 
sonic radius and collides with the stellar surface, pro- 
ducing a radiation pattern that rotates around the 
star with the Keplerian frequency at or near the sonic 
point vks- Because the rate of mass accretion over 
the whole neutron star surface remains approximately 
constant in time as the radiation pattern rotates, the 
total luminosity emerging from the system also re- 
mains approximately constant in time. For concise- 
ness we call the oscillations produced by rotation of 
this radiation pattern a beaming oscillation, because 
it is caused by the angular variation of the radiation 
field, even though the radiation pattern is unlikely to 
be a narrow beam. In a pure beaming oscillation, the 
total luminosity emerging from the source is indepen- 
dent of time. As described earlier, the accretion flow 
is also expected to produce a much weaker beaming 
oscillation at the stellar spin frequency i^spin- 

If the magnetic field of the neutron star is strong 
enough to channel the accretion flow near the stellar 
surface, the sonic-point model predicts that the re- 
sulting radiation pattern, which rotates with the star, 
will modulate inflow from the inner disk at the beat 
frequency i'ks — t^spim producing a luminosity oscil- 
lation with this frequency. In a pure luminosity os- 
cillation, the total luminosity of the system changes 
with time, but the angular pattern of the radiation 
field remains unchanged. In the sonic-point model, 
the only luminosity oscillation is the beat frequency 
oscillation.^ 

When the luminosity oscillation is present, it is 
necessarily modulated by the motion of the emis- 



^To clarify further the difference between luminosity and beam- 
ing oscillations, consider an observer whose line of sight to the 
neutron star is along the rotation axes of the disk and the 
star. X-rays from beams rotating with the gas in the disk or 
with the star will appear time-independent to such an observer 
and therefore will not generate any X-ray flux or color oscilla- 
tions. However, the luminosity oscillation at the beat frequency 
caused by the oscillation of the mass flux at this frequency will 
be visible. Such an observer will therefore see an oscillation 
at the beat frequency and its overtones, but not at any other 
frequency. 



29 



sion regions around the stellar surface, which gen- 
erates second-order beaming oscillations at the stel- 
lar spin frequency I'spin and the difference frequency 
2i/Ks — t'spin (see § 3.3 and Table 3). In general, higher- 
order oscillations are also generated at other (mostly 
much higher) frequencies. 

Scattering generally attenuates beaming oscilla- 
tions much more than luminosity oscillations, for the 
conditions relevant to LMXBs (Lamb 1988). The rea- 
son is that an anisotropic radiation pattern becomes 
approximately isotropic after only a few scatterings, 
whereas a luminosity oscillation is attenuated by scat- 
tering only if the time required for photons to escape 
from the scattering region is much larger than the 
period of the oscillation, which requires a large opti- 
cal depth for the oscillation frequencies and coronal 
dimensions of interest to us. 

The relative amplitude Aoo of a luminosity oscilla- 
tion with angular frequency uj = Itxv outside a spher- 
ical scattering cloud of radius i?c and optical depth r 
is related to the relative amplitude at the center of 
the cloud by the expression (Kylafis & Phinney 1989) 
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where x = [(3/2)u;i?c'''/c]^/^. This expression is accu- 
rate to better than 6% for a; > 1.3 (Kylafis & Phinney 
1989). The first term in parenthesis on the right side 
of equation (^3|) describes the amplitude attenuation 
caused by the spread of escape times from the cloud 
whereas the second term (which has been added in 
by hand) describes the amplitude of the oscillation 
produced by the photons that escape from the cloud 
without scattering. 

For a beaming oscillation caused by rotation of a 
narrow pencil beam, the amplitude of the oscillation 
outside the cloud is (Brainerd & Lamb 1987; Kylafis 
& Phinney 1989) 
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The factor multiplying Ao^boam in equation ( p^ de- 
scribes the tendency of scattering to isotropize the 
photon distribution. The amplitudes of oscillations 
produced by the broad radiation patterns or the ra- 
diation patterns with more than one lobe that we are 
concerned with here are reduced even more by scat- 
tering (Brainerd & Lamb 1987). 

Of all the beaming oscillations that may be gener- 
ated by the accretion flow, only the one at the sonic- 



point Keplerian frequency is likely to have a large 
enough amplitude at the neutron star to produce an 
oscillation that is strong enough, after attenuation by 
passage through the central corona, to be observed 
easily with current instruments. In contrast, the os- 
cillation at the sonic-point beat frequency, which is 
the only luminosity oscillation produced by the flow, 
should be strong enough to be observed even if its am- 
plitude at the neutron star surface is moderate. As a 
result, we expect that from the collection of possible 
frequencies of oscillations, the only oscillations that 
will appear strong far from the star will be those at 
the sonic-point Keplerian and beat frequencies. 

4. COMPARISON WITH OBSERVATIONS 

In this section we describe the observational impli- 
cations of the sonic-point model and compare them 
with observations of the kilohertz QPOs. We first 
list the important properties of these QPOs that any 
model must explain. These properties are (see § 1): 

(1) high (~ 300-1200 Hz) frequencies, which can vary 
by several hundred Hertz in a few hundred seconds 
(see Wijnands et al. 1998 and van der Klis 1995); 

(2) similar frequency ranges in stars with significantly 
different magnetic field strengths and accretion rates; 

(3) relatively high coherence (Q = vjbv up to ^ 100); 

(4) 2-60 keV rms amplitudes < 1% in the Z sources 
but much higher (up to ^ 15%) in the atoll sources; 

(5) the common occurrence of two (but never more 
than two) simultaneous kilohertz QPOs in a given 
source; (6) a frequency separation IS.v between the 
two kilohertz QPOs that is approximately constant 
in a given source; (7) approximate consistency of A;/ 
with the stellar spin frequency inferred from burst os- 
cillations; (8) the similarity of most atoll and Z source 
spin frequencies; (9) the frequently similar FWHM of 
the two QPO peaks seen simultaneously; (10) the in- 
crease of the frequencies of the kilohertz QPOs with 
increasing inferred accretion rate that is observed in 
many sources; and (11) a steep increase in the am- 
plitude of the higher-frequency of the two kilohertz 
QPOs with increasing photon energy in the 2-15 keV 
energy band in many sources. 

In § 4.1 we show that the stellar magnetic field 
strengths and electron scattering optical depths re- 
quired for the sonic-point mechanism to operate are 
consistent with the stellar magnetic field strengths 
and optical depths of the compact central corona in- 
ferred from previous modeling of the X-ray spectra 
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and < 100 Hz X-ray variability of the atoll and Z 
sources. In § 4.2 we compare the kilohertz QPO and 
neutron star spin frequencies expected in the sonic- 
point model with the kilohertz QPO and neutron star 
spin frequencies observed in the atoll and Z sources. 
We also discuss the expected and observed depen- 
dence of the QPO frequencies on accretion rate. In 
§ 4.3 describe the amplitude and coherence of the vari- 
ous QPOs expected in the sonic-point model, compar- 
ing them with the observed amplitudes and coherence 
of the kilohertz QPOs. In particular, we show that 
the sonic-point model explains naturally why at most 
two kilohertz QPOs have been detected in the power 
spectrum of each source. We also discuss the relative 
coherence of QPOs at different frequencies and the 
expected dependence of QPO amplitudes on accre- 
tion rate and photon energy. In § 4.4 we describe the 
inverse correlation between QPO amplitude and stel- 
lar magnetic field expected in the sonic-point model 
and show that there is already substantial evidence 
for such a correlation. 

4.1. Consistency virith Previously Inferred 
Properties of the Atoll and Z Sources 

The operation of the sonic-point mechanism for 

producing strong, coherent kilohertz QPOs requires 
that a number of conditions be satisfied: 

(1) In order to produce the higher- frequency kilo- 
hertz QPOs with the large amplitudes observed in 
some sources, the neutron star magnetic fields in these 
sources must be relatively weak, so that a substantial 
fraction of the accreting gas penetrates close to the 
star in a Keplerian disk flow (§ 3.1). 

(2) In order to produce the lower-frequency kilo- 
hertz QPOs, the neutron star magnetic field, although 
relatively weak, must be strong enough to channel 
some of the accreting gas near the stellar surface to 
produce bright spots that rotate with the star and 
modulate the accretion rate at the sonic-point beat 
frequency (§ 3.1). 

(3) In order to generate QPOs with the high co- 
herence observed in most sources, the disk flow at the 
sonic point must be geometrically thin (§ 3.4). 

(4) In order that only the QPOs at the sonic-point 

Keplerian and beat frequencies be detectable and that 
X-ray oscillations at the stellar spin frequency be very 
weak or undetectable at present sensitivities, electron 
scattering in the central corona must help to atten- 
uate the weak beaming oscillations at other frequen- 



cies and at the stcillar spin frequency. The central 
corona must therefore have an electron scattering op- 
tical depth >3 (§ 3.6). 

As we now explain briefly, all of these conditions 
follow naturally from the unified model of weak-field 
accreting neutron stars described in § 2.2. 

The magnetic field strengths in the 4U atoll sources 
are thought to be low enough ( < 10^ G) that at their 
inferred accretion rates (~ 0.01-0.1 Me) the cylindri- 
cal radius wq at which the Keplerian disk flow couples 
strongly to the stellar magnetic field is < 2 x 10^ cm. 
The magnetic field strengths of the neutron stars in 
the Z sources are thought to be a few times larger 
than in the 4U atoll sources, but the accretion rates 
are much larger as well, so wq is comparable to its 
value in the atoll sources. Thus, in both the 4U atoll 
sources and the Z sources, the gas in the Keplerian 
disk penetrates close to the star before any of it cou- 
ples strongly to the stellar magnetic field. We there- 
fore expect that a significant fraction of the accreting 
gas will remain in a Keplerian disk flow down to the 
sonic point and will continue in a disk flow very close 
to the stellar surface, as required in the sonic-point 
model. 

Evidence that the magnetic fields of neutron stars 
in LMXBs, while typically weak, are nevertheless 
strong enough in many sources to channel the flow 
near the star and hence to produce a QPO at the 
sonic-point beat frequency comes from power spec- 
tra constructed from RXTE observations of GX 5—1 
(van der Klis et al. 1996e), Sco X-1 (van der Klis 
1996a, 1996b, 1996d, 1997b), GX 17+2 (Wijnands et 
al. 1997c), and Cyg X-2 (Wijnands et al. 1998). In 
addition to two simultaneous kilohertz QPOs, these 
power spectra show horizontal branch oscillations, 
which appear to be magnetospheric beat-frequency 
oscillations (see § 2 and Psaltis et al. 1998). 

These power spectra indicate that in the Z sources, 
some of the gas in the disk couples strongly to the 
weak stellar magnetic field at two or three stellar radii 
and is funneled into hot spots that could modulate the 
mass flux from the sonic point at the sonic-point beat 
frequency by periodically irradiating the clumps at 
the sonic radius (see § 2). It is natural to expect that 
the magnetic fields of many of the atoll sources are 
only a little weaker than those in the Z sources and 
are therefore still strong enough to produce a QPO 
at the sonic-point beat frequency. This expectation 
is supported by the spectral modeling described in 
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§ 2.2. 

Accretion via a geometrically-thin Keplerian disk 
flow, which is required to explain the high coher- 
ence of the kilohertz QPOs, is expected in the atoll 
sources, because at their low accretion rates the inner 
part of the accretion disk is likely to be gas-pressure- 
dominated. In the Z sources, which are thought to 
be accreting at close to the Eddington critical rate, 
gas in the Keplerian disk flow that is not channeled 
by the stellar magnetic field is likely to be pinched by 
the stellar field into a geometrically thin disk. 

The final requirement listed above for the sonic- 
point model, namely that the optical depth of the 
hot central corona exceed ^ 3, follows directly from 
the spectral modeling discussed in § 2.2. 

In summary, the physical picture of the atoll and 
Z sources that was developed prior to the discovery 
of the kilohertz QPOs, based on their 2-20 keV X- 
ray spectra and 1-100 Hz X-ray variability, give the 
magnetic field strengths, accretion fiows, and electron 
scattering optical depths required by the sonic-point 
model. 

4.2. Kilohertz QPO and Neutron Star Spin 
Frequencies 

In the sonic-point model, the frequency of the 
higher-frequency kilohertz QPO (when two are present) 
is the Keplerian orbital frequency uks at the point 
in the disk flow where the inward radial velocity in- 
creases rapidly within a small radial distance, whereas 
the frequency of the lower-frequency kilohertz QPO 
is approximately the difference between z/ks and the 
stellar spin frequency t^spin- In this section we de- 
scribe the range of sonic-point Keplerian and neutron 
star spin frequencies and the variation of the sonic- 
point Keplerian and beat frequencies with accretion 
rate expected in the model and then compare these 
expectations with the observations. We discuss the 
expected amplitudes of oscillations at other frequen- 
cies in § 4.3. 

Similarity of the sonic-point Keplerian QPO fre- 
quencies in different sources. A key question is why 
the higher- frequency kilohertz QPOs in different sources 
always have frequencies between ~ 500 Hz and 
^ 1200 Hz, even though the neutron stars in these 
sources have time-averaged mass accretion rates that 
differ by more than a factor of 100 and magnetic fields 
that are thought to differ by more than a factor of 10. 
In the sonic-point model there are three main reasons 



why the Kcplcrian-frequency QPO is restricted to this 
frequency range: 

(1) As discussed in § 3.2, there is an upper bound 
on the radius Rs of the sonic point, which is set by 
the maximum fraction of the angular momentum of 
the accreting gas that can be removed by radiation 
from the stellar surface and is ~ 3i?nisj where i?nis 
is the radius of the marginally stable orbit. There 
is also a lower bound on R^, which is Ems ff the ra- 
dius R of the neutron star is smaller than R^s, or 
R otherwise. Hence, the Keplerian frequency at the 
sonic point is confined to a similar interval in the atoll 
and Z sources, regardless of their accretion rates and 
magnetic field strengths. As a result, if the neutron 
stars in the kilohertz QPO sources have typically ac- 
creted a few tenths of a solar mass and therefore 
have masses 1.7 Mq, these constraints on the sonic- 
point Keplerian frequency mean that the frequencies 
of their sonic-point Keplerian QPOs will fall between 
~400 Hz and - 1300 Hz. 

(2) Modeling of the 2-20 keV X-ray spectra and 
1-100 Hz power spectra of the atoll and Z sources in- 
dicates that the magnetic fields of these neutron stars 
are positively correlated with their time-averaged ac- 
cretion rates, i.e., sources with higher accretion rates 
appear to have stronger magnetic fields (see § 2.2 and 
Psaltis & Lamb 1998a, 1998b, 1998c). This has an 
important implication for the frequency range of the 
kilohertz QPOs, because the stronger the stellar mag- 
netic field, the larger the fraction of the gas in the 
disk that is likely to couple to the magnetic field at 
two or three stellar radii and be channeled out of the 
disk flow there. As a result, the stronger the mag- 
netic field, the smaller the mass flux Mgp through the 
disk at the sonic point. Hence, Ms,, differs by a much 
smaller factor in the Z and atoll sources than does the 
total mass flux M onto the neutron star. The sonic- 
point Keplerian frequency is governed more by 
Mgp than by M, so the positive correlation between 
magnetic fleld and accretion rate means that the fre- 
quency ranges of the sonic-point QPOs are likely to 
be more similar in a collection of neutron-star LMXBs 
than one would expect, based only on the different 
ranges of M in different systems. 

(3) In addition, the inner disk is expected to be 
thicker if Mgp is large. As a result, the radial optical 
depth through the disk at a given accretion rate is 
smaller, which allows radiation from the neutron star 
surface to penetrate further out into the disk than 
one would expect based only on the different values 
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of M in the different systems. This tendency also acts 
to make the range of sonic-point radii, and hence the 
frequency ranges of the Keplerian-frequency QPOs, 
more similar in different systems than consideration 
of the accretion rates alone would suggest (see § 3.2). 

Expected and inferred neutron star spin rates. — In 
the sonic-point model, the separation between the fre- 
quencies of the QPOs in a pair is approximately equal 
to the stellar spin frequency. The frequency separa- 
tions observed in the kilohertz QPO sources so far (see 
Table ^) all indicate spin rates of a few hundred Hertz, 
consistent with the spin rates expected in the magne- 
tospheric beat-frequency model of the HBO (see § 2, 
Lamb et al. 1985, and Ghosh & Lamb 1992). 

Oscillations have been observed during thermonu- 
clear X-ray bursts from five kilohertz QPO sources 
so far (again see Table Only a single oscillation 
has been observed from each source during a given 
burst, the oscillations in the tails of bursts appear 
to be highly coherent (see, e.g., Smith et al. 1997), 
and the frequencies are always the same in a given 
source. Furthermore, comparison of burst oscillations 
from 4U 1728—34 over about a year shows that the 
timescale for any variation in the oscillation frequency 
is >3,000 yr (see, e.g., Strohmayer 1997). The burst 
oscillations are thought to be caused by emission from 
one or two regions of brighter X-ray emission on the 
stellar surface (see Strohmayer et al. 1997b), produc- 
ing oscillations at the stellar spin frequency or its first 
overtone, respectively. 

The 363 Hz frequency of the burst oscillation ob- 
served in 4U 1728—34 is consistent with the sepa- 
ration frequency of its two simultaneous kilohertz 
QPOs. The 524 Hz and 581 Hz frequencies of the 
burst oscillations seen in KS 1731—260 and 
4U 1636—536 are probably twice the spin frequencies 
of these neutron stars (see Morgan & Smith 1996; 
Smith et al. 1997; Wijnands & van der Klis 1997; 
Zhang et al. 1996, 1997). If so, the spin rates of 
these stars are 262 Hz and 290 Hz, respectively. The 
frequencies of the fairly coherent oscillations seen at 
549 Hz in Aql X-1 (Zhang et al. 1998) and at 589 Hz in 
the so-far-unidentified source near the Galactic center 
(Strohmayer et al. 1996d) may also be twice the spin 
frequencies of these neutron stars. These spin rates 
are all consistent with those expected in the sonic- 
point model. 

Similarity of the neutron-star spin frequencies in 
different sources. — If the neutron stars in the kilo- 
hertz QPO sources have spin frequencies comparable 



to their equilibrium spin frequencies, then we expect 
their spin frequencies to be a few hundred Hertz. This 
can be seen from equation (|^) in § 2.2, which can be 
solved for the equilibrium spin frequencies at which 
continued accretion at the given rate leaves the spin 
frequency unchanged, with the result 
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(25) 

Here Wc is the critical fastness (Ghosh & Lamb 1979b). 
Hence, for tJc ~ 1, a 1.7 Af© atoll source with /io.27 = 
0.5 accreting from a GPD disk at a rate A/^ = 
0.001 Me has an equilibrium spin frequency of 230 Hz, 
which is similar to the 240 Hz equilibrium spin fre- 
quency of a 1.4 AfQ Z source with ^0,27 = 2 accreting 
from an RPD disk at a rate Mi = 0.5 Me. 

A recent analysis (Psaltis et al. 1998) of the prop- 
erties of the HBOs and kilohertz QPOs in a collection 
of five Z sources shows that they are consistent with 
many of the predictions of the magnetospheric beat- 
frequency model of the HBO and that the agreement 
is best if, as expected in this model, they are all in 
spin equilibrium (see also Ghosh & Lamb 1992). In 
particular, this analysis shows that the narrow range 
of spin frequencies observed in the Z sources is to be 
expected if they are in spin equilibrium. 

Variation of kilohertz QPO frequencies with mass 
accretion rate. — In the sonic-point model, the fre- 
quencies of the kilohertz QPOs are expected to rise 
as the accretion rate increases, within the bounds set 
by the orbital frequencies at the angular momentum 
loss radius and at the radius of the marginally sta- 
ble orbit (see above and § 3.2). In order to see why 
an increase is expected, suppose that the sonic radius 
is at Rsi for a given accretion rate. An increase in 
the accretion rate will cause the radial optical depth 
from the stellar surface to Rsi to increase, if all other 
physical quantities remain constant, and hence the 
sonic point will move inward to the smaller radius 
Rs2 at which the optical depth from the stellar sur- 
face is approximately the same as before, causing an 
increase in the orbital frequency at the sonic radius. 
The X-ray luminosity between bursts is produced al- 
most entirely by accretion and hence an increase in 
the accretion rate causes an increase in the luminos- 
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ity. We therefore expect a strong, positive correlation 
between the frequencies of the sonic-point Keplerian- 
and beat-frequency QPOs and the persistent X-ray 
luminosity. 

These basic ideas are illustrated by the general rel- 
ativistic accretion flow calculations of § 3.2. As shown 
in Figure ^, the sonic-point Kcplcrian and beat fre- 
quencies given by these calculations increase steeply 
with increasing mass accretion rate until the sonic 
point reaches Rms, at which point their frequencies 
stop changing. This flattening of the two frequency 
vs. accretion luminosity relations is one possible sig- 
nature of the existence of a marginally stable or- 
bit (see also § 5.4). Coherent Keplerian- and beat- 
frequency QPOs are unlikely to be produced if the 
sonic point moves inward, close to the stellar surface, 
because of the disruptive effect of the surface mag- 
netic field and the strong viscous shear layer that is 
expected to develop if the Keplerian flow interacts di- 
rectly with the neutron star surface. 

We emphasize that these calculations are intended 
to be illustrative rather than to reproduce the QPO 
frequency behavior of any particular source. If the 
mass of the neutron star is 1.7 Mq, rather than 
1.4 M0 as assumed in these calculations, the limit- 
ing frequency would be 1.3 kHz rather than 1.6 kHz 
(see Miller, Lamb, & Psaltis 1998). Again, if the ac- 
cretion disk is thicker near the star than is assumed in 
these calculations, more of the inner disk will be illu- 
minated and hence the angular momentum loss radius 
-Rami will be larger, causing the sonic-point Keplerian 
frequency at a given accretion rate to be lower. If 
the disk thickness changes with accretion rate, then 
the variation of the sonic-point Keplerian frequency 
with accretion rate will be different from the variation 
found in the calculations reported here, which treat 
the disk thickness as constant. Despite these uncer- 
tainties, the calculations described here show that an 
increase in the mass accretion rate leads naturally to 
an increase in the sonic-point Keplerian frequency. 

We stress that caution must be used in comparing 
the predicted variation of QPO frequencies with ac- 
cretion with the countrate measured by a given X-ray 
detector, because countrates are known to be poor in- 
dicators of the mass accretion rate and the accretion 
luminosity, at least for some sources at some times. 
The most likely reason is that a particular instru- 
ment measures the photon number flux over only a 
restricted energy range and that a change in the mass 
accretion rate typically causes a change in the shape 



of the X-ray spectrum as well as its normalization. 
Hence the detector countrate typically is not propor- 
tional to the accretion rate. 

In the Z sources, the mass accretion rate is known 
to be different when the EXOSAT countrate is the 
same but the source is on a different branch of the 
Z track (Hasinger & van der Klis 1989; Lamb 1989; 
Psaltis et al. 1995). Moreover, in the "Cyg-like" Z 
sources, the normal/flaring branch vertex, which is 
thought to correspond to a mass accretion rate ap- 
proximately equal to Me (Lamb 1989), occurs at dif- 
ferent countrates during different observations (Ku- 
ulkers 1995). 

The observations of Ford et al. (1997a) show that 
the properties of the accretion flow are not uniquely 
related to the X-ray countrate in the atoll sources. 
In these observations, the dependence of the kilo- 
hertz QPO frequencies on countrate seen in data on 
4U 0614-F091 taken with RXTE in 1996 August dif- 
fers, both in slope and in normalization, from the de- 
pendence of the QPO frequencies on countrate seen 
in data taken in 1996 April. In contrast, when the 
4U 06144-091 QPO frequencies are plotted against 
the energy of the spectral peak, which may be a 
good indicator of the accretion rate (Psaltis & Lamb 
1998c), the slope and normalization of the curves are 
the same for both observations (Ford et al. 1997b). 

Clearly, bolometric and other corrections are typ- 
ically important, so comparison of QPO frequency 
versus countrate data with the QPO frequency ver- 
sus accretion rate predictions of theoretical mod- 
els must be approached with caution. Despite this 
caveat, it is worth emphasizing that there is a strong, 
positive correlation between the frequencies of the 
kilohertz QPOs seen in 4U 1728-34, 4U 0614-^091, 
4U 1820-30, and KS 1731-260, and the countrate 
measured by RXTE, as expected in the sonic-point 
model if the countrate increases with increasing ac- 
cretion rate. 

It is interesting to speculate about the dependence 
of kilohertz QPO frequencies on luminosity that is to 
be expected if such QPOs are detected during a ther- 
monuclear X-ray burst. In this case the mass accre- 
tion rate and the radial optical depth are no longer 
tied to the luminosity, and hence the expected de- 
pendence of j/Ks on luminosity is different from what 
is expected during the periods between type I X-ray 
bursts. When radiation forces are extremely strong, 
such as at the peak of a type I X-ray burst that causes 
photospheric radius expansion, the radiation controls 
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the accretion flow out to large radii and we there- 
fore do not expect any QPO at the sonic-point Kep- 
lerian frequency or at the sonic-point beat frequency. 
However, QPOs at these frequencies may reappear in 
the tail of the burst. While the burst luminosity is 
several times the persistent accretion luminosity, we 
expect the sonic point to be farther from the star 
than when the luminosity is lower. Thus, we expect 
that if sonic-point QPOs are detected in the tail of 
a thermonuclear burst, their frequencies will increase 
as the luminosity declines. Note, however, that dur- 
ing the decay phase of such a burst, the luminosity 
changes rapidly. If the frequencies of the QPOs track 
this change, the QPO peaks in power spectra will be 
smeared out, unless the spectra are constructed from 
short ( <0.1 s) segments of data. 

So far, no kilohertz QPOs have been observed near 
the maximum of type-I X-ray bursts, which is consis- 
tent with what is expected in the sonic-point model. 
More sensitive searches will be required to determine 
if kilohertz QPOs are present during the decay phase 
of bursts. 

4.3. Kilohertz QPO Amplitudes and Coher- 
ence 

Amplitudes of oscillations at different frequencies. — 
A key question that any theory of the kilohertz QPOs 
must answer is why at most two strong QPOs have 
so far been seen simultaneously in the kilohertz QPO 
sources. In addressing this question, it is important to 
consider both the generation of oscillations near the 
neutron star and the effects of propagation of the ra- 
diation through the gas surrounding the neutron star. 
In order to be detectable far away, an oscillation must 
be strong at the source or have a small attenuation, 
or both. 

The QPOs at the sonic-point Keplerian and beat 
frequencies are the only QPOs expected to have mod- 
erately high amplitudes far from the neutron star be- 
cause the optical depth and luminosity oscillations 
that produce them generate relatively high ampli- 
tudes near the star, because these frequencies (and 
their overtones) are the only frequencies generated 
by the sonic-point mechanism in lowest order, and be- 
cause scattering of photons by the gas surrounding the 
neutron star selectively suppresses the already weak 
higher-frequency oscillations generated at higher or- 
ders. 

As discussed in § 3.3, to lowest (first) order, the 



sonic-point mechanism generates oscillations only at 
the sonic-point Keplerian and beat frequencies (see 
Table |^). As explained there, only very weak over- 
tones of the sonic-point Keplerian frequency are likely 
to be generated, because the angular distribution of 
the radiation from the bright impact footprints of the 
rotating density pattern produced by clumps is ex- 
pected to be very broad. 

The radiation pattern that rotates with the star 
is also expected to be broad, but at the sonic ra- 
dius, where it interacts with the orbiting gas to gener- 
ate the luminosity oscillation at the sonic-point beat 
frequency, it is not likely to be perfectly sinusoidal. 
Hence overtones of the sonic-point beat frequency vbs 
may be generated in the inward mass flux from the 
clumps at the sonic radius, producing overtones of z^bs 
in the X-ray flux from the star. As radiation from the 
stellar surface propagates outward through the part of 
the central corona that extends beyond the sonic ra- 
dius, the radiation pattern that rotates with the star 
is likely to be broadened further (see below) , reducing 
the amplitudes of oscillations at the stellar spin fre- 
quency and its overtones seen by a distant observer. 

In addition to the QPOs at the sonic-point Keple- 
rian and beat frequencies and their overtones, in sec- 
ond order the sonic-point mechanism also generates 
QPOs at 

t'spin and 2(^X3 ~ ^'spin and in third order at 

i^Ks — 2l'spin, i^Ks + i^spin, 3j^Ks — 2t'spin, and 3z^Ks " t'spin 

(again see Table |^). These high-order oscillations are 
expected to be very weak, even near the neutron star. 

The attenuation of an oscillation depends strongly 
on whether it is a beaming oscillation or a luminos- 
ity oscillation (see § 3.6). A beaming oscillation is an 
oscillation produced by rotation of an angular radia- 
tion pattern, like the beam of a lighthouse; in a pure 
beaming oscillation, the luminosity remains constant 
as the radiation pattern rotates. In contrast, a lumi- 
nosity oscillation is one produced by an oscillation of 
the luminosity of the source; in a pure luminosity os- 
cillation, the radiation pattern remains static as the 
luminosity oscillates. In the sonic-point model, the 
sonic-point Keplerian frequency is the only intrinsi- 
cally strong beaming oscillation, and the sonic-point 
beat frequency is the only luminosity oscillation pro- 
duced by the accretion flow. 

In passing through a scattering region, luminos- 
ity oscillations are attenuated only by time-of-flight 
smearing whereas beaming oscillations are also at- 
tenuated by the isotropization of the radiation pat- 
tern by the scattering. Hence beaming oscillations 
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are weakened more by propagation through a scatter- 
ing corona. Figure |o| compares the attenuation of 
luminosity and beaming oscillations produced at the 
center of a uniform, spherical, scattering cloud of ra- 
dius Rc = 3 X 10^ cm, as a function of the optical 
depth of the cloud. The size of this cloud is compara- 
ble to the dimensions of the small scattering coronae 
with optical depths ~ 3-5 that are thought to sur- 
round the neutron stars in the atoll and Z sources 
(Lamb 1989, 1991; see also § 2.2). The frequencies 
of the oscillations shown in Figure |l^ have been cho- 
sen to represent a hypothetical accreting neutron star 
with j^Ks — 1100 Hz and i^spin = 300 Hz. The beat- 
frequency luminosity oscillation is therefore at 800 Hz. 
The attenuation factors have been calculated using 
equations (|2^) and (p^). Figure |l^ shows that the 
beaming oscillation at the sonic-point Keplerian fre- 
quency must have a relatively high intrinsic amplitude 
in order to produce a strong, observable QPO for a 
distant observer, whereas the luminosity oscillation at 
the sonic-point beat frequency needs only to have a 
moderate amplitude. 

The amplitudes of the oscillations at the overtones 
of either luminosity or beaming oscillations seen by a 
distant observer are expected to be far smaller than 
the amplitudes seen at the fundamental frequencies, 
for several reasons. For example, beaming patterns 
with multiple lobes arc much more easily isotropized 
by scattering than are beaming patterns with single 
lobes (see, e.g., Brainerd & Lamb [1987]). More- 
over, for either beaming or luminosity oscillations, 
overtones of an oscillation are at higher frequencies 
than the fundamental and are thus more susceptible 
to time-of-flight smearing. 

As a specific example, we consider 4U 1728—34. 
As interpreted within the sonic-point model, the spin 
frequency of this neutron star is 363 Hz (see Table 3). 
During one observation, vks was 1045 Hz. Assum- 
ing that this neutron star is surrounded by a central 
corona with a radius Rc = 3x10^ cm and an optical 
depth r = 5, the amplitude of the luminosity oscilla- 
tion at the beat frequency vks — ^'spin = 682 Hz seen 
by a distant observer is about 85% of the amplitude 
at the neutron star. For comparison, the beaming os- 
cillation at j^Ks + t'spin — 1408 Hz , which is generated 
only in third order and is therefore likely to be very 
weak even near the neutron star, has an amplitude at 
infinity which is only about 20% of its amplitude at 
the star. The amplitudes of oscillations at the over- 
tones of the sonic-point Keplerian frequency are all 



reduced by factors ^J20. For these reasons, it is not 
surprising that the only strong QPOs seen are at the 
sonic-point Keplerian and beat frequencies. 

Although oscillations at frequencies other than the 
sonic-point Keplerian and beat frequencies are likely 
to be very weak, their amplitudes are surely not zero. 
We therefore expect that, if neutron-star LMXBs 
are observed with sufficient sensitivity, oscillations at 
other frequencies will be detected. In particular, a 
careful search near the spin frequency, at the first 
overtone of the beat frequency, and near the sum of 
the Keplerian and spin frequencies in sources with 
pairs of high-frequency QPO peaks may reveal very 
weak QPOs at these frequencies. 

The oscillation at the beat frequency is a lumi- 
nosity oscillation, but in the sonic-point model it is 
created by interaction of the weakly beamed radiation 
pattern that is rotating at the stellar spin frequency 
with the clumps of gas that are orbiting at or near 
the sonic point. It is therefore important to empha- 
size why the beaming oscillation at the stellar spin 
frequency is strong enough to modulate significantly 
the mass inflow rate from clumps at the sonic point, 
yet too weak to produce a significant peak in power 
spectra constructed by a distant observer. 

Previous modeling of the X-ray spectra (Lamb 
1989, 1991; Psaltis et al. 1995; Psaltis & Lamb 1998a, 
1998b, 1998c) and the 1-100 Hz X-ray variability of 
the atoll and Z sources (Lamb 1989, 1991; Miller & 
Lamb 1992) strongly indicates that a hot, central, 
Comptonizing corona surrounds the neutron star, ex- 
tending several stellar radii from its surface (see § 2.2). 
The sonic radius is well inside this corona, so the scat- 
tering optical depth from the stellar surface to the 
sonic point is much less than the scattering optical 
depth from the stellar surface to infinity. The finite 
size of the neutron star also diminishes the effect of 
attenuation on the amplitude of X-ray brightness os- 
cillations for gas orbiting close to the star. These 
factors make the brightness oscillation produced by 
the radiation pattern that rotates with the star much 
stronger at the sonic radius than it is far away. 

Dependence of QPO amplitudes on accretion rate. — 
In the sonic-point model, the ratio of the amplitude 
of the QPO at the sonic-point Keplerian frequency 
to the amplitude of the QPO at the sonic-point beat 
frequency generally depends on the strength of the 
neutron star magnetic field. If the magnetic field is 
too weak to channel gas even near the stellar surface, 
the amplitude of the oscillation at the beat frequency 
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will be too small to be detected (see § 3.1). However, 
the ratio of the amplitudes of the QPOs at different 
frequencies is also expected to depend on the mass 
accretion rate. 

For example, as the mass accretion rate rises, the 
electron density and optical depth of the gas around 
the neutron star also rise. Hence, the amplitude at 
infinity of the beaming oscillation at the sonic-point 
Keplerian frequency will fall faster relative to the 
amplitude of the luminosity oscillation at the sonic- 
point beat frequency, all else being equal; indeed, 
the QPO at the Keplerian frequency may become 
undetectable while the QPO at the beat frequency 
remains strong. This expectation is consistent with 
the observed behavior of the kilohertz QPO pair in 
4U 1728—34 (Strohmayer et al. 1996c), assuming that 
the increasing mass accretion rate causes an increase 
in the countrate (but see § 4.2). In this source, the 
amplitude of the higher-frequency QPO decreases rel- 
ative to the amplitude of the lower-frequency QPO as 
the countrate increases. 

Alternatively, if the sonic point moves far enough 
away from the stellar surface, as may happen at low 
accretion rates (and hence countrates), the optical 
depth through the flow from the stellar surface to 
the sonic radius may become large enough to reduce 
greatly the amplitude at the sonic point of the ra- 
diation pattern that rotates with the star, thereby 
suppressing the luminosity oscillation at i>bs- This 
is consistent with an observation of 4U 1728—34 in 
which the lower-frequency QPO is not observed when 
the frequency of the higher-frequency QPO is low 
(Strohmayer et al. 1996c). 

This analysis shows that if only a single high- 
frequency QPO peak is observed, it could be ei- 
ther the sonic-point Keplerian-frequency QPO or the 
sonic-point beat-frequency QPO, depending on the 
magnetic field of the neutron star and the mass accre- 
tion rate at the time. In this case, a secure identifica- 
tion can only be made by considering other properties 
of the QPO and comparing them with the properties 
of the QPOs seen in other observations of the same 
source. 

Dependence of QPO amplitudes on photon energy. — 
In the sonic-point model, the relative amplitudes of 
the Keplerian frequency and beat frequency QPOs 
are expected to increase steeply with photon en- 
ergy over the ~ 5-10 keV energy range (see § 3.5). 
This is because the optical depth is expected to os- 
cillate at both frequencies: the optical depth along 



the line of sight from the stellar surface oscillates 
at the Keplerian frequency, while the total optical 
depth of the scattering region oscillates at the beat 
frequency as the density of the accreting gas falling 
on the stellar surface oscillates at this frequency. As 
shown in § 3.5, a modest oscillation in the optical 
depth produces a QPO with a large relative ampli- 
tude at high photon energies. This is consistent with 
the steep increase of QPO amplitude with increas- 
ing photon energy observed in the higher-frequency of 
the two simultaneous QPOs in 4U 1636—536 (Zhang 
et al. 1996), 4U 0614-f091 (Ford et al. 1997b), and 
KS 1731-260 (Wijnands & van der Khs 1997), in 
the lower-frequency of the two simultaneous QPOs in 
4U 1608-52 (Berger et al. 1996; Mendez et al. 1998) 
and 4U 1728-34 (Strohmayer et al. 1996c), and in 
the kilohertz QPOs in the Z sources GX 5—1 (van der 
Klis et al. 1996e), GX 17+2 (Wijnands et al. 1997c), 
and Cyg X-2 (Wijnands et al. 1998). 

In § 3.5 we used a simple analytical model to de- 
rive an expression for the photon energy dependence 
of the rms amplitude of a QPO produced by oscilla- 
tions in the optical depth and in the injection rate 
of soft photons. Figure [ill compares the results of a 
more detailed numerical calculation performed using 
the algorithm of Miller & Lamb (1992) with ampli- 
tude data for the lower-frequency kilohertz QPO ob- 
served in 4U 1608—52 and the higher-frequency kilo- 
hertz QPO seen in 4U 1636—536. In performing this 
calculation we made the same assumptions as in § 3.5 
and assumed further that the injection rate of soft 
photons is constant in time, that the spectrum of 
the injected photons is a blackbody at temperature 
kT — 0.6 keV, that the electron temperature in the 
central corona is kT^ = 9 keV, and that the optical 
depth varies from r = 3 to r = 3.3 during an oscilla- 
tion. The model matches the data well. The geome- 
try of the actual upscattering region is undoubtedly 
much more complicated than assumed in this calcula- 
tion, but the excellent correspondence with the data 
and the ubiquity of the steep increase of amplitude 
with photon energy over the 5-10 keV energy range 
suggest that the model has many elements in common 
with the true physical situation. 

As a source moves in the X-ray color-color dia- 
gram and its spectrum changes (implying a change 
in the average optical depth and electron tempera- 
ture of the Comptonizing region) , we expect that the 
photon energy dependence of the rms amplitude of 
the oscillations will also change. Moreover, as dis- 
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cussed in §3.5, we expect that, in general, the depen- 
dence of the beat-frequency QPO ampUtude on pho- 
ton energy will be different from the dependence of 
the Keplerian-frequency QPO amplitude. This is con- 
sistent with the observations of 4U 0614-1-091 (Ford et 
al. 1997b) in which amplitude versus photon energy 
curve of the beat-frequency QPO has a minimum at 
^ 10 keV, whereas the amplitude versus photon en- 
ergy curve of the Keplerian-frequency QPO increases 
monotonically from 2 keV to 20 keV. Because the 
QPOs amplitude depend strongly on photon energy, 
to be meaningful a comparison of QPO amplitudes 
measured for different sources or at different times 
for the same source must consider the same range of 
photon energies. 

Coherence of the kilohertz QPOs. — In § 3.4 we 
showed that the high coherence of the Kcplerian and 
beat-frequency QPOs in the sonic-point model is pri- 
marily a consequence of the extremely sharp increase 
in the inward radial velocity near the sonic point, 
which maps a small range of orbital frequencies onto 
the stellar surface. There we also considered the deco- 
herence caused by destruction of clumps by turbulent 
dissipation within the disk flow, by advection to the 
stellar surface, and by decay as gas is stripped from 
them by the supersonic flow at the sonic point. We 
concluded that the QPO peaks can be as narrow as 
observed [yj^v ^ 30-200) for reasonable conditions 
in the accretion flow. 

The FWHM Jz^bf of the sonic-point beat fre- 
quency oscillation is expected to be comparable to 
the FWHM Ji'kf of the sonic-point Keplerian fre- 
quency oscillation, because the beat-frequency QPO 
is produced by the beat of a nearly periodic signal at 
frequency I'spin against the Keplerian frequency (see 
§ 3.4). We do not expect the FWHM of the two kilo- 
hertz QPO peaks to be identical, however, because 
there are processes that can affect the FWHM of one 
but not the other. For example, the FWHM of the 
beat-frequency peak depends in part on the range of 
radii over which radiation forces can affect the mass 
accretion rate significantly. This range can be either 
less than or greater than the range of radii over which 
the inward radial velocity increases rapidly, and hence 
iJi^BF can be either less than or greater than (5i^kf- 
This effect can also displace the centroid of the beat- 
frequency peak relative to the centroid of the Keple- 
rian frequency peak, so that the observed frequency 
difference t'Ks — ^Bs is close to, but not exactly equal 
to, the stellar spin frequency. 



4.4. Amplitudes of QPOs Produced by Stars 
with Different Magnetic Field Strengths 

In the sonic-point model, the larger the magne- 
tosphere, the smaller the fraction of the accreting 
gas that reaches the surface of the star without cou- 
pling to the magnetic field, and hence the smaller the 
amplitude of the kilohertz QPOs. Therefore, we ex- 
pect the rms amplitudes of the kilohertz QPOs to be 
roughly anticorrelated with the strength of the stellar 
magnetic field, if all other physical quantities remain 
fixed. As discussed in § 2.2, the "4U" atoll sources 
are thought to have the weakest magnetic fields, the 
"GX" atoll sources and the "Sco-like" Z sources are 
thought to have somewhat stronger fields, and the 
"Cyg-like" Z sources are thought to have the strongest 
fields (Psaltis & Lamb 1998a, 1998b, 1998c). Hence, 
we expect kilohertz QPOs to be common and strong 
in the "4U" atoll sources, but weak or even unde- 
tectable with current instruments in the "GX" atoll 
sources and the Z sourcesF] This is indeed the case, 
as is evident from Table |lf If, as indicated by the 
spectral modeling discussed in § 2, the magnetic field 
of Cir X-1 is weak, we expect it to exhibit sonic-point 
QPOs when it is in its low state. We also expect that 
kilohertz QPOs will be undetectable by current in- 
struments in any sources that have strong magnetic 
fields and therefore produce strong, periodic oscilla- 
tions at their spin frequencies. 

The relation between kilohertz QPO amplitude 
and magnetic field strength may be made semi-quanti- 
tative using the spectral calculations (see § 2.2) of 
Psaltis et al. (1995) and Psaltis & Lamb (1998a, 
1998b, 1998c). These calculations indicate that LMXBs 
containing neutron stars with weaker magnetic fields 
generally have larger hard X-ray colors (e.g., the ratio 
of the 7-20 keV countrate to the 5-7 keV countrate). 
Based on this physical picture, we expect the rms am- 
plitudes of the kilohertz QPOs produced by the sonic- 
point mechanism to be higher in sources with larger 



hard colors. This trend is evident in Figure 12, which 



shows the rms amplitudes versus the hard X-ray col- 
ors of four sources observed with the PGA detector 
onboard the RXTE satellite. If other sources with 
kilohertz QPOs follow this same trend, this will be 
strong support for the sonic-point model. 



^ After this specific prediction was made in the originally sub- 
mitted and circulated version of the present paper, kilohertz 
QPOs with very low amplitudes were detected in all six of the 
originally identified Z sources; see Table hi. To date, kilohertz 
QPOs have not been detected in any of the GX atoll sources. 
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As discussed in § 2.3, the observed anticorrelation 
between kilohertz QPO amplitude and magnetic field 
strength is also strong evidence that the magneto- 
spheric beat-frequency mechanism not the correct ex- 
planation of these oscillations. 

5. IMPLICATIONS FOR NEUTRON STARS 
AND DENSE MATTER 

In the sonic-point model, the higher-frequency QPO 
in a kilohertz QPO pair has a frequency equal to the 
orbital frequency of gas near the sonic point; the rel- 
atively high coherence of this QPO is a consequence 
of the fact that the gas clumps that generate it are 
in nearly circular orbits. In this section we show that 
the inferred existence of a nearly circular orbit around 
a neutron star with a frequency in the kilohertz range 
can be used to derive interesting new upper bounds on 
the mass and radius of the star and constraints on the 
equation of state of the dense matter in all neutron 
stars. For simplicity we discuss first the case of a non- 
rotating star around which the radial component of 
the radiation force is negligible. We then consider the 
changes in the mass and radius constraints caused by 
frame-dragging when the star is spinning and by the 
radial component of the radiation force. We show that 
rotation at ~ 300 Hz typically increases the bound on 
the mass by ~ 20% and the bound on the radius by a 
few percent. In contrast, the radiation force reduces 
the upper bounds on the masses and radii by a few 
percent or less in the atoll sources but perhaps by 
much larger percentages in the Z sources. 

Observations of kilohertz QPOs may be able to 
establish the existence of an innermost stable circular 
orbit around some neutron stars (see § 4.2). If this 
can be accomplished, it would be the first evidence 
concerning a prediction of general relativity in the 
strong-field regime. If furthermore the frequency of a 
particular kilohertz QPO can be securely established 
as the orbital frequency at the radius of the innermost 
stable circular orbit and if the spin frequency of the 
neutron star can be determined, then the frequency of 
the QPO can be used to fix the mass of the neutron 
star for each assumed equation of state, tightening 
the constraints on the properties of dense matter and 
possibly ruling out many currently viable equations 
of state. We discuss how this can be done and the 
evidence that would signal detection of a QPO with 
the orbital frequency of the marginally stable orbit. 



5.1. Nonrotating Star 

Suppose that the frequency of the higher-frequency 
QPO in a kilohertz pair is j^qpo2 and that, as in the 
sonic-point model, i'qpo2 is the orbital frequency of 
gas in a nearly circular Keplerian orbit around the 
neutron star. Assume for now that the star is not ro- 
tating and is spherically symmetric. Then the exte- 
rior spacetime is the Schwarzschild spacetime. In this 
spacetime, the orbital frequency (measured at infin- 
ity) of gas in a circular orbit at Boyer-Lindquist radius 
r around a star of mass M is (see eq. [|oj) 

z.^(M,r) = (l/27r)(GAf/r3)i/2 . (26) 

Here and below the superscript zero indicates that 
the relation is that for a nonrotating (j — 0) star. If 
the mass of the star is known, equation ( p6| ) may be 
solved for the orbital radius i?oib where the Keplerian 
frequency is i^QP02, with the result 

RU{M, z^qpo2) = (GAf/4^2^2p^^)i/3 (27) 

Conversely, if the orbital radius of the gas is known, 
equation (^6|) may be solved for the mass of the star 
that gives an orbital frequency equal to j^qpo2, with 
the result 

M°(i?orb, i^QP02) = (4^VG)i?;^,bi'QP02 • (28) 

If, as is so far the case for the kilohertz QPO sources, 
neither M nor i?orb are known, equations (|2^) and ( p8|) 
do not determine i?orb or M but do establish a rela- 
tion between them. In the radius-mass plane, this 
relation is a curve that begins at the origin and rises 
up and to the right. 

As a specific example, the dashed curve marked 
M'^iRorh) in Figure [l^ shows the relation given by 
equation ( p8|) for i'qpo2 — 1220 Hz, a value of i^qpo2 
observed in 4U 1636—536 (W. Zhang, personal com- 
munication). Hence, if 4U 1636—536 were not rotat- 
ing, the mass of the neutron star in this source and the 
orbital radius of the gas clumps producing the QPO 
during this observation would have to correspond to 
some point along this dashed curve. 

Consider now the constraints on the neutron star 
mass and radius that follow from the frequency of the 
higher-frequency QPO. Obviously, in order for gas to 
be in orbit, it must be outside the star. This means 
that for an orbit of given radius i?orb, the mass M of 
the star must be greater than A/°(i?orb, f^QP02); alter- 
natively, for a star of given mass M, the radius R of 
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the star must be less than R'^^^{M, i'qpo2)- Thus, the 
point that represents the mass and radius of the neu- 
tron star must he above the curve Af°(i?orb, t^QP02) 
in the radius-mass plane. The larger the value of 
fQP02, the higher the curve, so the most stringent — 
and hence the relevant — constraint on the mass and 
radius of the neutron star in a particular source is 
given by the highest value of i^qpo2 ever observed in 
that source, which we denote i^qpo2- For example, 
1220 Hz is the highest value of i'qpo2 seen so far in 
4U 1636—536, so this is the current value of i^Qpo2 
for this source (in fact, this is the highest value of 
t'QP02 seen so far in any source). Consequently, if 
the neutron star in 4U 1636—536 were not rotating, 
the point in the radius-mass plane that represents its 
mass and radius would have to lie above the curve 
M = M°(i?orb, 1220Hz), which is the dashed curve 
shown in Figure |l^. 

So far, the only information we have used to con- 
strain the mass and radius of the star is the fre- 
quency of the higher-frequency QPO in a kilohertz 
pair. However, we have the additional information 
that the coherence of the higher-frequency QPO is 
high (t'Qpo/'^'^QPO ^ 100) in many sources. Let us 
assume that the higher-frequency QPO in a pair 
is generated by the same mechanism in all sources 
that show such pairs. Then the high coherence of 
these higher-frequency QPOs imposes additional con- 
straints on the neutron star's mass and radius be- 
cause, in order to produce a QPO with such high co- 
herence, the gas that generates the QPO must be in 
a nearly circular orbit. Hence _Rorb(-^^, '^qpo2) must 
be greater than the radius Rms{M) of the innermost 
stable circular orbit, because gas inside R^s spirals 
quickly inward to the stellar surface. 

For a nonrotating star, the radius of the innermost 
stable circular orbit is a function only of the mass of 
the star and is given by R'^g{M) — 6GM/c^. Invert- 
ing this relation gives 



and 



M°iRn.s) - (cV6G)i?„ 



(29) 



In the radius-mass plane, relation (|29| ) is a straight 
line of slope -1-1 through the origin and is the dotted 
line marked M°(i?ins) in Figure |l^. This line inter- 
sects M*'(_Roib, i'qpo2) mass and radius values 



= 2.2(1000 Hz /i/Qpoa) Mg 



R 





max 



c{V2A7riy^P02) ^ 



= 19.5 (1000 Hz /j/^po2) km. (31) 

If M were larger than M^^^i the orbital radius 
-^orb(-^' ^QP02) of the gas generating the QPO would 
be less than R'^^{M), so M^^^ is an upper bound on 
the mass of the star; M^^^ is inversely proportional 
to i'QP02 but independent of the star's radius and 
the orbital radius of the gas clumps that are produc- 
ing the higher- frequency QPO. Similarly, i?^^^ is an 
upper bound on the radius of the star; is also in- 

versely proportional to i'qpo2 but independent of the 
star's mass and the orbital radius of the gas clumps 
that are producing the higher-frequency QPO. 

The heavy horizontal line plotted in Figure 13 
shows the upper bound on the mass of a nonro- 



tating star for 



^QP02 



= 1220 Hz. This upper 



(30) 



bound is 1.8 Mq, so the mass of the neutron star in 
4U 1636—536 would have to be less than this if it were 
not rotating. For i^qpo2 — 1220 Hz , i?J^ax is 16.0 km, 
so the radius of the neutron star in 4U 1636—536 
would have to be be less than this if the star were 
not rotating. 

Suppose now that the frequency i'qpo2 ^ P^'^" 
ticular QPO is securely identified as the orbital fre- 
quency of gas in the innermost stable circular orbit 
around a particular neutron star. Then i?oib = ^^ms 
for this QPO, so the representative point of the orbit 
is at the intersection of the diagonal Af'^(i?,„s) line 
and the Af''(_Rorb, i'Qpo2) curve. The mass of the star 
is therefore M^^^. Hence identification of a QPO 
frequency with the frequency of the innermost stable 
circular orbit immediately determines the mass of the 
star. The radius of the star is not determined by such 
an identification, but must still be less than -R^ax- 

As a specific example, suppose that the 1220 Hz 
QPO observed in 4U 1636—536 is securely identified 
as the orbital frequency of gas in the innermost sta- 
ble circular orbit around this neutron star. Then 
the mass of this neutron star would be determined 
as 1.8 Af©, if it were not rotating. 

These arguments apply to stars with arbitrary spin 
rates as well as to nonrotating stars, although the ex- 
pressions for i?orb(A/, i'qpo2) B.iid i?,„s(Af) are differ- 
ent for a rotating star and depend on the star's spin 
rate as well as its mass. They may be summarized as 
follows: 
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(1) Stellar radii R > Roih{M, i^qpo2) ^-^e excluded, 
because there is no Keplerian orbit with frequency 
J/QPQ2 outside a star with such a large radius; the 
value of RmsiM) is irrelevant. 

(2) If i? < i?orb(M,t^^P02) but RnUM) > 
Rovh{M,u^pQ2)j there is a Keplerian orbit with fre- 
quency i'qpo2 outside the star, but any oscillation 
produced by gas in this orbit would have a coherence 
much lower than that observed. 

(3) If i? < ii'orb(M,J^^P02) and i?,n.(M) < 
-RorbC^f, i^QP02)i there is a Keplerian orbit with fre- 
quency J^Qpo2 outside the star and the oscillation pro- 
duced by gas in this orbit can have the required high 
coherence. 

(4) If the frequency i^qpo2 of a particular QPO is 
securely identified as the orbital frequency of gas in 
the innermost stable circular orbit around a particular 
neutron star, the mass of the star is M^^^\ its radius 
must be less than i?max- 

For a nonrotating star with i'qpo2 = 1220 Hz , the 
combinations of stellar mass and radius allowed by 
condition (3) are the points in the radius-mass plane 
above the curve M = M"(i?orb, 1220 Hz ) but below 
M = M^g^j,. This region is outlined in Figure ^ by 
the heavy solid line. The point representing the mass 
and radius of the neutron star in 4U 1636—536 would 
have to lie in this region if the star were not rotat- 
ing. The allowed region collapses to the heavy solid 
horizontal line if the QPO frequency is identified as 
the orbital frequency of gas in the innermost stable 
circular orbit so that condition (4) applies. 

Figure |lj compares the mass-radius relations for 
nonrotating neutron stars given by five equations of 
state for neutron-star matter ranging from soft to 
hard with the regions of the radius-mass plane al- 
lowed if i'qpo2 is a Keplerian orbital frequency, for 
nonrotating stars and three values of j^qpo2- order 
to make possible comparisons with previous studies 
of neutron star properties (see, e.g, Pethick & Raven- 
hall 1995), we show the mass-radius curve given by 
the early Friedman-Pandharipande-Skyrme (FPS) re- 
alistic equation of state (Friedman & Pandharipande 
1981; Lorenz, RavenhaU, & Pethick 1993). The FPS 
equation of state uses a different approach, but is 
similar to the softest equations of state permitted by 
modern realistic models of the nucleon-nucleon inter- 
action (Pandharipande, Akmal, & RavenhaU 1998). 
These equations of state all give maximum gravita- 
tional masses of about 1.8 Af© for nonrotating stars. 



As an example of the mass-radius curves given by 
later realistic equations of state, we show the mass- 
radius curve predicted by the UU equation of state 
(Wiringa, Fiks, & Fabrocini 1988). Although it is 
based on older scattering data, the UU equation of 
state is similar to the recent A18 -f- UIX' -t- Svi, equa- 
tion of state (Akmal, Pandharipande, & RavenhaU 
1998), which is based on the most modern scattering 
data. Like the A18 -I- UIX' + 5v\^ equation of state, 
the UU equation of state gives a maximum mass of 
about 2.2 Mq for a nonrotating neutron star. 

As an example of the mass-radius curves predicted 
by the relatively stiff equations of state typically given 
by mean field theories, we include the mass-radius 
curve for the mean-field equation of state of Pand- 
haripande & Smith (1975b; L in the Arnctt & Bow- 
ers [1977] survey). The maximum mass of a nonro- 
tating star constructed using equation of state L is 
2.7M0. We also show the mass-radius curve given by 
the very early tensor interaction (TI) equation of state 
of Pandharipande & Smith (1975a; M in the Arnett & 
Bowers [1977] survey) and the Reid soft-core equation 
of state of Pandharipande (1971; A in the Arnett & 
Bowers [1977] survey). The maximum masses of non- 
rotating stars constructed using equations of state M 
and A are 1.8 Mq and 1.65 Mq, respectively. Equa- 
tions of state A, L, and M are no longer of interest 
to nuclear physicists and arc included here primarily 
to facilitate comparison with previous work on mass- 
radius constraints. 

For a particular source, QPO frequency, and equa- 
tion of state, the allowed portion of the mass-radius 
relation is the segment within the pie-slice shaped re- 
gion analogous to the region in Figure |l^ outlined 
by the heavy line or, if the QPO frequency is identi- 
fied as the frequency of the innermost stable circular 
orbit, simply the horizontal heavy line. If the mass- 
radius relation given by a particular equation of state 
intersects the allowed region defined by the highest 
QPO frequency seen in a given source, that partic- 
ular equation of state is viable. If that equation of 
state is furthermore the correct equation of state, the 
mass and radius of the neutron star in the source must 
correspond to one of the points along the segment of 
the mass-radius relation that lies within the pie-slice 
shaped region, so the mass and radius of the star are 
bounded from above and from below; for most equa- 
tions of state, only a narrow range of radii is allowed. 

For example, the highest-frequency QPO so far 
seen in any source is the 1220 Hz QPO observed 
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in 4U 1636—536 (W. Zhang, personal communica- 
tion). Hence, if equation of state M is the correct 
equation of state, then the mass and radius of the 
neutron star in 4U 1636—536 would have to satisfy 
1.7 Mq < M < 1.8 Mq and 11.6 km < i? < 16.0 km 
if it were nonrotating. If equation of state M is the 
correct equation of state and 1220 Hz is identified 
as the frequency of the innermost stable circular or- 
bit, then the mass and radius of the neutron star in 
4U 1636—536 would be determined as about 1.8 Mq 
and between 11 km and 12 km, if the star were not 
rotating. 

An equation of state that gives a mass-radius re- 
lation that does not intersect the region allowed for 
a given source is ruled out for that source. We stress 
that because the equation of state of the matter in 
neutron stars is expected to be essentially the same 
in all such stars, an equation of state that is incon- 
sistent with the properties of any neutron star is ex- 
cluded for all neutron stars. For example, observation 
of a 1500 Hz Keplerian orbital frequency in any source 
would rule out equations of state L and M. 

So far our discussion of constraints on the neutron 
stars in the kilohertz QPO sources has assumed that 
they are not rotating. As we discuss in the next sub- 
section, the changes in these constraints caused by the 
spin of the star are likely to be small in most sources. 

5.2. Effects of Stellar Rotation 

Rotation affects the structure of the star for a given 
mass and equation of state and the spacetime exte- 
rior to the star. Hence the mass-radius relation, the 
orbital frequency at a given radius, and the radius 
Rms of the innermost stable circular orbit are all af- 
fected. As a result, the bounds on the mass and ra- 
dius of the star and the constraints on the equation 
of state implied by observation of a kilohertz QPO 
of a given frequency are affected. Obviously, the size 
of these effects depends on the spin rate of the star. 
Treating these effects accurately will be particularly 
important if the frequency of a kilohertz QPO is ever 
securely identified as the Keplerian frequency at the 
marginally stable orbit around a neutron star (see 
§ 5.4), because the ability to rule out some equations 
of state may depend on it. 

The parameter that characterizes the importance 
of rotational effects is the dimensionless quantity 
j = cJ/GM'^, where J and M are the angular mo- 
mentum and gravitational mass of the star. The 



value of j that corresponds to a given observed spin 
frequency depends on the neutron star mass and 
equation of state, and is typically higher for lower 
masses and stiffer equations of state. For example, 
at the 363 Hz spin frequency inferred for 4U 1728—34 
(Strohmayer et al. 1996c), a 1.4 M0 neutron star with 
the softer equation of state A has j ~ 0.1, whereas a 
1.4 neutron star with the very stiff equation of 
state L has j « 0.3 (see Table ^). 

As discussed in § 5.1, if the frequency of the higher- 
frequency QPO in a kilohertz QPO pair is an or- 
bital frequency, then the region of the radius-mass 
plane allowed for the neutron star in that source is 
bounded below by the curve i?orb(A/), which gives 
the orbital radius at which the Keplerian frequency is 
equal to the highest observed frequency '^qpo2 of the 
higher-frequency QPO. Given j, it is straightforward 
to compute the orbital radius for any given orbital 
frequency as a function of M. However, it is not j 
but instead the star's spin frequency I'spin that is de- 
termined by the observations. Hence, to obtain the 
relevant i?oib(A/) curve one must vary j in such a 
way that i^spin is kept constant as M is varied. More- 
over, j depends not only on I'spin and AI but also 
on the equation of state. Thus the relevant i?oib(M) 
curve depends not only on the star's spin rate, but 
also on the equation of state assumed. Therefore, for 
rotating stars (unlike static stars), one cannot present 
a single i?orb(-M) curve that constrains all equations 
of state, even for a star with a given spin frequency. 

As also discussed in § 5.1, the region of the radius- 
mass plane allowed for the neutron star in a given 
source is bounded above by the horizontal line M = 
Mmax, which is the mass at which the curve i?orb(Af) 
intersects the curve Rais{M); the radius R of the star 
is bounded above by i?orb(Mmax)- The two curves 
Rms{M) and Roih{M) are both affected by rotation 
of the star, so Mmax depends on the stellar spin rate. 

Determining the mass at which i?orb(Af) = Rms{M) 
is equivalent to determining the mass at which VK{Rms) 
is equal to i'qpo2- Here we focus on the latter condi- 
tion. The allowed mass of a rotating star is bounded 
above if, as M is increased at constant t'spin, the or- 
bital frequency vxiRras) crosses the QPO frequency 
Vqpo2 from above, once and only once. As we show 
below, this is the case for slowly rotating stars. How- 
ever, this may not be the case for some neutron star 
equations of state and spin rates. (For a Kerr black 
hole with fixed spin frequency, fK(-Rms) first decreases 
and then increases with increasing mass.) 
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Fortunately, most of the neutron stars that ex- 
hibit kilohertz QPOs appear to have spin frequen- 
cies in the range 250-350 Hz (see § 4.2), which is 
low enough that a first-order treatment of rotational 
effects is adequate. To see this, note that neutron 
stars with masses ^ 1.5-2 have moments of in- 
ertia / ^ (1-3) X 10''^ gcm~^ , so stars with spin fre- 
quencies < 350 Hz have j-values < 0.3 (see Table |). 
The lowest-order changes in the structure of a rotat- 
ing star are O(j^) so, to first order in j, the mass- 
radius relation and moment of inertia of a rotating 
star are the same as for a nonrotating star of the 
same mass (Hartle & Thorne 1968). Thus the error 
made by neglecting higher-order terms is < 10% for 
spin frequencies < 350 Hz . 

We now compute the mass and radius constraints 
imposed on a slowly rotating neutron star by observa- 
tion of a kilohertz orbital frequency. The calculation 
is simplified by the fact that to first order in j, the 
spacetime outside a uniformly and steadily rotating 
relativistic star is the same as the Kerr spacetime for 
the same M and j (Hartle & Thorne 1968). 

Consider first the effect of stellar rotation on the 
motion of gas orbiting the star. To first order in j, 
the orbital frequency (measured at infinity) of gas in 
a prograde Keplerian orbit at a given Boyer-Lindquist 
radius r is 

i^Kir, M,j) « [1 - j(GAf/rc2)3/2] „o^{r, M) (32) 
and the radius of the marginally stable orbit is 

RUM, 3) « [1 - j(2/3)'/'] i?L(M) , (33) 

where and are the Keplerian frequency and 
radius of the marginally stable orbit for a nonrotat- 
ing star. Expressions (|3^) and (|3|) are first-order 
expansions of the exact expressions for these quanti- 
ties given by Bardeen, Press, & Teukolsky (1972) for 
the Kerr spacetime. Hence, to first order in j, the fre- 
quency of the prograde orbit at i?ms around a star of 
given mass M and dimensionless angular momentum 
j is (cf. Kluzniak, Michelson, & Wagoner 1990) 



l-j(l/6)^/^ 



1 



(2/3)1/2 1^1,^^ 



2210 (1 4^ O.75j)(Af0/Af) Hz 



(34) 



where Vy^ is the Keplerian orbital frequency at the 
radius of the innermost stable circular orbit for a non- 
rotating star. Therefore, the net effect of the star's 



rotation is to increase the frequency of the prograde 
orbit at Rms- 

As explained above, the region of the radius-mass 
plane allowed for a given star is not the region allowed 
for constant j but is instead the region allowed for 
constant I'spin- To demonstrate that the mass and ra- 
dius of the star are bounded from above, it is sufficient 
to show that VK{Rras) decreases with increasing M at 
fixed Vspim or equivalently, that (dz^K,ms/rf-Af)i/api„ < 
for all M for the equation of state under considera- 
tion. Now to first order in j. 
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(36) 



Note that the derivatives on the right side of equa- 
tion ( |35| ) are to be evaluated at j = 0. In order to 
show that (rfi^K,ms/'^-Af)i/a i„ is negative for slowly ro- 
tating stars, it is sufficient (but not necessary) to show 
that (dlnl/dln Af)i/^ is always less than 2. This is 
the case for all the equations of state tabulated by 
Cook et al. (1994) and is the case even for incompress- 
ible matter, which is unphysically stiff. (For a star 
made of incompressible matter, / cx Afi?^ cx M^^^, 
so (dlnl/dlnM) = | < 2.) Thus, for slowly rotat- 
ing stars VKiR-ms) decreases with increasing mass for 
constant t'spin, so the masses and radii of such stars 
are bounded above. 

Computation of the mass and radius constraints is 
straightforward but depends on the stellar spin rate 
and equation of state. The upper bounds on the mass 
and radius are given implicitly by 



and 



M„ 



Rn 



[l + 0.75j(i^,pi„)]A/0, 



(37) 



« [1 + 0.20j(i^spi„)]<ax , (38) 

where J(^'spin) is the value of j for the observed stellar 
spin rate at the maximum allowed mass for the equa- 
tion of state being considered and M^^^ and i?^ax 
are the maximum allowed mass and radius for a non- 
rotating star (see eqs. |3^ and |3|]). Expressions (37) 
and ( |3^ ) show that the bounds are always greater for 
a slowly rotating star than for a nonrotating star, re- 
gardless of the equation of state assumed. 
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Figure ^ illustrates the effects of stellar rotation 
on the region of the radius-mass plane allowed for a 
given star for spin rates ~ 300 Hz , like those inferred 
for the kilohertz QPO sources, and i^Qpo2 = 1220 Hz , 
the frequency of the highest-frequency QPO so far 
observed in 4U 1636—536, which is also the highest- 
frequency QPO so far observed in any source. Our 
calculations show that the mass of the neutron star in 
4U 1636—536 must he less than ~ 2.2 and its ra- 
dius must be less than ~ 17 km. As explained above, 
the precise upper bounds depend on the equation of 
state assumed. 

For rapidly rotating stars, j is not small compared 
to unity and the structure of the star depends appre- 
ciably on its rotation rate. Derivation of bounds on 
the mass and radius of a given star for an assumed 
equation of state therefore requires construction of a 
sequence of stellar models and spacetimes for differ- 
ent masses using the assumed equation of state, with 
Jospin as measured at infinity held fixed. The maxi- 
mum and minimum possible masses and radii allowed 
by the observed QPO frequency can then be deter- 
mined. 

5.3. Effects of the Radial Radiation Force 

The luminosities of the Z sources are typically 
~ 0.5-1 Lb, where Le is the Eddington luminosity 
(see § 2). Hence, in the Z sources the outward accel- 
eration caused by the radial component of the radia- 
tion force can be a substantial fraction of the inward 
acceleration caused by gravity. Therefore, in these 
sources the radial component of the radiation force 
must be included in computing the Keplerian orbital 
frequency near the neutron star and taken into ac- 
count when constraints on the mass and radius of 
the star are derived using the procedures discussed 
in § 5.1 and § 5.2. 

The radially outward component of the radiation 
force reduces the orbital frequency at a given radius. 
For example, if the star is spherical and nonrotating 
and emits radiation uniformly and isotropically from 
its entire surface, the orbital frequency (measured at 
infinity) of a test particle at Boyer-Lindquist radius r 



(1 - 3GAf/rc2)i/2 L 1 ^'^ 



(l-2GM/rc2) Le 



sence of radiation forces. Thus, the Boyer-Lindquist 
radius of a circular orbit with a given frequency is 
smaller in the presence of the radial radiation force 
and the constraints on the mass and radius of the 
star are therefore tightened. 

For the atoll sources, which have luminosities 
L < 0.1 Lb, the change in the Keplerian frequency is 
at most ~ 5%. For the Z sources, on the other hand, 
which have luminosities L fa Le, the change may be 
much larger, although the change in the sonic-point 
Keplerian frequency may be smaller than would be 
suggested by a naive application of equation (|39|), if 
a substantial fraction of the radiation produced near 
the star is scattered out of the disk plane before it 
reaches the sonic point. A more detailed analysis of 
the effect of the radial radiation force on the con- 
straints on the mass and radius of the star will be 
reported elsewhere. 

5.4. Signatures of the Innermost Stable 
Circular Orbit 

As explained in § 4.2, observations of kilohertz 
QPO sources may be able to demonstrate the ex- 
istence of an innermost stable circular orbit around 
some of the neutron stars in these sources. If so, this 
would be the first confirmation of a strong-field pre- 
diction of general relativity and an important step 
forward in our understanding of strong-field gravity. 

If the frequency of a kilohertz QPO produced by 
an LMXB can be securely established as the orbital 
frequency at the radius of the innermost stable cir- 
cular orbit and the spin frequency of the star can 
be determined, the mass of the neutron star in that 
source can be determined for each assumed equation 
of state. Depending on the range of allowed equations 
of state and masses, this could have profound conse- 
quences for our understanding of the equation of state 
of neutron stars.]] For example, establishing a mass of 
2.0 M0 for a slowly rotating neutron star would rule 
out eight of the twelve currently viable equations of 
state considered by Cook et al. (1994). 

Given the profound consequences that would fol- 
low from identifying the frequency of a kilohertz QPO 



(39) 

where L is the luminosity of the star measured at in- 
finity and t^K (0) is the Keplerian frequency in the ab- 



^FoUowing submission and circulation of the original version of 
the present paper, Kaaret, Ford, & Chen (1997) and Zhang, 
Strohmayer, & Swank (1997) considered the implications if the 
frequencies of some of the kilohertz QPOs that have already 
been observed are the orbital frequencies of innermost stable 
circular orbits. 
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with the orbital frequency of the innermost stable or- 
bit in any source, it is very important to establish 
what would constitute strong, rather than merely sug- 
gestive, evidence of such a detection. Observations 
that would signal detection of the innermost stable 
orbit include the following: 

QPO frequency signature. — The strongest evidence 
that a QPO with the orbital frequency of the in- 
nermost stable orbit has been detected would be re- 
producible observation of a fairly coherent, kilohertz 
QPO with a frequency that first increases steeply with 
accretion rate but then, at a high frequency, becomes 
nearly constant as the accretion rate continues to in- 
crease (see § 3.3, § 4.2, and Fig. 9). 

The QPO frequency will approach a constant be- 
cause the sonic point in the flow cannot retreat closer 
to the star than the radius of the marginally stable 
orbit and coherent QPOs are not expected from the 
rapidly inspiraling gas inside the sonic point. Al- 
though the count rate above which the frequency be- 
comes constant may vary, the constant frequency it- 
self should not vary, because the orbital frequency of 
the marginally stable orbit depends only on the mass 
and rotation rate of the neutron star, which remain 
almost constant over many years. 

In the sonic-point model, the frequency of the 
lower-frequency QPO in a kilohertz QPO pair is the 
beat frequency and therefore should also increase 
steeply with accretion rate at first but then become 
approximately constant at the same accretion rate at 
which the frequency of the higher-frequency QPO be- 
comes constant. 

QPO am.plitude signature. — A second signature that 
the marginally stable orbit has been detected would 
be reproducible observation of a decrease in the am- 
plitude of the lower-frequency QPO in a kilohertz 
QPO pair or a simultaneous decrease in the ampli- 
tudes of both QPOs in a pair, at a QPO frequency 
that is always the same in a given source. 

The amplitude of the QPO at the beat frequency is 
expected to decrease once the sonic point has moved 
inward to the marginally stable orbit because in the 
sonic-point model, the QPO at the beat frequency is 
generated by the drag force exerted by the radiation 
coming from the stellar surface that reaches the sonic 
point. Therefore, a strong QPO is expected at the 
sonic-point beat frequency only if the drag force ex- 
erted by the radiation is dynamically important at the 
sonic point. If, however, the sonic point has moved in- 



ward to the innermost stable orbit and remains there 
as the accretion rate continues to increase, the opti- 
cal depth from the stellar surface to the sonic point 
will continue to rise and radiation drag will become 
less and less important there, causing modulation of 
the inflow from the sonic point by the radiation force 
to weaken. Thus, a decrease in the amplitude of the 
QPO at the sonic-point beat frequency with increas- 
ing accretion rate would signal the approach of the 
sonic-point to the radius of the innermost stable or- 
bit (see § 3.1). 

If radiation forces play an important role in creat- 
ing or amplifying the clumps that produce the QPOs 
(see § 3.1), the amplitude of the Keplerian-frequency 
QPO may also decrease with increasing accretion 
rate, once the radius of the sonic point has reached 
the radius of the innermost stable orbit (see § 3.2). 

Possible QPO coherence signature. — A possible sig- 
nature that the orbit of the gas that is generating the 
QPO has receded inside the radius of the marginally 
stable orbit would be a steep drop in the coherence 
of both QPOs in a kilohertz QPO pair (or in the co- 
herence of the Keplerian-frequency QPO, if the beat- 
frequency QPO is not visible), at a certain critical 
frequency, as the frequencies of the QPOs increase 
steadily with accretion rate. This would occur if ra- 
diation forces are able to generate clumping in the 
flow at a radius inside the sonic point. If this occurs 
at all, clumps are likely to be produced at a range of 
radii and to last only a short time, so any oscillations 
that may be generated are likely to have low coher- 
ence (see § 3.4). The critical frequency would be the 
orbital frequency at the marginally stable orbit and 
hence should always be the same in a given source. 

Approach of the sonic radius to the radius of the 
innermost stable orbit can be distinguished from ap- 
proach of the sonic radius to the radius of the star be- 
cause coherent Keplerian- and beat-frequency QPOs 
may continue in the fist case but are very unlikely in 
the second. Kilohertz QPOs are unlikely to be gen- 
erated if the sonic point moves close to the stellar 
surface because of the disruptive effect of the stellar 
magnetic field and the the viscous shear layer that 
is expected to develop if the Keplerian flow interacts 
directly with the stellar surface (see § 4.2). 

6. SUMMARY AND CONCLUSION 

The sonic-point model explains naturally the most 
important features of the kilohertz QPOs observed 
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in the atoll and Z sources, including their frequen- 
cies, large amplitudes, and high coherence. It also 
explains the frequent occurrence of two simultaneous 
QPOs, the observed steep increase of kilohertz QPO 
amplitudes with increasing photon energy in the ~ 5- 
10 keV energy range, and the anticorrelation of kilo- 
hertz QPO amplitudes with the strength of the stellar 
magnetic field inferred from spectral models. An at- 
tractive feature of the sonic-point model is that the 
magnetic fields, accretion rates, and scattering optical 
depths that it requires are completely consistent with 
those inferred previously from observations and mod- 
eling of the X-ray spectra and lower-frequency X-ray 
variability of the kilohertz QPO sources. 

The sonic-point model leads to several general ex- 
pectations about the kilohertz QPOs: 

(1) Kilohertz QPOs with the properties seen in the 
atoll and Z sources should not be observed in black- 
hole LMXBs, because collision of the accretion flow 
with the stellar surface plays an essential role in the 
sonic-point model. 

(2) The rms amplitudes of the kilohertz QPOs 
should be anticorrelated with the strength of the neu- 
tron star's magnetic field. From this it follows that 
kilohertz QPOs should be very weak or undetectable 
with current instruments in the GX group of atoll 
sources and in the Cyg-like Z sources (which in- 
clude CygX-2, GX5-1, and GX 340-hO), and un- 
detectable with current instruments in any sources 
that have strong magnetic fields and therefore pro- 
duce strong, periodic oscillations at their spin fre- 
quencies. (After this specific prediction was made 
in the original version of this paper, kilohertz QPOs 
with very low amplitudes were detected in all six of 
the originally identified Z sources; see § 4.4.) Another 
consequence of expectation (2) is that the positive 
correlation of rms amplitude with hard color shown 
by the four sources plotted in Figure |l^ should be 
found to be general among sources showing kilohertz 
QPOs. 

(3) // kilohertz QPOs are detected in the tails of 
X-ray bursts, their frequencies will tend to be lower 
when the accretion rate is higher. 

(4) Weak oscillations should eventually be detected 
at overtones of the sonic-point beat and Keplerian fre- 
quencies and perhaps at the stellar spin frequency, but 
oscillations at other frequencies, such as t^Ks + t'spiiu 
should be extremely weak. 

At a lower level of certainty, the sonic-point model 



suggests that (1) the dependence of QPO amplitude 
on photon energy will change as the X-ray spectrum 
of the source changes, (2) the amplitude of the higher- 
frequency QPO in a kilohertz QPO pair will drop rel- 
ative to the amplitude of the lower-frequency QPO at 
high luminosities, and (3) either the lower-frequency 
QPO or the higher-frequency QPO in a pair may be 
undetectable even when the other QPO is detectable. 

In closing, we emphasize that measurement of Ke- 
plerian frequencies in the kilohertz range provides in- 
teresting new upper bounds on the masses and radii 
of the neutron stars in the kilohertz QPO sources and 
important new constraints on the equation of state of 
the matter in all neutron stars, as we have shown in 
§ 5. As we demonstrated there, if the neutron star in 
the atoll source 4U 1636—536 has a spin frequency of 
^ 290 Hz, as indicated by the difference between the 
spin frequencies of its two high-frequency QPOs, the 
1220 Hz QPO observed in this source constrains its 
mass to be less than about 2.2 Mq and its radius to 
be less than about 17 km; the precise bounds depend 
on the equation of state assumed. 

If at some future time we are able to demonstrate 
the existence of an innermost stable circular orbit 
around one or more neutron stars using observations 
of kilohertz QPOs, this would be the first confirma- 
tion of a prediction of general relativity in the strong- 
field regime and a major advance in our understand- 
ing of strong-field gravity. If the frequency of a kilo- 
hertz QPO is securely established as the orbital fre- 
quency at the radius of the innermost stable circular 
orbit in a particular system by, for example, observing 
the signatures discussed in § 5, and if, in addition, the 
spin frequency of the neutron star can be determined, 
then the frequency of the QPO will fix the mass of 
the neutron star for each assumed equation of state, 
providing a better understanding of the properties of 
dense matter and possibly ruling out many currently 
viable equations of state. 
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Fig. 1. Schematic drawings of the accretion flows 
expected around atoll and Z sources. Top: Side view 
of an atoll source with a dipole magnetic field ~ 10^ G 
and a mass accretion rate ~ 0.01 Me- The star's spin 
axis and the direction of its magnetic moment are in- 
dicated by the arrows labeled O and ji. Gas in the 
inner part of the Keplerian disk is clumpy and pen- 
etrates very close to the stellar surface before some 
of it is channeled by the star's magnetic field. Colli- 
sion of the channeled gas with the stellar surface pro- 
duces slightly brighter spots that rotate with the star. 
The light shading indicates the hot gas in the mag- 
nctosphcre and the central corona that surroimds the 
neutron star. Bottom: Side view of a Z source with 
a magnetic field ~ 10^-5 x 10^ G and a mass accre- 
tion rate ^ Me , showing how some of the gas in the 
disk is channeled out of the disk by the star's mag- 
netic field at ~ 2-3 stellar radii but some continues to 
flow inward in a Keplerian disk flow that penetrates 
close to the stellar surface. Collision of the channeled 
gas with the stellar surface produces brighter spots 
that rotate with the star. The light shading again 
indicates the hot gas in the magnetosphere and the 
central corona that surrounds the neutron star. The 
arrows indicate the cooler, approximately radial in- 
flow outside the central corona that is present in the 
Z sources. 

Fig. 2. — Parameter space of dipole magnetic mo- 
ment fi and mass accretion rate (in units of the Ed- 
dington mass accretion rate Me that produces an Ed- 
dington luminosity) for neutron star LMXBs. The 
axis on the right shows the equivalent surface mag- 
netic field at the pole for a neutron star with a radius 
of 10 km. The thin solid lines, plotted using equa- 
tion (1), indicate the magnetic moment for which, at 
the given accretion rate, the Keplerian frequency is 
500 Hz (top line) and 1100 Hz (bottom line) where 
a substantial fraction of the gas in the disk first 
couples to the stellar field and is channeled out of 
the disk. The break at M = 0.3 Me illustrates the 
effect of the expected transition from gas-pressure- 
dominated fiow at low accretion rates to radiation- 
pressure-dominated flow at high accretion rates. The 
lines in this M range are dotted to indicate that they 
are not accurate. The solid lines, plotted using equa- 
tion (2) for electron temperatures kT^, = 5 keV (top 
line) and fcTg = 15 keV (bottom line), divide the M- 
Bd plane into the regions where the dominant source 
of soft photons is cyclotron emission in the neutron 
star magnetosphere (above the line) or blackbody 



emission from the neutron star surface (below the 
line). The labels "4U", "GX", and "Z" indicate the 
approximate regions in the M-B^ plane occupied, re- 
spectively, by 4U atoll sources, GX atoll sources, and 
Z sources. These curves were calculated assuming a 
neutron star mass M = 1.4 Mq, an effective cyclotron 
photosphere of radius i?cyc = 1-5 x 10® cm, a neutron 
star radius of 10® cm, and a spectrum truncated at 
cyclotron harmonic number n = 15, appropriate for a 
scattering optical depth t — 6. For masses, radii, or 
optical depths different from these values, the curves 
can shift by as much as 50% in magnetic field at a 
given accretion rate. 

Fig. 3. — A neutron star with a dynamically neg- 
ligible magnetic field, accreting via a Keplerian disk 
that penetrates close to the star. The star and disk 
are viewed along the rotation axis of the disk, which 
is rotating counterclockwise in this view. The pan- 
els show, in Boyer-Lindquist r, (f) coordinates, (a) the 
spiral trajectory followed by a single element of gas 
as it falls supersonically from the sonic radius to the 
stellar surface and (b) the spiral pattern of higher gas 
density formed by gas streaming inward along spiral 
trajectories with the shape shown in panel (a) from 
a region of denser gas (a "clump" ) orbiting near the 
sonic radius. The gas trajectories and resulting den- 
sity pattern were computed in full general relativity 
assuming that the gas is exposed to radiation from the 
star at a radius of 9M and that the star has a radius 
of 4M, is nonrotating, radiates isotropically, and has 
a luminosity measured at infinity of 0.005 i^- The 
surface density of the disk flow is much smaller in- 
side the sonic radius (lighter shaded region) than out- 
side (darker shaded region), because of the sharp in- 
crease in the inward radial velocity at the sonic radius. 
Gas falling inward from the sonic radius along spi- 
ral trajectories collides with the neutron star around 
its equator, producing an X-ray emitting equatorial 
ring, which is indicated by the grey ring around the 
star. The white arc at the stellar surface indicates 
the bright, arc-shaped "footprint" where the denser 
gas from the clump collides with the stellar surface 
and produces a beam of X-rays (white dashed lines) 
that rotates around the star at the sonic-point Kep- 
lerian frequency. The footprint generally moves with 
respect to the stellar surface. 

Fig. 4. — Time sequence of four snapshots of a neu- 
tron star with a dynamically negligible magnetic field 
accreting via a Keplerian disk that penetrates close 
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to the star, showing schcmaticany how the QPO at 
the sonic-point Keplerian frequency is produced. The 
viewpoint, gas trajectories, and density pattern are 
the same as in Fig. 3. The drop in the surface density 
of the disk flow at the sonic radius, the bright ring 
of emission around the stellar equator, the pattern of 
denser infalling gas inside the sonic radius, the bright 
footprint where the denser gas collides with the stel- 
lar surface, and the beam of radiation coming from 
the footprint are all indicated in the same way as in 
Fig. 3. The clump is shown advancing in its orbit by 
90° from one panel to the next. The spiral pattern 
of higher density gas, its bright footprint, and the 
resulting beam of X-rays all rotate around the star 
with a frequency equal to the orbital frequency of the 
clump, so they are also shown advancing by 90° from 
one panel to the next. The footprint generally moves 
with respect to the stellar surface. 

Fig. 5. — Time sequence of six snapshots of a spin- 
ning neutron star with a weak magnetic field accreting 
via a prograde Keplerian disk that penetrates close 
to the star, showing schematically how the QPO at 
the sonic-point beat frequency is produced. The disk 
and the star are viewed along their common rotation 
axes and are rotating counterclockwise in this view. 
The beam of radiation produced by the collision with 
the stellar surface of the gas channeled by the weak 
stellar magnetic field is indicated by the faint white 
beam emerging from behind the star; this beam ro- 
tates with the star. To make the sequence of events 
easier to visualize, the snapshots show the sequence of 
events as seen in a frame corotating with the star, so 
the beam of radiation that rotates with the star is al- 
ways pointing in the same direction (here, the vertical 
direction). The angular velocity of the clump, which 
is orbiting near the sonic radius, is greater than the 
angular velocity of the star, so with time the clump 
advances counterclockwise relative to the beam of ra- 
diation that rotates with the star. In this sequence, 
the clump advances relative to the star by 90° from 
one snapshot to the next. The drop in the surface 
density of the disk flow at the sonic radius, the bright 
ring of emission around the stellar equator, the pat- 
tern of denser infalling gas inside the sonic radius, the 
bright footprint where the denser gas collides with 
the stellar surface, and the beam of radiation coming 
from the footprint are all indicated in the same way 
as in Fig. 3. The events occurring in each panel are 
described in the text. 



Fig. 6. — Solid curve: Angular momentum per unit 
mass Ufj, (as measured at infinity in units of the grav- 
itational mass M of the star) of an element of gas 
in a circular Keplerian orbit around a nonrotating 
neutron star, as a function of the Boyer-Lindquist ra- 
dial coordinate of the orbit (in units of M). Dashed 
curve: for a rotating star with dimensionless an- 
gular momentum j =0.1, which is close to the j-value 
given by the best modern neutron-star matter equa- 
tions of state for a neutron star with M = 1.4 Mq and 
a spin frequency of 300 Hz as measured at infinity (see 
§ 5.2). The top axis shows the radius in kilometers for 
M = IAMq. The radius i?ms of the innermost stable 
circular orbit is the radius at which the specific angu- 
lar momentum is a minimum; for j = 0, Rms = 6M, 
whereas for j = 0.1, i?ms = 5.7 M. 

Fig. 7. — Solid curve: Fraction rjdow of the specific 
angular momentum of an element of gas in Keple- 
rian circular orbit at Boyer-Lindquist radial coordi- 
nate r (measured in units of the stellar mass M) that 
must be removed in order for the gas to fall from r 
to the radius Rms of the innermost stable circular or- 
bit. The curve shown is for a static star {j = 0); the 
curves for slowly rotating stars would be little differ- 
ent. Dashed curve: Estimate of the largest fraction 
^rad(inax) of the specific angular momentum of an el- 
ement of gas that can be removed by radiation coming 
from the surface of a nonrotating, isotropically radi- 
ating, spherical star of radius 5M (see text). The top 
axis shows the radius in kilometers for M = IAMq. 

Fig. 8. Numerical results for the structure of 
the accretion flow near the neutron star in the fully 
general relativistic model of gas dynamics and radia- 
tion transport in the inner disk described in the text, 
(a) The inward radial velocity of the gas in the 
disk measured by a local static observer, (b) The an- 
gular velocity of the gas in the disk measured by an 
observer at infinity, (c) The surface density Tt^o of the 
disk flow measured by an observer comoving with the 
flow, (d) The radial optical depth from the stel- 
lar surface through the disk flow to the radius shown 
on the horizontal axis. The radial coordinate is the 
Boyer-Lindquist radius in units of the stellar mass M; 
the radius of the star is 5M. The four curves in each 
panel are labeled with the assumed accretion rate Mi 
through the inner disk, measured in units of the ac- 
cretion rate Me that would produce an accretion lu- 
minosity at infinity equal to the Eddington critical 
luminosity. 
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Fig. 9. Sonic-point Kcplcrian and beat frc;qucncic;s 
in Hertz as functions of the accretion luminosity in 
units of the Eddington critical luminosity, computed 
using the accretion flow model of § 3.2. The neutron 
star spin rate is assumed to be 300 Hz. The other 
assumptions are the same as in Fig. 8. The sonic- 
point Keplerian frequency i/ks (solid line) increases 
steeply with increasing accretion luminosity until it 
reaches z^(iims) (dashed horizontal line), the orbital 
frequency at the innermost stable circular orbit, at 
which point stops increasing. The sonic-point 
beat frequency i^bs (dotted line) increases steeply 
with increasing accretion luminosity until it reaches 
^'K(^^ms) — t'spin, at which point it too stops increasing. 

Fig. 10. — Sample attenuation factors Aoc/Aq for 
the relative amplitudes of the X-ray brightness oscil- 
lations produced by a hypothetical accreting neutron 
star with a spin rate i/gpin = 300 Hz and a sonic-point 
Keplerian frequency uks = 1100 Hz , at the center of 
a uniform, spherical scattering cloud with a radius of 
3 X 10^ cm. The sonic-point beat frequency z^bs is 
assumed to be j^ks — t^spin and is therefore 800 Hz. 
The attenuation factors for the oscillations at i^ks 
(dashed line) and at the upper sideband frequency 
i^Ks + Jospin = 1400 Hz (dotted line), which are as- 
sumed for the sake of illustration to be pure beaming 
oscillations produced by a rotating pencil beam, are 
>5 for scattering optical depths r>5. In contrast, 
the attc;nuation factor for the oscillation at i^bs (solid 
line), which is a pure luminosity oscillation, is much 
smaller. 

Fig. 11. — Measured amplitudes of the high- 
frequency QPOs seen in 4U 1608—52 (open circles) 
and 4U 1636—536 (filled circles) as a function of pho- 
ton energy. Data for 4U 1636—536 were kindly pro- 
vided by W. Zhang (1997, personal communication) 
and reflect corrections made after the report by Zhang 
et al. (1996) was published. The dotted curve shows 
the variation of the rms amplitude with photon en- 
ergy due to the optical depth variations expected in 
the sonic-point model (see text). This curve is not 
a fit; instead, for illustrative purposes we have as- 
sumed (consistent with the unified model of neutron 
star LMXBs) that the input spectrum is a blackbody 
of temperature kT = 0.6 keV and that the Comp- 
tonizing corona has a temperature kT = 9 keV and 
an optical depth that varies from r = 3 to t = 3.3. 

Fig. 12. — Measured rms amplitudes of the high- 



frequency QPOs seen in Sco X-1, 4U 1735-444, 
4U 1636-536, and 4U 0614+091 plotted against the 
PCA X-ray hard colors (defined as the ratio of the 
counts in the 7-20 kcV bin to the counts in the 5 7 
keV bin) of these sources. The photon energy ranges 
used in computing these amplitudes were 2-20 keV for 
Sco X-1 (van der Klis et al. 1996d), 0.24 18.37 kcV 
for 4U 1735-444 (Wijnands et al. 1996), 7-20 keV 
for 4U 1636-536 (van der Klis et al. 1996c), and 2- 
20 keV for 4U 0614-^091 (van der Klis et al. 1996c). 
The observed steep increase of QPO amplitudes with 
photon energy (Fig. 11) means that a direct, quan- 
titative comparison of QPO amplitudes in different 
sources can only be made if they are all computed 
using the same photon energy ranges. For example, 
the QPOs in GX 5—1 arc not detected below 10 keV, 
but between 10 keV and 50 keV the QPO amplitude 
is as high as ~7% (van der Klis et al. 1996e). It is 
not straightforward to compare this with the reported 
~1% amplitude for Sco X-1 in the 2-20 keV band. 
The range of hard colors for Sco X-1 is that near the 
soft vertex. The range of hard colors for 4U 1735—444 
has been artificially extended to reflect possible sys- 
tematic errors in the PCA response caused by gain 
changes. In neutron star LMXBs, the X-ray hard 
color is found to be greater for sources with weaker 
magnetic fields (Psaltis et al. 1995; Psaltis & Lamb 
1998), so the correlation evident in this figure is strik- 
ing confirmation that the amplitudes of these QPOs 
are lower for sources with stronger magnetic fields. 
The PCA colors were kindly provided to us by Rudy 
Wijnands. 

Fig. 13. — Radius-mass plane, showing how bounds 
on the mass and radius of a nonrotating neutron 
star in a QPO source with fQpo2 = 1220 Hz can 
be constructed. The dashed curve that begins at 
the origin and rises up and to the right is the re- 
lation M - M"(7?orb: fQP02)' foi" ''QP02 = 1220 Hz ; 
this curve gives the mass for which the orbital fre- 
quency at the radius shown on the horizontal axis 
is 1220 Hz. The diagonal dotted line is the re- 
lation M = M°(i?ms); this line gives the mass for 
which the radius of the innermost stable circu- 
lar orbit is equal to the radius shown on the hor- 
izontal axis. The horizontal solid line is the mass 
M^^^ at which M°{R^s) is equal to M°(i?orb, i^QP02) 
and hence R'^^^M) is equal to ^orb(-^' ^qp02)> for 
^QP02 = 1220 Hz . The hatched region shows the 
combinations of stellar mass and radius that are 
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excluded. The region outlined by the heavy solid frequency and the mass of the star (see text), 

line shows the combinations allowed by the fre- 
quency and coherence of the QPO (see text). Only 
masses less than M^^^^^ ~ 1.8 and radii less than 
-^max = 16.0 km are allowed. The allowed region col- 
lapses to to the heavy solid horizontal line if the QPO 
frequency is identified as the orbital frequency of gas 
in the innermost stable circular orbit (see § 5.1). The 
steps in the construction are the same for a rotating 
star but the bounds are different (see § 5.2). 

Fig. 14. — Comparison of the mass-radius relations 
for nonrotating neutron stars given by five represen- 
tative equations of state for neutron-star matter with 
the regions of the mass-radius plane allowed for non- 
rotating stars (see Fig. 13 and text) in sources with 
three values of t^Qpo2: the highest observed frequency 
of the higher-frequency QPO in a pair. The light solid 
curves show the mass-radius relations given by equa- 
tions of state A (Pandharipande 1971), FPS (Lorenz, 
Ravenhall, & Pethick 1993), UU (Wiringa, Fiks, & 
Fabrocini 1988), L (Pandharipande & Smith 1975b), 
and M (Pandharipande & Smith 1975a). Each al- 
lowed region is labeled by the value of fQpo2 assumed 
in constructing it. As in Figure 13, the hatched region 
shows the combinations of stellar mass and radius 
that are excluded by j^Qpo2=1220 Hz. The region 
bounded by the heavy line is the region that would 
be allowed for the neutron star in 4U 1636—536 if it 
were not rotating (see text). The star's probable spin 
rate affects the mass-radius relations hardly at all, 
but enlarges the allowed region by ~ 20% (see § 5.2 
and Fig. 15). 

Fig. 15. — Approximate constraints imposed on 

the mass, radius, and equation of state of neutron 
stars by the 1220 Hz QPO frequency observed in 
4U 1636—536, when first-order effects of the stellar 
spin are included. The QPO frequency is assumed 
to be the Keplerian frequency of gas in a prograde 
orbit. The light solid lines are the mass-radius rela- 
tions for the same equations of state as in Figure 14. 
Allowed combinations of mass and radius for the in- 
dicated values of the dimensionless angular momen- 
tum j are those within the corresponding pie-shaped 
regions bounded by the heavy dashed lines. For a 
rotating star, accurate determination of the allowed 
region requires computation of the bounding curves 
with fspin held equal to the observed value rather 
than with j held constant and therefore depends on 
the equation of state assumed as well as the QPO 
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Table 1 

Known High-Frequency QPOs and Burst Oscillations'* 



■'CoiiiploU; as ol 1998 Ala.}- 31. Burst oscillation. 


Source 


Type 


Frequencies (Hz) 


RMS Amplitude (%) 


References 


4U 0614+091 


atoll 


400-600 
500-1145 

327 
630, 727 


6-15 
16.5, 15.8 


Ford et al. 1996, 1997a, 1997b 
van der Klis ct al. 1996c 
Mendez et al. 1997 


4U 1608-52 


atoll 


650-890 
940-1125 
570-800 


5-14 


Berger et al. 1996 
Mendez et al. 1998 
Yu et al. 1997 


4U 1636-536 


atoll 


840-920 
1150-1220 
580* 


6.0 
6.6, 6.1 


Zhang et al. 1996, 1997a 
van der Klis et al. 1996c 
Wijnands et al. 1997b 


4U 1728-34 


atoll 


637-716 

500-1100 
363* 


5.2-6.9 

5.5-8.1 
1.5-5 


Strohmayer et al. 1996a, 1996b, 
1996c, 1997a 


KS 1731-260 


atoll 


524 

900, 1170-1207 


12 

4-5 


Morgan & Smith 1996 

Smith et al. 1997 

Wijnands & van der Klis 1997 


4TI 1735—444 


atoll 


1149 


3.1 


WiinanHs pt fll 1 QQfi 


4U 1820-30 


atoll 


546 796 
1066 


3.2-5.0 


Smale et al. 1996, 1997 


Aql X-1 


atoll 


750-830 

549* 




Zhang et al. 1998a 


Cyg X-2 


Z 


730 -1020 
490-530 


3-5 


Wijnands et al. 1998 


GX 5-1 


Z 


567-895 
325-448 


2.0-6.7 


van der Klis et al. 1996e 


GX 17+2 


z 


470-780 

645-1087 


3-5 


van der Klis et al. 1997a 

Wijnands et al. 1997c 


GX 340+02 


z 


247 -625 
625 -820 


2.5 
2-5 


Jonker et al. 1998 


GX 349+2 


z 


712 
978 


1.2 
1.3 


Zhang, Strohmayer, & Swank 1998b 


Sco X-1 


z 


570-830 
870-1130 


0.9-1.2 
0.6-0.9 


van der Klis et al. 1996a,1996b 
1996d, 1997b 


Unknown 


unknown 


589' 


2-4 


Strohmayer et al. 1996d 
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Table 2 
Location of the Sonic Point 



Stellar radius Strong radiation drag force Weak radiation drag force 



R < i?ms Sonic point at -Rami (Case la) Sonic point at iims (Case lb) 

R > i?nis Sonic point at i?ami (Case 2a) No sonic point in Keplerian 

disk flow (Case 2b) 



Table 3 

Amplitudes of frequencies generated by accretion flow 



Frequency First Order Second Order Third Order 



l^Bs Bi — — 

VKs Ki — 

2t'Bs B2 

^UKs K2 

I'spin jKiBi 

2t'Ks — I'spin 2-^1-^1 

VKs — 2l'spin \K\B2 

t'Ks + i'spin — 

3fKs — 2l'spin jKiB2 

3uks — fgpin 
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Table 4 

Inferred Spin Frequencies of Neutron Stars with Kilohertz QPOs*^ 



^Complete as 


of 1998 May 31. Here Az/ 


is the difference 


betwcxai tfic frequencies of the kilohertz QPOs seen 


simultaneously in the persistent X-ray 


omission of the 


source and i^jurst is the frequency of the brightness 




oscillation seen during type I (thermonuclear) X-ray bursts from the source. 


Kilohertz QPO 


I^spin (Hz) 


'^spin (Hz) 


References 




from A;/ 


from //hiusi 




4U 0614+091 


~330 




Ford et al. 1997a, 1997b 


4U 1608-52 


230-290 




Mendez et al. 1998 


4U 1636-536 


~290 


~580 


Zhang et al. 1996, 1997 


4U 1728-34 


~360 


~363 


Strohmayer et al. 1996a, 1996b, 1996c 


KS 1731-260 


~260 


~520 


Smith et al. 1997, Wijnands & van der Klis 1997 


4U 1735-444 








4U 1820-30 


~275 




Smale et al. 1997 


Aql X-1 




~550 


Zhang et al. 1998a 


Cyg X-2 


~345 




Wijnands et al. 1998 


GX 5-1 


~325 




van der Klis et al. 1996e 


GX 17+2 


~295 




Wijnands et al. 1997c 


GX 340+0 


~325 




Jonker et al. 1998 


GX 349+2 


~266 




Zhang et al. 1998 


Sco X-1 


~250-300 




van der KHs et al. 1996a, 1996b, 1996d, 1997b 


Unknown 




~590 


Strohmayer et al. 1997a 



Table 5 

DiMENSIONLESS ANGULAR MOMENTA FOR t'soin = 363 Hz 



EOS 


Mass {Mq) 


j = cJ/GM^ 


A 


1.40 


0.13 




1.66 


0.10 


FPS 


1.40 


0.16 




1.80 


0.11 


L 


1.40 


0.28 




2.70 


0.17 


M 


1.40 


0.31 




1.80 


0.13 



56 



Atoll Source 




Radius (km) 




Radius (M) 



Radius (km) 




0.09 



— I 1 1 1 1 1 1 1 r 

L b 



0.08 - 



^ 0.07 

O 

> 

^ 0.06 

3 

< 



0.05 - 



0.04, 



"1 1 1 r 



0.02 

0.03 
- - 0.04 
0.05 




J I I I I I I I I I I I I L 



6 7 
Radius (M) 



8 



T 




Radius (M) 



1700 



1600 - 



1500 - 



900 - 



800 



ISCO 




0.01 0.02 0.03 0.04 0.05 0.06 




J I \ ^ I I f I -— f^^ 

6 8 10 12 14 



Optical Depth 



25 



"1 1 1 1 1 1 1 1 1 1 1 1 1 r 



20 - 



4U 1608-52 
#4U 1636-536 
Model C) 



O C) - 



15 

:3 



< 10 



5 - 







J I L 



J I I I I I I I L 







10 20 
Photon Energy (keV) 



30 



stellar Magnetic Fields 

"1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



4U 0614 + 091 



4U 1636-536 



4U 1735-444 



Sco X-1 
I — i — I 

I I I I I I I I I I I I I I L 

,8 1.6 1.4 1.2 

PCA (7-20 keV)/(5-7 keV) 




Radius (km) 




Radius (km) 



3 




